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SUMMARY 


This  report  was  developed  as  a part  of  a program  under- 
taken by  the  Manufacturing  Technology  Directorate  of  Picatinpy 
Arsenal  in  recognition  of  the  need  for  expanded  design  information 
pertaining  to  structures  subjected  to  the  blast  environment  due  to 
accidental  explosions  within  ammunition  facilities.  The  purpose 
herein  is  to  provide  facility  designers  with  procedures,  and  the 
computer  program  written  to  implement  them,  for  determining  the 
gross  motions  of  protective  structures  on  their  supporting  soils. 
The  report  is  intended  to  supplement  the  design  methods  of  the 
tri -service  design  manual  "Structures  to  Resist  the  Effects  of 
Accidental  Explosions"  (TM  5-1300). 

Gross  motions  of  the  structure  subjected  to  high  inten- 
sity blast  loads  are  determined  using  a technique  which  embodies 
the  rigid  body  approximation  to  the  structure's  response  in  con- 
junction with  a discrete  element  representation  of  the  supporting 
soil.  Non-linear  soil  behavior  is  considered  in  the  analysis  and 
incorporated  into  the  computer  program.  The  analysis  is  directed 
primarily  towards  cubicle-type  structures  but  has  application  to 
other  structure  configurations. 

A classification  of  various  soils  is  provided  together 
with  representative  values  of  critical  soil  properties  for' use  in 
the  computer  program. 

The  computer  program  is  presented  in  the  form  of  program 
documentation,  coded  input  card  formats,  input  deck  structures, 
FORTRAN  listing  for  the  COC  6600  computer,  and  sample  problems. 

Detailed  procedures  are  given  for  the  computation  of  the 
time  history  of  the  blast  loads  on  the  structure,  and  the  struc- 
tural design  of  the  foundations.  Numerical  examples  are  included 
to  illustrate  the  procedures.  Conclusions  and  recommendations  are 
also  presented. 
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CONCLUSIONS  AND  RECOMMENDATIONS 


This  report  presents  analysis  techniques  for  predicting 
gross  motions  of  structures,  and  design  criteria  and  procedures 
for  determining  the  size  of  the  foundation  base  slab  necessary 
to  resist  the  soil  pressure  build-up  beneath  the  structure 
during  the  response  of  the  structure  to  the  blast  loads. 

It  Is  recommended  that  procedures  and  the  computer  pro- 
gram presented  in  this  report  be  utilized  in  the  design  of 
blast-resistant  protective  structures  for  facilities  engaged  in 
the  manufacture,  maintenance,  modification,  inspection  and 
storage  of  explosive  materials. 
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SECTION  1 


INTRODUCTION 


1 .1  Background 

The  design  of  facilities  for  the  manufacture,  maintenance, 
modification,  inspection  and  storage  of  explosive  materials  uti- 
lizes special  procedures  and  criteria  in  order  to  avoid  mass  det- 
onations and  explosive  propagation  in  the  event  of  an  accidental 
explosion  and  to  ensure  protection  for  personnel  and  equipment. 

The  basic  design  docuiiient  in  this  area  is  the  tri-service  manual 
"Structures  to  Resist  the  Effects  of  Accidental  Explosions"  (TM  5- 
1300).  This  manual  presents  the  procedures  to  quantitatively 
evaluate  the  ability  of  reinforced  concrete  structures  to  resist 
the  effects  of  a detonation  of  high  explosives.  It  provides  pro- 
cedures and  the  necessary  data  for  the  design  of  laced  reinforced 
concrete  elements  which  are  required  to  resist  close-in  explosions. 

Briefly,  the  methods  in  the  manual  treat  the  design  of 
each  blast-resistant  element  in  a structure  individually.  These 
methods  are  based  on  the  premise  that  the  supports  along  the 
periphery  of  an  element  are  completely  fixed  against  translation 
and  rotation  and  are  capable  of  fully  developing  the  strength  of 
the  element.  Generally,  this  approach  is  adequate  for  the  design 
of  nx)st  protective  structures. 

In  some  design  situations,  however,  an  additional  con- 
sideration, the  motion  of  the  protective  structure  or  barrier  on 
the  supporting  soil,  must  be  included  in  the  design  process. 

This  always  occurs  in  the  design  of  cantilever  wall  and  single 
cell  barriers  which  are  isolated  from  surrounding  structures 
(see  Figures  1 and  2).  These  structures  rely  completely  on  the 
supporting  soil  to  provide  the  needed  resistance  to  the  over- 
turning and  translational  motions.  Some  protective  barriers 
located  within  explosive  storage  and  manufacturing  facilities 
also  fall  into  this  category.  An  example  of  this  is  the  canti- 
lever wall  barrier  of  Figure  3 which  is  placed  in  an  existing 
facility.  In  this  case,  the  existing  foundation  slab  does  not 
have  sufficient  strength  or  rigidity  to  support  the  wall.  There- 
fore, a thick  foundation,  not  integral  with  the  base  slab  of  the 
structure,  is  provided  for  the  wall.  Another  example  is  the 
exterior  blast  wall  of  the  structure  shown  in  Figure  4.  The  entire 
resistance  to  overturning  is  provided  by  the  foundation,  as  the 
side  walls,  shown  in  Figure  4,  are  not  blast-resistant  and  there- 
fore will  fail.  MS  a result,  the  v^all  can  be  expected  to  experi- 
ence large  overturning  motions  under  the  action  of  the  blast. 

The  excessive  height  of  the  wall  will  also  add  to  the  severity  of 
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the  overturning  motions.  In  this  case,  buttress  walls  are  added 
to  restrain  the  exterior  blast  wall.  In  these  situations,  the 
overturning  and/or  translational  motions  of  the  protective  struc- 
tures or  barriers  may  result  in  the  propagation  of  the  explosion 
or  injury  to  personnel  and,  therefore,  these  motions  must  be  con- 
sidered in  the  design. 

In  contrast  to  the  exterior  wall  shown  in  Figure  4,  the 
interior  blast-resistant  dividing  wall  of  the  same  structure  is 
not  susceptible  to  overturning  as  it  is  restrained  by  the  founda- 
tion and  the  side  walls  (in  the  undamaged  portion  of  the  struc- 
ture). In  the  multi-cell  barrier  of  Figure  5,  an  explosion  in 
one  cell  is  confined  by  the  blast-resistant  walls  and  the  overall 
structure  is  restrained  from  overturning  by  the  massive  walls  and 
foundation  slab.  In  these  situations,  the  motion  of  the  protec- 
tive structure  or  barrier  on  the  supporting  soil  is  not  a critical 
factor  and,  therefore,  need  not  be  considered  in  the  design. 

In  this  report,  the  procedures  of  Hi  5-1300  have  been  ex- 
tended to  include  analytical  techniques  to  evaluate  the  motions 
of  a structure  on  its  supporting  soil  and  methods  and  criteria  to 
design  those  elements  which  prevent  the  structure  from  overturning 
and,  if  important,  sliding. 

The  procedures  and  the  computer  program  presented  in  this 
report  were  developed  by  the  Manufacturing  Technology  Directorate 
of  Picatinny  Arsenal,  with  the  assistance  of  Aiimann  & Whitney,  Con- 
sulting Engineers,  as  part  of  the  overall  Picatinny  Arsenal  Safety 
Engineering  Support  Program  for  the  U.S.  At*my  Armament  Command. 

1 .2  Objectives 


The  primary  objective  of  this  report  is  to  present  pro- 
cedures and  the  computer  program  written  to  implement  them  for 
determining  the  gross  motions  of  structures  subject  to  the  blast 
effects  of  high  explosive  detonations.  These  procedures  are  in- 
tended to  supplement  the  design  methods  of  TM  5-1300, 

Secondary  objectives  include: 

1.  The  presentation  of  procedures  for  determining  the 
blast  load  history  on  the  structure. 

2.  The  presentation  of  representative  values  of  critical 
soil  properties  for  use  in  the  computer  program. 

3.  The  presentation  of  criteria  and  p»'Ocedures  for  de- 
signing foundations  of  protective  structures  subjected 
to  large  overturning  and/or  translational  motions. 
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1.3 


Format  of  the  Report 


The  report  is  divided  into  three  main  parts.  The  first 
part,  consisting  of  Sections  2,  3 and  4,  is  devoted  to  a descrip- 
tion of  the  analysis  technique  utili2ed  to  compute  the  response 
of  structures  to  time-dependent  loadings  and  to  the  definition 
of  the  environments  to  which  protective  structures  are  subjected. 
In  Section  2,  the  concept  of  rigid  body  analysis  is  introduced 
together  with  the  equations  of  motion  of  the  structure.  Section 
3 discusses  quantitative  procedures  for  computing  the  blast  out- 
put of  the  explosives,  while  Section  4 presents  the  analytical 
methods  utilized  to  simulate  the  behavior  of  soils  under  the 
dynamic  motions  o^  foundations.  Included  in  the  latter  section 
are  recommended  values  for  critical  soil  properties  to  be  used 
in  the  computer  program. 

The  second  part  consists  of  Sections  5 and  6 and  relates 
to  the  use  of  the  computer  program  written  to  implement  the  pro- 
cedures presented  in  earlier  sections.  Section  5 presents  the 
capabilities  of  the  computer  program  and  a detailed  description 
of  the  input  requirements  including  input  card  formats  and  deck 
structures  for  the  various  options  of  the  program.  Section  6 des- 
cribes in  detail  the  computer  program  printed  output  and  includes 
a discussion  on  the  interpretation  and  utilization  of  the  results 
as  related  to  the  design  of  foundations  for  protective  structures. 

The  third  part  consists  of  the  four  appendices  presented 
at  the  end  of  the  report.  The  first  three  appendices  present 
various  quantitative  procedures  utilized  in  the  performance  of 
dynamic  analyses  of  protective  structures  along  with  numerical 
examples  which  illustrate  the  use  of  these  procedures.  Appendix 
A contains  the  p»'ocedures  for  deteraining  tiie  impulse  loads  on 
foundation  pads  of  cantilever  wall  barriers,  resulting  from  the 
detonation  of  an  explosive  charge  located  outside  the  periphery 
of  the  structure.  Appendix  B presents  the  procedure  for  computing 
the  arrival  time  and  duration  of  blast,  pressures  on  a protective 
structure  subjected  to  a close-in  explosion.  Appendix  C contains 
criteria  and  ptxjcedures  for  designing  the  foundations  and  other 
elements  which  stabilize  the  structure  against  overturning. 
Appendix  D contains  the  FORTRAN  listing  of  the  computer  program 
for  the  CDC  6600  computer.  Included  in  this  appendix  are  samples 
of  the  punched  card  input  data  decks  and  the  printed  output  of 
the  program. 

In  order  to  simplify  the  overall  presentation,  most  of 
the  directly  applicable  material  fvw  TH  5-1300  has  not  been 
repeated  in  this  report.  As  far  as  possible,  applicable  equa- 
tions, design  charts  and  tables  and  commentary  material  have  beer, 
included  herein  by  refet'ence. 
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SECTIGN  2 


METHOD  OF  ANALYSIS 

2.1  Introduction 


The  method  of  analysis  fo.'  computing  the  gross  motions  of 
a protective  structure  subjected  to  the  blast  effects  of  a high 
explosive  detonation  is  presented  in  this  section.  The  response 
of  the  structure  to  the  action  of  the  blast  is  determined  by  per- 
forming a time-history  dynamic  analysis.  The  analysis  considers 
the  effect  of  the  supporting  soil  on  the  response  of  the  structure. 
A rigorous  approach  to  this  problem  is  exceedingly  difficult  and 
time  consuming  and,  in  fact,  is  not  warranted  since  only  the  gross 
motions  of  the  stru«-ti.re  are  required.  Therefore,  the  problem  is 
reduced  to  the  simplest  form  possible  but  still  retains  the  non- 
linear complexity  of  the  soil.  This  is  accomplished  utilizing 
the  rigid  body  approxi'  ^tion.  Briefly,  this  method  treats  the 
structure  as  a rigid  body  consisting  of  infinitely  stiff  elements. 
The  blast  loads  and  element  inertial  loads  are  assumed  to  act 
through  the  center  of  gravity  of  the  structure.  This  reduces  the 
response  of  the  structure  to  the  motion  of  the  center  of  gravity 
and,  therefore,  greatly  simplifies  the  computation. 

The  rigid  body  approach  is  well  suited  to  the  problem  at 
hand.  Generally,  protective  structures  consist  of  massive  laced 
concrete  wall  elements  supported  by  equally  massive  foundation 
slabs  resting  on  relatively  soft  soils.  Consequently,  the  indi- 
vidual elements  will  usually  reach  their  respective  peak  res- 
ponses well  before  any  significant  motion  of  the  foundation  has 
occurred.  As  a result,  there  is  usually  little  interaction  be- 
tween the  individual  responses  of  the  elements  and  the  gross 
response  of  the  overall  structure.  Therefore,  the  rigid  body 
analysis  will  generally  produce  reliable  estimates  of  the  gross 
motions  of  the  structure  without  a severe  overestimate  of  the 
response  of  the  individual  elements.  The  section  that  follows 
provides  a more  detailed  discussion  of  the  mechanics  of  the 
method  and  presents  the  equations  of  motion  of  the  structure. 

2.2  Equations  of  Motion 

The  structure  is  considered  to  be  a rigid  body  constrained 
to  move  parallel  to  the  x-y  plane  (see  Figure  6).  This  condition 
is  easily  attained  if  we  consider  the  x-y  plane  to  be  a plane  of 
symmetry  of  the  structure  and  the  loading  on  the  structure  to  be 
symmetrical  about  this  plane.  The  mass  and  rotary  inertia  of  the 
entire  structure  are  lumped  at  the  center  of  gravity.  Blast  pres- 
sures acting  on  the  individual  elements  as  well  as  the  soil  pres- 
sures acting  on  the  foundation  are  transposed  to  the  center  of 
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gravity.  In  this  way,  the  response  of  the  structure  is  reduced 
to  three  degrees  of  freedom.  These  are  u,  v and  o which  are, 
respectively,  the  translations  in  the  x and  y axes  and  the  rota- 
tion about  the  z axis. 

The  soil  is  represented  analytically  as  a series  of  hori- 
zontal and  vertical  elements  which  resist  the  motion  of  the  struc- 
ture. The  soil  elements  exhibit  non-linear  behavior.  The  total 
resistance  developed  in  the  soil  at  any  timeduring  the  response 
is  dependent  on  the  displaced  configuration  of  the  foundation. 

A detailed  discussion  of  the  soil  model  is  presented  in  Section  4. 

The  equations  of  motion  for  this  system  are  listed  below: 
mu  + Rfj  = H(t)  (2.1) 

mv  + Ry  = V(t)  + mg  (2.2) 

10  + Rq  = M(t)  (2.3) 

wherein: 

u,  V,  0 = horizontal,  vertical  and  rotational 
accelerations  of  the  structure 

m = mass  of  structure 

I = mass  moment  of  inertia  of  the  structure 
about  the  z axis  at  the  center  of  gravity 

g = accele'^ation  of  gravity 

H(t)  = resultant  of  horizontal  blast  loads  at 
time  t 

V(t)  = resultant  of  vertical  blast  loads  at 
time  t 

H(t)  = moment  of  resultant  blast  loads  about 
the  z axis  at  the  center  of  gravity 
at  time  t 

R|^  = horizontal  resistance  of  soil 

Ry  = vertical  resistance  of  soil 

Rq  =•  moment  of  horizontal  and  vertical  soil 
resistance  forces  about  the  center  of 
gravity. 


The  solutions  for  u,  v and  6 are  readily  obtained  by 
nunerically  integrating  the  equations  of  motion  utilizing  the 
constant  velocity  procedure.  This  is  a procedure  by  which  the 
differential  equations  of  motion  are  solved  step  by  step,  start- 
ing at  zero  time,  when  the  displacement  and  velocity  are  presum- 
ably known.  The  time  scale  is  divided  into  discrete  intervals, 
and  one  progresses  by  successively  extrapolating  the  displacement 
from  one  time  station  to  the  next.  This  procedure  is  well  suited 
to  the  solution  of  non-linear  problems  since  it  allows  for  the 
inclusion  of  multi-linear  resistance  functions. 
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SECTION  3 
BLAST  LOADING 


3.1  Introduction 


This  section  presents  the  procedures  for  computing  the 
load  history  for  protective  structures  subjected  to  close-in  ex- 
plosions. The  elements  of  the  structure  experience  relatively 
high  initial  pressures  which  decay  rapidly  to  zero.  Definition 
of  the  loads  for  use  in  design  of  the  individual  elements  is 
usually  limited  to  the  computation  of  the  average  impulse  as  the 
elements  attain  their  peak  responses  long  after  the  load  acting 
on  them  has  decayed  to  zero.  However,  the  analysis  method  of 
Section  2,  which  employs  numerical  integration  techniques,  re- 
quires a complete  definition  of  the  load  history.  The  load  his- 
tory is  characterized  by  a peak  pressure,  an  idealized  pressure 
function  and  a duration.  The  quantities  required  to  determine 
the  load  history  are  the  average  impulse,  the  time  of  arrival  of 
the  shock  front  and  the  duration  of  the  blast  pressures  for  each 
element  of  the  structure.  The  idealized  pressure  function  uti- 
lized is  an  initially  peaked  triangle  defined  by  a peak  pressure 
and  a duration.  Extensive  studies  of  this  simplified  loading 
function  have  shown  it  to  yield  relatively  accurate  results  in 
most  applications.  The  peak  pressures  are  computed  using  the 
aforementioned  data.  The  total  force  on  each  surface  is  then 
calculated  as  a function  of  time  and  transposed  to  the  center 
of  gravity  of  the  structure.  Figure  7 illustrates  the  manner  in 
which  the  blast  load  history  is  calculated.  A more  detailed  des- 
cription of  the  load  history  computation  is  provided  in  Section 
5.3.3. 


The  procedures  for  computing  the  data  required  to  deter- 
mine the  load  h>tory  on  the  structure  are  discussed  in  the  sec- 
tions that  follow. 

3.2  Computation  of  Average  Impulse 

The  computation  of  the  impulse  loads  on  the  structure  is 
outlined  in  Reference  1.  This  manual  provides  a systematic  pro- 
cedure to  determine  the  impulse  loads.  The  method  is  extremely 
tedious  and  requires  extensive  extrapolation  and  interpolation  of 
the  data.  To  facilitate  the  computation,  a computer  program  was 
developed  (Reference  2)  which  uses  the  quantity  and  location  of 
the  explosive  and  the  geometry  of  the  structure  to  determine  the 
impulse  load.  Like  the  impulse  data  given  in  TM  5-1300,  the 
"Impulse  Program"  can  be  used  for  computing  the  impulse  loads  on 
the  walls  and  foundations  of  protective  structures. 
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The  impulse  data  of  Reference  1 is  given  in  the  form  of 
charts  which  express  the  average  scaled  impulse  (impulse  per  cube 
root  of  charge  weight)  as  a function  of  the  parameters  which  de- 
fine the  structure  configuration,  charge  size  and  charge  location. 
The  charts  contain  the  impulse  loads  for  cantilever  wall  barriers 
and  two-,  three-  and  four-wall  cubicles.  The  impulse  charts  and 
the  "Impulse  Program"  were  developed  to  compute  the  loads  when 
the  charge  is  located  within  the  periphery  of  the  cell.^  The  case 
of  a charge  located  outside  of  a cell  was  not  included  in  the 
charts  of  Reference  1 nor  considered  in  the  development  of  the 
"Impulse  Program".  This  situation  rarely  occurs  in  the  design  of 
cells,  but  it  is  not  uncoirmon  in  the  design  of  cantilever  wall 
barriers.  In  this  situation,  the  charts  of  Reference  1 and  the 
"Impulse  Program"  can  be  utilized  to  determine  the  impulse  on  the 
cantilever  wall  but  additional  data  is  required  to  determine  the 
impulse  on  the  foundation. 

In  order  to  provide  a means  for  determining  the  impulse 
loads  on  the  foundations  of  cantilever  wall  barriers,  additional 
impulse  charts  were  prepared  with  a modified  version  of  the  "Im- 
pulse Program".  The  charts  express  the  average  scaled  impulse 
as  a function  of  the  dimensions  of  the  structure,  charge  weight 
and  charge  location.  These  additional  impulse  charts  are  pre- 
sented in  Appendix  A. 

Use  of  the  impulse  charts  may  require  interpolation  of  the 
data  in  many  cases.  A procedure  for  interpolating  the  chart  pa- 
rameters is  presented  in  Appendix  A.  The  procedure  is  based  on 
the  methods  presented  in  Section  4-10  of  Reference  1.  Included 
in  the  appendix  is  a sample  problem  which  illustrates  the  use  of 
the  procedure. 

3.3  Computation  of  Arrival  Times  and  Load  Duration 

A large  body  of  empirical  data  has  been  collected  on  the 
time  of  arrival  of  the  blast  front  and  the  duration  of  the  posi- 
tive phase  of  the  pressure  wave.  It  has  been  determined  that 
most  blast  parameters  are  scaled  by  the  cube  root  of  the  charge 
weight  for  a given  explosive.  These  parameters  are  thus  pre- 
sented in  a scaled  form,  i.e.,  t/\/W'/3^  is  the  time  of  ar- 
rival) and  t(j/W^'^  (tg  is  the  duration  of  the  positive  phase) 
and  plotted  as  functions  of  the  scaled  distance,  Z = 
where  is  the  radial  distance  from  the  charge.  Such  correlation 
of  data  is  found  in  Figure  4.5  of  Reference  1 and  reproduced  in 
Figure  B.l  of  this  report.  To  detemine  the  arrival  time  or 
duration  at  a point  of  interest,  one  enters  the  curves  of  Figure 
B.l  with  the  given  scaled  distance  and  reads  the  parameter  of 
interest  from  the  appropriate  curve. 
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The  procedure  for  computing  the  arrival  time  of  the  blast 
wave  and  the  duration  of  the  load  on  those  elements  of  the  struc- 
ture directly  exposed  to  the  blast  is  based  on  the  methods  (Sec- 
tion 4)  and  empirical  data  (Figure  4.5)  presented  in  Reference  1. 
The  procedure  begins  with  the  definition  of  the  arrival  time  of 
the  blast  wave  on  an  element  ("t«"  in  Figure  7)  as  the  time  re- 
quired for  the  wave  to  arrive  at  the  point  on  the  element  nearest 
to  the  explosive.  An  estimate  of  the  load  duration  ("to"  in  Fig- 
ure 7)  is  then  obtained  by  computing  an  average  time  for  the  wave 
to  fully  engulf  the  element  and  adding  this  quantity  to  the  aver- 
age of  the  load  durations  at  those  points  on  the  element  farthest 
from  the  explosive  (see  Figure  B.2).  The  times  of  arrival  and 
load  durations  for  the  points  of  interest  on  the  element  are  ob- 
tained from  Figure  B.l. 

Appendix  B presents  an  outline  of  the  method  and  a sample 
application  of  it. 
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SECTION  4 


SOIL  - STRUCTURE  INTERACTION 

4.1  Introduction 

This  section  presents  the  soil -structure  interaction 
model  for  simulating  the  non-linear  behavior  of  soils  subjected 
to  the  dynamic  motions  of  foundations.  Included  in  the  discus- 
sion is  the  relationship  between  the  principal  features  of  the 
model  and  the  actual  behavior  of  the  soil.  Following  this  is  a 
system  for  classification  of  various  soils  and  a tabulation  of 
their  respective  properties  for  use  in  the  interaction  model. 

The  model  is  used  in  the  dynamic  analysis  of  protective 
structures.  In  the  analysis,  the  model  utilizes  the  motions  of 
the  foundation  to  determine  the  resisting  forces  in  the  soil. 
These  forces  are  then  substituted  into  the  equations  of  motion 
presented  in  Section  2. 

4.2  Soil  Structure  Interaction  Model 

The  significant  physical  characteristics  of  the  soil 
medium  that  are  incorporated  into  the  interaction  model  are: 

1.  The  effect  of  a continuous  supporting  medium 
beneath  the  foundation. 

2.  The  resistance  to  the  downward  and  horizontal 
motions  of  the  foundation.  Both  linear  and 
bilinear  resistance  functions  a»'e  included 
in  the  model. 

3.  The  lack  of  rebound  experienced  when  the  founda- 
tion moves  upward. 

4.  The  effect  of  friction  between  the  foundation 
and  the  soil . 

The  effect  of  the  continuous  medium  is  simulated  by  rep- 
resenting the  soil  as  a series  of  discrete  element  pairs  attached 
to  the  foundation  at  equally  spaced  intervals.  Each  pair  con- 
sists of  a horizontal  and  vertical  element  as  shown  in  Figure  8. 
Tlie  resistance  developed  in  each  elerent  at  any  time  during  the 
response  is  dependent  upon  the  displacesiient  of  the  foundation  at 
the  attachment  point  of  the  element.  Generally,  10  to  15  pairs 
of  elanents  are  necessary  to  produce  an  adequate  representation 
of  a continuous  medium. 
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Resistance  to  the  downward  motion  of  the  foundation  is 
developed  in  the  soil  by  the  elastic  deformation  of  the  grains 
and  the  inter-granular  friction  generated  as  the  grains  slip  past 
each  other  and  fill  the  voids.  Resistance  to  horizontal  transla- 
tion of  the  foundation  results  from  the  elastic  deformations  of 
the  grains  until  the  structure  begins  to  slide,  after  which  the 
motion  is  resisted  by  dynamic  friction.  A more  detailed  dis- 
cussion of  friction  and  sliding  of  the  foundation  is  presented 
in  a subsequent  paragraph. 

Conventional  analysis  methods  utilize  equivalent  springs 
to  simulate  the  soil  resistance.  Applications  of  this  method  are 
presented  in  References  3 and  4.  Equations  for  determining  equiva- 
lent spring  constants  for  the  soil  are  presented  by  Whitman  in 
Reference  5.  These  equations  were  derived  utilizing  elastic  half 
space  theory.  The  equations  for  the  vertical  and  horizontal 
spring  constants  are: 

Vertical:  Ky  = - u)  (4.1) 

Horizontal:  Kj.  = 2(1  + u)Gs  ^ (4.2) 

wherein: 

Ky,  = total  spring  constant  for  vertical 
and  horizontal  translation 

G = shear  modulus  for  soil 

u = Poisson's  ratio  for  soil 

B = length  of  foundation,  alonn  axis  of 
rotation  for  rocking  or  normal  to 
direction  of  horizontal  force 

L = length  of  rectangular  foundation  in 
plan  of  rotation  for  rocking  or  in 
direction  of  horizontal  force 

6j.  = Influence  coefficient  for  horizontal 
spring  constant 

8,  = Influence  coefficient  for  vertical 
spring  constant 

Using  the  sprino  constants  given  by  Equation  4.1  and  4.2,  the 
elastic  constants  of  the  individual  soil  eleiments  are  determined 
by: 
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ky  = Ky/NS 
k,  = yNS 


(4.3) 

(4.4) 


wherein 


NS  = number  of  soil  elements  used  in 
analysis. 

The  elastic  resisting  forces  are: 

NS 

Ru  = I k„(u  - h'sin  o)  (4.5) 
1=1 

NS 

Ry  = r ky(v  + x^sin  e) 


NS 

+ T^ky(Vi)sx 

NS 

R = r kw(v  + XiSin  o)(x^ros  e - h'sin  o) 
1=1  ^ 

NS 

- T.  kv(u  - h'sin  o)(h'cos  o + x<sin  o) 
i=l 

(4.7) 

The  variables  in  these  expressions  are  defined  below  and  illus- 
trated in  Figure  6. 

u,  V,  8 = translation  of  structure  in  the 

X and  y axes  and  the  rotation  about 
the  z axis 

h'  = vertical  distance  from  center  of 
gravity  to  soil -structure  interface 

X,-  = horizontal  distance  from  center  of 
gravity  to  soil  element  attachment 
point  on  foundation 
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(V-j)sT  = static  deflection  under  weight  of 

structure  at  soil  element  attachment 
point 

(u  - h'sin  o)  = horizontal  displacement  of  the 
foundation 

(v  + x^sin  o)  = vertical  displacement  of  soil 
element  attachment  point. 

The  relationships  presented  in  the  preceding  paragraph 
apply  when  the  soil  is  being  loaded.  When  the  foundation  mo\,^s 
upward  away  from  the  soil,  only  the  energy  expended  in  elastically 
deforming  the  grains  is  recovered.  The  energy  recovered  is  a 
strall  portion  of  the  total  energy  expended  in  compressing  the 
soil,  with  the  result  that,  in  effect,  the  soil  rebound  is  cut 
off.  This  lack  of  rebound  is  exhibited  in  both  granular  and  clay 
soils.  Clay  soils  exhibit  more  elastic  behavior  but  the  major 
portion  of  the  deformation  is  still  due  to  slip  between  the  par- 
ticles. The  non-rebounding  behavior  of  the  soil  is  simulated 
analytically  by  simply  disconnecting  the  soil  element  when  the 
attachment  point  on  the  foundation  attains  its  peak  dov'/nward 
displacement. 


The  preceding  two  paragraphs  have  presented  a discussion 
of  the  loading  and  unloading  features  of  the  interaction  model. 

In  sunmary,  the  model  behaves  elastically  when  loaded  and  ex- 
hibits no  rebound  when  unloaded.  Figure  9 depicts  this  behavior. 
The  upper  portion  of  the  figure  shows  the  loading  cycle  to  be 
linear.  In  the  lower  portion  of  the  figure,  a bilinear  loading 
cycle  is  depicted.  This  latter  curve  is  a more  accurate  repre- 
sentation of  the  actual  behavior  of  the  soil.  As  the  loading 
increases,  the  voids  are  filled  and  a greater  portion  of  the  re- 
sistance results  futim  the  elastic  defonnation  of  the  grains 
thereby  increasing  the  stiffness  of  the  soil.  Generally,  the 
linear  nx)del  is  adequate  for  the  solution  at  hand,  but  when  the 
appropriate  data  is  available,  the  bilinear  model  should  be 
utilised  in  the  dynamic  analysis. 


As  noted  previously,  the  resistance  to  the  hcrisontal 
motion  of  the  foundation  is  developed  by  the  elastic  deforration 
of  the  grai.ns  until  sliding  cocraiences.  This  resistance  can  be 
simulated  using  a linear  function.  At  some  point,  the  resistajvce 
exceeds  the  friction  force  on  the  foundation  and  sliding  cotnuences 
and  continues  until  the  foundation  attains  its  peak  horizontal  dis- 
placement. At  this  point,  the  foundation  regains  contact  with  the 
soil  and  proceeds  to  displace  in  the  opposite  direction.  During 
this  stage,  as  in  the  initial  stage  of  the  motion,  th;£  resistance 
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Figure  10.  Horizontal  resistance  vs.  deflection. 
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Figure  H.  Horizontal  displacement  history. 


24 


is  developed  by  elastically  deforming  the  grains.  This  effect  is 
illustrated  in  Figure  10  for  a constant  normal  force.  From 
Point  A to  Point  B,  the  resistance  is  linear  until  the  friction 
force  is  exceeded  at  Point  B.  Sliding  coimiences  at  Point  B and 
continues  until  the  foundation  reaches  its  peak  horizontal  dis- 
placement at  Point  C,  at  which  time  the  structure  regains  contact 
with  the  soil.  From  Point  C to  Point  D,  the  resistance  is  again 
linear.  After  Point  n,  sliding  in  the  opposite  direction  occurs 
between  the  slab  and  soil.  After  Point  E,  elastic  soil  forces  are 
redeveloped  until  the  structure  comes  to  rest  at  Point  F where  a 
permanent  displacement  will  be  realized.  Figure  11  depicts  the 
time  history  of  the  horizontal  motion  of  the  foundation  and  iden- 
tifies the  key  points  when  the  effect  of  friction  enters  into  the 
horizontal  resistance. 

4.3  Soil  Properties 

Definitive  determination  of  representative  values  for 
critical  soil  properties  for  use  in  a design  study  is  highly 
problematic.  Ideally,  soil  samples  should  be  collected  at  the 
construction  site  and  tested  unde“  conditions  representative  of 
the  anticipated  operating  movement.  Soil  spring  constants  would 
then  be  evaluated  from  the  test  results  taking  ,nto  account: 

1.  The  effect  of  partial  embedment  of  the  footing. 

2.  The  dependence  of  the  spring  constants  upon  the 
initial  static  stress  as  well  as  upon  the  magnitude 
of  the  dynamic  stress  increment. 

3.  The  distribution  of  stresses  over  the  contact  area 
between  the  foundation  and  the  soil,  and 

4.  The  dependence  of  the  spring  constants  upon  the 
size  of  the  contact  area  which,  in  turn,  depends 
upon  the  variation  of  the  soil  molecules  with 
depth  or  with  the  presence  of  a layered  soil 
structure. 

Under  more  realistic  circumstances,  it  is  likely  that 
the  soil  data  available  to  the  designer  are  the  result  of  a mini- 
mum of  shallow  test  borings  with  visual  descriptions  of  the  soils 
encountered  accompanying  values  of  "N“  as  determined  from  standard 
penetration  lests  (N  is  the  number  of  blows  of  a 140-pound  haniner 
dropping  30  inches  required  to  drive  a sampler  of  2-inch  outer 
diameter  1 foot  into  the  soil).  It  is  necessary  to  correlate  this 
type  of  Information  with  the  modulus  of  elasticity  required  for 
analysis. 
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Tables  1 and  2 summarize  a correlation  of  modulus  of 
elasticity  with  soil  type  and  "N"  value.  The  tabulated  moduli 
were  extracted  from  Reference  6.  The  values  presented  were  de- 
rived from  repeated  load  tests  and  modified  utilizing  observed 
resonant  frequencies  in  the  material.  These  values  are  approxi- 
mate and  should  only  be  used  where  the  time  or  resources  for  the 
determination  of  the  soil  parameters  by  a more  rigorous  method 
are  not  available. 

With  the  data  correlation  of  Tables  1 and  2,  the  designer, 
at  best,  has  only  a gross  estimate  of  the  actual  soil  properties. 
The  soil  properties,  especially  the  modulus  of  elasticity,  affect 
both  the  stability  of  the  structure  and  the  strength  requirements 
of  the  foundation.  A gross  overestimate  of  the  soil  modulus  could 
result  in  a structure  that  overturns  or  experiences  large  hori- 
zontal motions  under  the  action  of  the  blast  while  a gross  under- 
estimate of  the  soil  modulus  could  result  in  a failure  of  the 
foundation.  Consequently,  in  the  absence  of  more  reliable  data 
(such  as  described  in  the  initial  paragraph  of  this  section),  the 
structure  should  be  designed  for  the  range  of  soil  properties 
(specified  in  Tables  1 and  2)  for  the  particular  type  of  soil  at 
the  construction  site. 
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SECTION  5 


COMPUTER  PROGRAM 


5.1  Introduction 


This  section  presents  the  computer  program  used  to  Imple- 
ment the  procedures  (for  performing  dynamic  analyses  of  protective 
structures)  developed  in  Sections  3,  4 and  5 of  this  report. 

The  program  is  essentially  a special  purpose  pv'ogram,  de- 
signed specifically  for  performing  dynamic  analyses  of  structures 
subjected  to  close-in  explosions.  To  this  end,  the  bulk  of  the 
input  data  required  for  the  dynamic  analysis  is  computed  inter- 
nally by  the  program.  Items  such  as  the  inertial  properties  of 
the  structure  (weight,  moment  of  inertia,  etc.),  the  time  history 
of  the  blast  loads  acting  on  the  structure  and  the  elastic  con- 
stants and  spacing  of  the  soil  elements  beneath  the  foundation 
are  computed  by  the  program.  These  computations  are  performed 
for  the  most  common  types  of  cantilever  wall  and  multi -wall  bar- 
riers such  as  those  shown  in  Figure  12. 

The  inclusion  of  these  computations  in  the  program  greatly 
diminishes  the  quantity  of  input  data  entered  on  punched  cards. 

In  addition,  the  program  input  is  tailored  to  facilitate  the  data 
preparation  for  the  most  common  types  of  protective  barriers  such 
as  those  shown  in  Figure  12.  The  range  of  applicability  of  the 
program  can  be  extended,  however,  by  overriding  the  built-in  com- 
putational routines.  In  this  way,  the  program  can  be  used  to 
analyze  structures  of  any  configuration. 

5.2  Program  Capabilities 

The  computer  program  perfo'^ns  a rioid  body  analysis  of  a 
protective  structure  subjected  to  che  blast  pressures  of  an  ex- 
plosive detonation.  The  results  of  the  computeri r.ed  analysis 
consist  of  the  displacetn'^nt,  velocity  and  acceleration  time  his- 
tories of  the  structure  The  bearing  pressures  beneath  the 
foundation  are  also  determined.  Options  are  included  for  the 
computation  of  the  shears  and  moments  for  the  foundation  slab 
design  of  cantilever  wall  barriers  and  the  peak  response  time  for 
the  back  wall  element  (see  Figure  14)  designed  to  the  incipient 
failure  or  post-failure  fragment  conditions. 

The  computer  program  contains  three  optional  modes  of 
operation.  They  are: 
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1.  Normal  Option:  The  "Normal  Option"  mode  is  used  for 
analyzing  the  most  cotmion  types  of  protective  struc- 
tures encountered  in  explosive  manufacturing  and 
storage  facilities.  It  contains  internal  routines 
which  calculate  the  inertial  properties  of  the  struc- 
ture and  the  time  history  of  the  blast  loads  acting 
on  the  structure.  To  compute  the  load  history,  the 
explosive  data  for  TNT  were  incorporated  into  the 
program.  The  data  may  be  utilized  for  other  ex- 
plosives by  applying  a TNT  equivalence  for  the  par- 
ticular explosive  to  the  charge  weight  used. 

2.  Special  Loading  Options  The  "Special  Loading  Option" 
is  used  to  accommodate  rectangular  or  trapezoidal 
load  histories  in  the  analysis. 

3.  General  Structure  Option;  The  "General  Structure 
Option'*  is  used  to  extrnd  the  applicability  of  the 
computer  program  to  structures  of  any  arbitrary 
configuration. 

5.3  Input  Requirements 

5.3.1  General  Input 

The  input  to  the  program  consists  of  punched  cards  which 
contain  the  following  data: 

1.  The  configuration  and  dimensions  of  the  structure. 

2.  The  properties  of  the  soil  for  use  in  the  inter- 
action model  described  in  Section  4. 

3.  The  quantities,  locations  and  blast  outputs  (average 
impulses  on  the  elements  of  the  structure)  of  the 
explosives. 

Since  the  computer  program  has  been  designed  specifically 
for  the  most  common  types  of  protective  structures,  the  quantity 
of  the  input  containing  the  configuration  and  the  dimensions  of 
the  structure  and  the  blast  output  of  the  explosives  has  been 
greatly  diminished.  Computational  flexibility  is  retained, 
however,  by  including  provisions  for  overriding  those  features 
of  the  program  that  apply  to  standard  configurations  of  protec- 
tive structures.  In  this  way,  data  applicable  to  any  type  of 
structure  can  be  entered.  However,  the  quantity  of  data  re- 
quired for  this  purpose  is  much  greater  than  what  is  required 
for  normal  program  applications. 
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As  noted  previously,  a computational  routine  has  been  in- 
stalled in  the  program  for  the  purpose  of  calculating  the  peak 
response  time  of  the  back  wall  element.  To  utilize  this  routine, 
additional  input  describing  the  structural  details  of  the  element 
is  required.  Utilization  of  this  routine  as  well  as  the  input 
required  for  it  is  optional. 

A description  of  the  various  input  quantities  required  is 
contained  in  the  following  section. 

5.3.2  Configuration  of  the  Structure 

As  applied  to  the  most  common  types  of  protective  struc- 
tures, the  structural  configuration  is  described  in  terms  of  the 
number  and  size  of  the  blast-resistant  elements  making  up  the 
structure.  The  input  required  for  the  "Normal  Option"  mode  of  the 
program  is  designed  specifically  for  the  protective  structure  con- 
figurations shown  in  Figure  12.  Included  in  the  figure  are  the 
dimensions  and  parameters  necessary  to  completely  define  the  struc- 
tural geometry.  The  program  uses  the  dimensions  illustrated  in 
Figure  12  to  compute  the  following  quantities  for  use  in  the 
dynamic  analysis  of  the  structure: 

1.  Weight  and  moment  of  inertia  of  the  structure. 

The  computation  is  restricted  to  structures 
constructed  of  nom«l  weight  concrete  weighing 
150  pounds  per  cubic  foot. 

2.  Location  of  the  center  of  gravity  of  the  struc- 
ture as  shown  in  Figure  14. 

3.  The  areas  of  the  surfaces  on  the  structure 
subjected  to  blast  pressures.  The  computation 
is  limited  to  surfaces  normal  to  the  plane  of 
motion  of  the  structure. 

4.  The  location  of  the  centroids  of  the  loaded 
surfaces  relative  to  the  center  of  gravity 
of  the  structure.  Only  the  centroids  of  the 
surfaces  normal  to  the  plane  of  motion  of  the 
structure  are  computed. 

5.  The  plan  area  of  the  foundation. 

In  applications  of  the  program  to  non-standard  configura- 
tions, the  quantities  listed  above,  with  the  exception  of  Item  5, 
have  to  be  computed  beforehand  and  entered  using  the  optional  in- 
put formats  provided  for  this  purpose.  This  optional  form  of 
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entering  the  data  is  the  "General  Structure  Option".  The  pertinent 
data  relating  to  the  dimensions  of  the  foundation  is  still  required 
and  must  be  entered  using  the  regular  input  card  for  the  structural 
geometry. 

5.3.3  Properties  of  the  Soil 

The  soil  properties  are  utilized  with  the  soil  structure 
interaction  model  described  in  Section  4.  The  quantities  required 
are  the  modulus  of  elasticity,  Poisson's  ratio  and  the  appropriate 
friction  constant  for  the  soil.  Soil  properties  for  various  soils 
are  provided  in  Section  4.3.  Shape  factors  for  the  foundation  (as 
defined  in  Reference  5)  are  also  required.  This  data  is  utilized 
by  the  program  to  compute  the  elastic  spring  constants  for  the 
loading  cycle  of  the  interaction  model. 

The  inclusion  of  non-linear  behavior  is  accomplished  by 
entering  two  different  moduli  of  elasticity  and  the  normal  stress 
at  which  the  modulus  changes.  For  bilinear  soil  behavior,  the 
interaction  model  utilizes  a bilinear  loading  cycle  in  the  analysis 
with  spring  constants  being  computed  for  both  moduli  of  elasticity. 

5.3.4  Blast  Output  of  Explosives 

In  applications  of  the  program  to  standard  structural  con- 
figurations, the  blast  output  of  the  explosives  is  described  in 
terms  of  charge  weights  and  locations  and  the  average  unit  impulses 
on  all  loaded  surfaces.  This  limited  description  of  the  explosive 
blast  output  is  all  that  is  required  when  utilizing  the  "Normal 
Option"  mode  of  the  program.  The  program  uses  this  data  to  compute 
the  load  history  on  the  structure.  The  computation  proceeds  in  two 
stages. 


In  the  first  stage,  the  following  quantities  are  computed 
utilizing  the  procedure  presented  in  Section  3.3,  the  weights  and 
locations  of  the  explosive  charges  entered  on  the  input  cards,  and 
the  TNT  explosives  data  contained  in  the  prograf;: 

1.  The  smallest  load  duration  on  any  surface. 

2.  The  arrival  time  of  tlie  blast  wave  on  every  surface. 

3.  The  arrival  time  of  the  blast  plus  the  average 
duration  of  the  loading  on  every  surface. 

The  above  quantities  are  computed  for  each  explosive  charge 
and  stored  internally  in  the  computer. 
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In  the  second  stage,  the  total  blast  load  history  on  the 
structure  is  computed.  The  loadings  on  the  structure  produced  by 
each  explosive  charge  are  computed  individually.  In  the  computa- 
tion cycle  for  each  explosive  charge,  the  load  on  each  surface  is 
computed  in  the  following  manner.  First,  the  program  utilizes  the 
impulse  data  entered  on  the  input  cards  and  the  quantities  (Items 
2 and  3)  listed  on  the  preceding  page  to  determine  the  peak  pres- 
sure acting  on  the  surface.  The  pressure  is  computed  using  a 
triangular  (linear)  pressure-time  history.  The  peak  load  on  the 
surface  is  then  computed  by  multiplying  the  peak  pressure  by  the 
area  of  the  surface.  After  this,  the  load  history  is  determined 
by  dividing  the  triangular  load-time  function  into  a series  of 
small  equal  time  increments.  In  order  to  insure  that  the  total 
impulse  is  considered  in  the  analysis,  the  program  digitizes  the 
load  history  at  a time  increment  equal  to  1/20  of  the  shortest 
load  duration  on  any  one  surface  of  the  structure.  As  the  load- 
time history  is  digitized,  the  time  history  of  the  moment  of  the 
load  around  the  center  of  gravity  of  the  structure  is  also  com- 
puted. The  load-  and  moment-time  histories  are  then  added  to  the 
total  force  vector  for  the  structure. 

Application  of  the  program  to  non-standard  protective 
structures  requires  that  the  data  listed  on  the  preceding  page 
(Items  1,  2 and  3)  be  computed  and  entered  on  the  alternate  input 
formats  provided.  This  alternate  form  of  entering  the  data  is 
the  "Special  Loading  Option"  which  can  also  be  utilized  to  accom- 
modate rectangular  and  trapezoidal  pressure-time  functions  in  the 
analysis. 


5.3.5  Structural  Details  of  the  Back  Wall  Element 

This  input  consisting  of  structural  design  data  relating 
to  laced  reinforced  concrete  walls  is  entered  for  the  purpose  of 
computing  the  maximum  response  time  of  the  back  wall  element. 

In  this  report,  the  back  wall  element  is  defined  as  the  principal 
blast-resistant  wall  perpendicular  to  the  x-y  plane  of  motion  of 
the  structure  as  shown  in  Figure  14.  The  response  time  computa- 
tion is  performed  using  the  procedures  outlined  in  Chapters  5 and 
6 of  Manual  TN  5-1300. 

Generally,  the  response  time  is  of  interest  to  the  user 
in  judging  the  adequacy  of  the  rigid  body  approximation  for  com- 
puting the  gross  response  of  the  structure.  A small  element  res- 
ponse time  in  comparison  to  the  gross  response  time  of  the 
structure  indicates  that  the  rigid  body  approach  is  adequate  to 
the  task.  In  cases  where  the  response  times  of  individual  ele- 
ments are  nearly  equal  to  that  of  the  gross  structural  motion, 
the  rigid  body  method  will  yield  conservative  results.  This  is 
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especially  true  when  the  overall  structure  responds  well  before 
the  back  wall  element  responds.  For  these  cases  which  are  not 
cormon,  a more  sophisticated  approach  would  be  required. 

The  computation  is  performed  for  elements  designed  to 
incipient  failure  or  to  the  post-failure  fragment  conditions. 

In  the  event  of  the  latter,  the  response  time  of  the  element  is 
used  as  an  upper  time  limit  on  the  rigid  body  response  as  no 
further  computations  are  required  after  the  wall  has  failed. 

5.4  Computer  Usage  and  Restrictions 

The  program  is  written  in  FORTRAN  IV  for  the  CDC  6600 
computer.  The  FORTRAN  coding  for  the  program  can  be  found  in 
Appendix  0. 

The  size  restrictions  imposed  by  the  dimensional  con- 
straints of  the  program  are  summarized  below; 


Item 

Maximum  Number 

Number  of  charges 

20 

Number  of  wall  elements 

4 

Number  of  soil  elen^nts 

15 

Number  of  force-time  stations 

1,000 

Number  of  output-time  stations 

1,000 

Number  of  loaded  surfaces 

5 

The  results  of  the  analysis  are 

stored  internally  in  the 

computer  until  the  numerical  integration  is  com|ileted  after  v^ihich 
they  are  retrieved  and  printed  ouL  There  is  no  limitation  on 
the  number  of  integration  time  steps  that  can  beutitizvxi  in  Ihe 
analysis.  However,  there  is  a limitation  on  the  amount  of  data 
that  can  be  stored  internally  in  the  computer.  This  limitation, 
established  by  the  dimensional  constraints  of  the  program,  deter- 
mines the  maximum  number  of  output  time  stations  allowed.  There- 
fore, caution  must  be  exercised  in  specifying  the  number  of 
integration  time  steps  to  be  skipped  between  output  time  stations. 
If  an  insufficient  number  of  integration  time  steps  av'e  skipped 
between  output  time  stations,  the  quantity  of  data  stored  in  the 
computer  will  exceed  the  dimensional  limitations  of  the  program. 
This  will  result  in  premature  termination  of  the  calculation. 
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5.5 


User's  Manual 


5.5.1  Introduction 

The  input  to  the  computer  program  is  presented  in  coded 
card  format  in  Sections  5.5.2  and  5.5.3. 

Section  5.5.2  contains  the  input  data  forms  required  for 
the  "Normal  Option"  of  the  proqram.  The  "Normal  Option"  mode  of 
the  proqram  is  applicable  only  to  the  structural  confiourations 
shown  in  Figure  12.  The  proqram,  in  this  mode,  considers  only 
the  triangular  (linear)  pressure-time  function  in  computing  the 
load  history. 

Section  5.5.3  contains  the  additional  data  forms  required 
for  the  "General  Structure"  and  "Special  Loading  Options".  Sec- 
tion 5.5.4  describes  the  input  terminator  cards.  The  composition 
of  the  data  decks  ^or  the  various  options  available  in  the  prooram 
is  presented  in  Section  5.5.5.  Following  this  is  Section  5.5.6 
which  contains  user  instructions  for  running  sevo'^ftl  computer 
analyses  consecutively  (Multiple  Job  Procession). 

5.5.2  Input  Data  Forms:  Normal  Option 

There  are  6 types  of  cards  used  to  specify  data  for  the 
“Normal  Option"  mode  of  the  program.  Each  type  of  card  is  des- 
cribee below  in  terms  of  data  fomiat,  defip.ition  and  fie-d  allo- 
cations. The  numbers  above  the  graphic  representation  of  e"'ch 
cord  identify  the  last  column  in  each  field  of  thi  ard,  I'n 
fields  designated  “I",  the  quantity  must  he  ri^. ht  adjostid  to 
the  last  column  in  the  field.  No  '®cimal  point  is  r_.,jired  for 
"I"  fonratted  input  In  fields  designated  “F"',  a deci'aal  point 
is  required;  however,  the  number  con  be  located  anywhere  vfithin 
the  field.  A plus  sign  for  a positive  quantity  is  not  required 
and  will  abort  the  execMtion  of  the  program.  Minus  signs  for 
negative  quantities  must  be  placed  ■>  Jie  first  blank  column  to 
the  left  of  the  nuniher.  The  cards  ? numbered  (Gat'd  Type  1,  2, 
etc.)  according  to  the  order  in  which  they  are  read  in  by  the 
program. 


Card  Type  1 - Structure  ilescription  Card  (Required) 


7f 

STRUCTURE  DESCRIPTION  CARO 

■'  y y / j 

This  card  may  contain  alphanumeric  inforwation  in  Columns 
2-71  that  will  be  printed  at  the  top  of  each  page  of  the 
output. 
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Card  Type  2 - Problem  Specification  "ard  (Required) 


5 

10 

15 

30 

35 

40 

45 

50 

55 

60 

65 

NP 

N 

NS 

m 

ICI 

ICAl 

s 

‘■2 

N 

H 

W 

E. 

L 

N 

F 

a 

("l”  FORMAT  - ALL  FIELDS) 


NP  = number  of  charges 
N = number  of  walls 
NS  = number  of  soil  elements 
NUM7M  = number  of  integration  time  steps 

ICl  = 0 output  limited  to  a tabulation  of  maximum 
displacements,  soil  pressure  and  foundation 
shear  and  moment 

= 1 program  prints  displacements,  soil  resistance 
forces  and  soil  bearing  pressures  at  every 
output  time  station 

ICAl  = 0 Normal  Option:  mass  of  the  structure,  loca- 
tion  of" the  center  of  gravity,  areas  of  loaded 
surfaces  and  the  location  of  their  centroids 
relative  to  the  center  of  gravity  are  calcu- 
lated f>’om  structure  geometry  entered  on 
Card  Type  4.  Pressui'e-time  histories  are 
calculated  using  structure  geometry  input  on 
Card  Type  4 and  charge  data  input  on  Card 
Type  5. 

~ 1 General  Structure  Option:  Internal  coinputa- 

tfon  of  aforementioned  data  is  bypassed  and 
the  inforretion  entered  using  C^’rd  Types  6 
and  8 through 


KDE12  = constant  used  for  changing  integration  time  step; 
normally,  the  tisae  step  used  in  the  coi^utation 
is  1/20  of  the  smallest  load  duration  on  any  one 
surface  of  the  structure.  If  this  option  is  ex- 
ercised, the  tine  step  will  be  multiplied  by  this 
factor.  The  integration  time  step  is  altered 
only  after  the  pressure  ©n  the  structure  has 
decayed  to  zero. 


NUMPT  = number  of  integration  time  steps  skipped  between 
output  time  stations.  NOTE:  NUMTM/NUMPT  £ 1000. 

NWALL  = 1 the  program  calculates  the  maximum  response 
time  of  back  wall  element  for  incipient 
failure  or  post-failure  fragments. 

NVEL  = 1 the  acceleration  and  velocity  of  the  struc- 
ture is  printed  at  every  output  time  station. 

NFDN  = 1 the  maximum  moment  on  the  foundation  at  the 
wall  face  and  the  maximum  shear  and  corres- 
ponding moment  at  a specified  distance  from 
the  wall  face  are  computed  by  the  program. 

The  specified  distance  used  in  the  computation 
is  15  percent  of  the  length  of  the  foundation 
overhang.  The  width  of  a wall  haunch  is  in- 
cluded in  the  computation,  if  one  is  present. 
This  option  applies  only  to  cantilever  wall 
type  barriers. 

NLOAD  = 0 Normal  Option:  loading  on  the  structure  will 
be  calculated  internally  using  the  structure 
geometry  entered  on  Card  Type  4 and  the 
charge  data  entered  on  Card  Type  5. 

NLOAD  > 0 Special  Loading  Option:  internal  computation 

of  the  loading  will  be  bypassed  and  the  load- 
ing data  will  be  entered  on  Card  Types  6 and 
12.  The  value  of  "NLOAD"  is  used  as  the 
number  of  loaded  surfaces  on  the  structure. 
When  this  option  is  not  used  in  conjunction 
with  the  "General  Structure  Option",  enter 
zero  for  the  parameter  "ICAl". 

Card  Type  3 - Soil  Properties  Card  (Required) 


10 

20 

30 

40 

50 

60 

70 

El 

E2 

SMTP 

/8x 

m 

(V'  FORMAT -ALL  FIELDS) 


El  = modulus  of  elasticity  (psi)  of  first  portion  of 
bilinear  stress-strain  curve 


Card  Type  3 (continued) 

E2  = modulus  of  elasticity  (psi)  of  second  portion  of 
bilinear  stress-strain  curve.  For  linear  stress- 
strain  curves,  enter  for  "E2"  the  value  of  "El". 

SMTP  = stress  (psi)  at  which  modulus  of  elasticity 
changes  (see  Figure  9,  page  22).  For  linear 
stress-strain  curves,  enter  a value  of 
10000.0  for  SMTP. 

f^  = coefficient  of  friction  between  soil  and  base 
of  structure  for  non-cohesive  soils  or  adhesion 
constant  (psf)  for  cohesive  soils 

u = Poisson's  ratio  for  soil 

6 = influence  coefficient  for  calculating  soil 
spring  constants: 

(3^  is  the  coefficient  for  the  horizontal 
spring  constant 

62  is  the  coefficient  for  the  vertical  spring 
constant 

Card  Type  ^ - Structure  Geometry  Card  (Required)  (See 
Figure  12) 


10 

20 

30 

40 

50 

60 

70 

r ■' 

j TW 

C9 

L 

HW 

SSB 

HB 

HAUNH 

("f"  format- all  FIELDS) 


TW  = back  wall  thickness  (in) 

C9  = ratio  of  foundation  tliickness  to  back  wall 
thickness  (TS/TW) 

L - length  of  hack  wall  (ft) 

HW  = height  of  back  wall  (ft) 

SSB  = ratio  of  length  of  loaded  area  of  foundation 
to  total  length  of  foundation  (Ljr/B) 
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HB  = ratio  of  height  of  back  wall  to  length  of 
foundation  (HW/B) 

HAIINH  = width  of  back  wall  haunch  (in) 

Card  Type  5 - Charge  Data  Cards  (Required  for  Normal 
Option  only)  (See  h’gure  13) 

Normal  Option:  computation  of  pressure- 
time  histories  is  performed  internally 
by  the  program  utilizing  the  data  on  these 
cards  and  the  TNT  explosive  data  contained 
in  the  program 


^ 

20 

30 

40 

50 

60 

70 

< 

(T 

W 

ifad) 

i^(2) 

i„(3) 

1 

h 

("f"  format- all  FIELDS) 


Ra  = distance  from  center  of  charge  to  rear  face  of 
back  wall  (ft) 

W = charge  weight  (lbs) 

i^(l)  = unit  blast  impulse  on  back  wall  (psi-ms) 

ib(2)  = unit  blast  impulse  on  foundation  (psi-ms) 

ib(3)  = unit  blast  impulse  on  roof  (psi-ms) 

I ~ minimum  distance  fr/m  charge  to  an  adjacent 
wall  (ft) 

h « height  of  charge  above  foundation  (ft) 

NOTE:  One  "Charge  Data  Card"  is  required  for  each  ex- 
plosive charge. 
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Card  Type  7 - Back  Wall  Structural  Details  Card 

(Optional;  may  be  used  with  Normal  and 
Special  Loading  Options) 

This  option  is  used  if  the  maximum  response  time  of  the 
back  wall  element  is  required.  The  peak  response  time 
is  calculated  for  the  incipient  failure  or  the  post- 
failure fragment  conditions.  In  the  event  of  the  post- 
failure fragment  condition,  this  quantity  is  used  as  an 
upper  time  limit  on  the  response  of  the  structure. 

The  dimensions  and  parameters  entered  on  this  card  are 
defined  and  illustrated  in  Chapter  5 of  Reference  1. 
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20 
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40 
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Py 

Ph 

X 

y 

^ds 

(KLM)u 

("f"  FORMAT- all  FIELDS) 


lifj  = summation  of  unit  blast  impulses  on  back  wall 
element  (psi-ms) 

d^  = distance  between  the  centroids  of  the  tension 
and  compression  reinforcement  (in) 

Py  = reinforcement  ratio  vertical 

p^  = reinforcement  ratio  horizontal 
X or  y = yield  line  location  (in) 

fds  = dynamic  steel  design  stress  (psi) 

(KLM)jj  = plastic  load  mass  factor 

Vf  = velocity  of  post-failure  fragments  (in/ms). 

This  quantity  is  required  for  elements  designed 
to  the  post-failure  fragment  condition  only. 

NOTE: 

1.  If  the  structure  being  analyzed  is  a cantilever 
wall,  enter  only  Ei[j,  d^.,  Py,  f^s* 

2.  To  use  this  option,  a value  of  1 must  be  entered 
for  the  parameter  "NWALL"  on  Card  Type  2. 
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5.5.3  Input  Data  Forms:  General  Structure  and  Special 
Loading  Options 

There  are  6 additional  types  of  cards  (Card  Types  6 and 
8 through  12)  that  are  required  for  the  "General  Structure  Option" 
and  two  additional  types  of  cards  (Card  Types  6 and  12)  that  are 
required  for  the  "Special  Loading  Option".  These  cards  are  re- 
quired in  addition  to  Card  Types  1 through  4.  Each  type  of  card 
is  described  on  the  following  pages  in  the  manner  specified  in 
Section  5.5.2.  Card  Types  6 and  12  are  utilized  in  lieu  of  Card 
Type  5 ("Charge  Data  Card")  to  enter  the  loading  data.  Card  Types 
8 through  11  are  used  to  enter  the  inertial  properties  and  the 
geometry  of  the  structure.  The  data  on  Card  Types  6 and  9 through 
12  is  used  by  the  program  to  compute  the  load  history  on  the 
structure. 

Card  Type  6 is  used  to  enter  the  arrival  time  of  the 
blast  wave  on  the  structure  and  the  smallest  duration  of  the 
loading  produced  by  any  one  explosive  charge  on  any  surface. 

Card  Type  8 is  used  to  enter  the  inertial  properties 
of  the  structure  and  the  location  of  the  center  of  gravity  of  the 
structure.  The  dimensions  required  to  locate  the  center  of  grav- 
ity are  shown  in  Figure  14. 

Card  Types  9 through  11  contain  geometry  data  pertaining 
to  the  loaded  surfaces  of  the  structure.  As  discussed  in  Sec- 
tion 2,  the  analysis  treats  the  structure  as  a rigid  body,  con- 
strained to  move  in  one  plane;  therefore,  only  loading  and  ge- 
ometry data  pertaining  to  surfaces  normal  to  the  plane  of  motion 
of  the  structure  have  to  be  included  in  the  input  data  deck. 

Card  Type  9 is  used  to  enter  the  areas  of  those  surfaces 
directly  exposed  to  the  blast  pressures. 

Card  Type  10  is  used  to  enter  the  location  of  the  centroid 
of  a loaded  surface  relative  to  the  center  of  gravity  of  the 
structure. 

Card  Type  11  is  used  to  enter  the  horizontal  and  vertical 
components  of  a unit  vector  normal  to  a loaded  surface.  The  vec- 
tor defined  on  this  card  must  be  directed  towards  tfie  surface. 

Areas  for  all  loaded  surfaces  can  be  entered  on  Card 
Type  9 whereas  one  each  of  Card  Types  10  and  11  is  required  for 
each  loaded  surface.  Card  Types  10  and  11  are  entered  in  two 
separate  groups,  that  is,  one  group  contains  all  Card  Types  10 
and  the  other  contains  all  Card  Types  11. 
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Card  Type  1J?  contains  loading  data  that  includes  the  aver- 
age impulse  on  the  surface,  the  arrival  time  of  the  blast  wave  on 
the  surface  and  the  sum  o^  the  arrival  time  of  the  blast  wave  plus 
the  load  duration  on  the  surface.  The  surface  loading  data  is 
entered  in  groups.  Each  group  contains  the  surface  loadings  pro- 
duced by  one  explosive  charge.  In  each  group,  one  "Surface  Load- 
ing Data  Card"  (Card  Type  12)  is  entered  for  loaded  surface. 

Figure  16  on  page  51  illustrates  the  input  data  deck  for 
the  "Special  Loading  Option".  The  figure  shows  several  groups  of 
data  cards,  each  labeled  with  a charge  number  and  the  encircled 
number  12.  Each  group  contains  several  cards,  with  each  card  con- 
taining the  loading  data  (average  impulse,  arrival  time  and  load 
duration)  for  one  loaded  surface.  For  the  "Special  Loading  Op- 
tion", the  surface  loading  data  cards  must  be  ordered  (within  each 
group)  as  follows: 

Order  of  Data  Card  Within  Group  Surface 

First  Card  Back  Wall 

Second  Card  Floor  or  Foundation 

Third  Card  Roof  (if  present) 

If  the  data  cards  are  not  in  the  correct  sequence,  the  results  of 
the  analysis  will  be  erroneous. 

Figure  17  on  page  52  illustrates  the  input  data  deck  for 
the  'General  Structure  Option".  The  figure  shows  several  groups 
of  data  cards  (one  group  each  for  Card  Types  9,  10  and  11  with 
the  card  numbers  encircled)  in  addition  to  the  groups  containing 
the  loading  data  cards  (Card  Type  12).  Although  the  groups  must 
be  entered  in  the  same  order  as  shown  in  the  figure,  within  each 
group  a fixed  order  for  entering  tlie  data  for  the  loaded  surfaces 
is  not  required.  However,  the  order  in  v«ihich  the  data  is  entered 
must  be  the  same  for  each  group;  that  is,  if  the  areas  entered  on 
Card  Type  9 are  for  the  roof,  floor  and  back  wall  in  that  order, 
then  the  surface  data  entered  in  the  remaining  gj'oups  (Card  Types 
10,  11  and  12)  must  also  be  in  the  same  order.  Changing  the 
order  of  the  input  in  each  group  will  cause  erroneous  analysis 
resul  ts . 
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Card  Type  6 - Time  Step  Data  Card  (Required  for  both 
options) 


10 

20 

TASM 

TOLG 

(“f"  format -all  FIELDS)' 

TASM  = time  of  arrival  of  the  blast  wave  on  the  struc- 
ture. This  quantity  is  determined  by  computing, 
as  described  in  Section  3.3,  the  time  for  the 
blast  wave  from  any  charge  to  traverse  the  least 
distance  from  that  charge  to  a point  on  the 
structure.  The  smallest  quantity  computed  for 
any  explosive  charge  is  the  time  of  arrival  of 
the  blast  wave  on  the  structure.  This  value  is 
used  as  the  first  integration  time  station  in 
the  analysis. 

TOLG  = smallest  duration  of  loading  produced  by  any  one 
explosive  charge  on  any  surface  of  the  structure. 


Card  Type  8 - Mass  Data  Card  (Required  for  General 

Structure  Option  only)  (See  Figure  14). 
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("f“  format  - ALL  FIELDS) 


WT  = weight  of  structure  (lbs) 

I = mass  moment  of  inertia  (Ibs-sec^/in) 

XB  = horizontal  distance  in  inches  from  center  of 
gravity  to  rear  face  of  back  wall  element 

YB  “ vertical  distance  in  inches  from  center  of 
gravity  to  top  of  the  foundation  slab 

XR  = horizontal  distance  in  inches  from  center  of 
gravity  to  the  left  end  of  the  foundation 

NOTE:  The  signs  of  the  dimensions  entered  on  this  card 
must  be  consistent  with  the  sign  convention  of 
the  x-y  coordinate  system  shown  in  Figure  14. 
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Card  Type  9 - Loaded  Surface  Area  Card  (Required  for 
General  Structure  Option  only) 


10 


20 


AFACE(i)UfACE(2) 


30 


40 


AFACE(3)|AFACE(4) 


AFACECS) 


( F FORMAT -ALL  FIELDS) 


AFACE(i)  = area  of  surface  in  square  inches.  The  area 
is  to  be  entered  for  each  loaded  surface 
on  the  structure. 

Card  Type  10  - Loaded  Surface  Centroid  Card  (Required 
for  General  Structure  Option  only) 


10 

20 

XCOROd) 

YCORD(I) 

("f"  format- all  FIELDS)’ 

XCORO(i)  = horizontal  distance  in  inches  from  the 
center  of  gravity  to  the  centroid  of  the 
surface 

YCORD(i)  = vertical  distance  in  inches  from  the  center 
of  gravity  to  the  centroid  of  the  surface 

NOTE: 


1.  One  "Loaded  Surface  Centroid  Card"  is  required 
for  each  loaded  surface  of  the  structure. 

2.  The  signs  of  the  dimensions  must  be  consistent 
with  the  sign  convention  of  Wie  x-y  coordinate 
system  shown  in  Figure  14. 

Card  Type  11  - Loaded  Surface  iiormal  Vector  Cav'd 

(Required  for  General  Structure  Option 
only) 


10 
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UVECTCI) 

UVECT(21 

i-y. . 

("F**  FORMAT  - ALL  FIELDS) 
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YCORD(W) 


Surface  oraos  referred  to 
ore  perpendiculor  to  x*y  pione. 


Figure  14,  Center  of  gravity  location:  Card  Type  8. 


UVECT(l)  = horizontal  component  of  a unit  vector  normal 
to  the  surface 

UVECT(2)  = vertical  component  for  a unit  vector  normal 
to  the  surface 

NOTE: 


1.  One  “Loaded  Surface  Normal  Vector  Card"  is  re- 
quired for  each  loaded  surface  of  the  structure. 

2.  The  signs  of  the  vector  components  must  be  con- 
sistent with  the  sign  convention  of  the  x-y 
coordinate  system  shown  in  Figure  14. 

Card  Type  12  - Surface  Loading  Data  Card  (Required 
for  both  options) 
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(“f“  format  - all  FIELDS) 


i{j  = unit  blast  impulse  on  surface  (psi-ms) 
t;^  = arrival  time  of  blast  wave  on  surface  (sec) 

t^O  “ arrival  time  of  blast  wave  plus  duratimi  of 
loading  on  surface  (sec) 

P2/P1  = ratio  of  final  pressure  to  initial  pressure 
on  surface 

5.5>4  Input  Teminator 

Following  the  input  data  cards  are  the  input  teminator 
cards  which  signal  the  program  that  no  further  data  is  to  be  en- 
tered. Tne  program  will  then  execute  a noimvsl  exit.  The  input 
terminator  cards  are: 

1.  One  (1)  blank  card 

2.  One  (1)  card  with  a negative  integer  in 
Colisans  1 - 5. 
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(^KWALL  STRUCT®  i 
DETAILS  CARD 


CHARGE  DATA 
CARDS  ^ 


STRUCTURE 
GEOMETRY  CANO 


INPUT  TERMINATOR 
CARDS 


j^lL^PROPERTIES  @ 

fwOBLEM  SPEC.  (2 

[card  ^ 

) 

title  CARO  0 

note: 

ENCIRCLED  NUMBERS  REFER 
TO  CARD  IDENTIFICATION 
NUMBERS  IN  USER’S  MANUAL 


Figure  15  Input  data  deck:  Konnal  Option. 


. m 


The  “Special  Loading  Option"  is  used  to  acconsno- 
date  rectanqular  or  trapezoidal  load  histories 
in  the  analysis.  It  is  also  used  in  conjunction 
with  the  "General  Structure  Option"  to  accoitmo- 
date  structures  of  any  arbitrary  configuration. 
This  option  requires  considerably  snore  input 
data  than  the  "No»*Bial  Option". 


The  data  deck  for  the  "Special  Loadino  Option" 
is  sh»v*n  in  Figut’e  16.  The  back  wall  element 
response  con^Jutation  can  also  be  included  in 
this  Option,  if  so  desired. 


3 . Ge*"  ^ra  1 S true ture  Option 


The  “General  Structure  Option"  is  utilized  for 
structures  of  any  arbitrary  configuration.  It 
requires  the  greatest  asaunt  of  input  of  all  the 
options  available  in  the  program.  In  a>:  iition 
to  the  load  history  d-5ta  necessary  for  *ae 
"Special  loading  Option",  all  of  thr  data  pertain- 
ing to  the  inertial  properties  and  the  geometry 
of  the  structure  are  also  required.  The  data 
deck  for  the  "General  Structure  Option"  is  showr? 
in  Fi§u»'0  17. 


$.5.6  ffeiltipTe  dob  Processing 


Several  pi-oMers  can  be  processed  in  one  ce?sfHiter  run  by 
Simply  stacking  the  input  data  decks  for  each  proble?^  sne  after 
the  other.  The  input  terrainator  is  tlien  pUee4  after  the  last 
data  deck. 


Figure  16.  Input  data  deck:  Special  Loading  Option. 
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SECTION  6 


COMPUTER  PROGRAM  OUTPUT 


6.1  Introduction 

This  section  presents  a description  of  the  computer  pro- 
gram printed  output.  Included  also  is  a discussion  on  the  utili- 
zation of  the  output  with  emphasis  on  that  portion  related  to  the 
design  of  foundation  slabs. 

Samples  of  the  printed  output  are  presented  in  Appendix  D. 

6.2  Description  of  Output 

6.2.1  Summary  of  Output 

The  computer  program  output  consists  of  three  parts: 

(1)  a surmary  of  the  input  parameters  used  in  the  analysis;  (2) 
the  loading  on  the  structure;  and  (3)  the  response  of  the  struc- 
ture to  the  applied  loads. 

The  printing  sequence  has  been  designed  to  provide  the 
input  summary  as  a single  block  followed  by  a series  of  data 
blocks  containing  the  blast  loads  on  the  structure.  These  data 
are  arranged  to  present  the  loads  associated  with  each  explosive 
charge  separately.  Descriptive  data  about  the  charge  are  printed 
first.  This  is  followed  by  the  pressures  and  forces  acting  on 
each  loaded  face  of  the  structure  in  sequence.  The  second  charge 
is  then  described  and  the  loads  acting  on  the  loaded  faces  due  to 
this  charge  are  printed.  Finally,  the  structure's  response  to 
the  load  history  is  given. 

6.2.2  Summary  of  Input 

The  input  summary  includes  the  program  control  parameters, 
the  basic  dimensions  of  the  structure  and  the  elastic  properties 
of  the  soil.  Both  the  horizontal  and  vertical  spring  constants 
for  the  soil  are  presented.  Where  the  soil  behavior  is  bilinear, 
both  moduli  of  elasticity  are  given.  In  addition,  the  weight, 
mass,  mass  moment  of  inertia,  the  significant  dimensions  (illus- 
trated in  Figure  14)  locating  the  structure  center  of  gravity  and 
the  integration  time  step  used  in  the  analysis  are  tabulated. 
Presentation  of  the  input  data  in  this  form  facilitates  checking 
and  tabulation  of  the  problems. 
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6.2.3  Blast  Loads  on  Structure 

The  loads  produced  by  each  explosive  charge  acting  on  each 
wall  are  presented  in  sequence.  Where  a multi-charge  arrangement 
occurs,  a separate  set  of  load  histories  is  presented  for  each 
charge  after  a description  of  that  charge.  The  charge  descriptive 
output  includes  the  weight  of  explosive  and  its  location  relative 
to  the  back  wall,  floor  and  side  walls  of  the  structure.  This 
data  along  with  the  average  impulse  acting  on  each  wall  is  omitted 
when  the  "Special  Loading  Option"  is  exercised. 

The  data  set  for  each  loaded  face  includes  the  average 
impulse,  the  surface  area,  the  location  of  the  centroid  of  the 
area  with  respect  to  the  structure  center  of  gravity,  the  initial 
and  final  pressures  on  the  surface  and  the  times  at  which  they 
occur.  Time  histories  of  the  load  and  the  moment  produced  by  the 
load  relative  to  a horizontal  axis  through  the  center  of  gravity 
are  also  presented. 

Following  the  above  data,  a summation  of  all  the  horizontal 
and  vertical  loads  acting  on  the  structure  together  with  the  re- 
sultant moment  with  respect  to  the  center  of  gravity  is  presented. 

6.2.4  Response  of  the  Structure 

When  the  response  time  option  is  exercised,  the  input  data 
associated  with  the  design  of  the  back  wall  element  is  printed 
followed  by  the  peak  response  time  of  this  wall. 

In  general,  the  output  contains  displacement,  velocity  and 
acceleration  time  histories  of  the  structure  resolved  into  hori- 
zontal and  vertical  components  at  the  center  of  gravity  and  an 
angular  component  about  an  axis  which  passes  through  the  center  of 
gravity  and  is  perpendicular  to  the  plane  of  motion.  Similar  data 
are  presented  for  the  foundation.  The  foundation  quantities  are 
found  in  a column  headed  "XFDN". 

The  format  of  the  printout  presents  several  groups  of 
tabulations.  The  first  tabulation  includes  the  displacements  of 
the  center  of  gravity,  horizontal  displacement  of  the  foundation 
and  the  soil  "Resisting  Forces".  These  latter  are  summations  of 
the  forces  in  the  horizontal  and  vertical  soil  elements  and  the 
moments  of  these  forces  about  the  center  of  gravity  of  the  struc- 
ture. The  second  tabulation  contains  the  acceleration-  and 
velocity-time  histories  of  the  structure.  The  horizontal  velocity- 
time history  of  the  foundation  is  also  included  under  the  heading 
"XFDN". 
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The  third  data  group  consists  of  the  variation,  with  time, 
of  the  bearing  pressures  in  the  soil.  As  discussed  previously, 
the  soil  is  represented  in  the  analysis  as  a series  of  discrete 
elements  attached  to  the  foundation  at  equally  spaced  intervals. 
The  quantities  printed  are  the  bearing  pressures  at  the  soil  ele- 
ment attachment  points.  A value  of  zero  for  the  bearing  pressure 
indicates  that  the  foundation  at  that  point  is  moving  upward  away 
from  the  soil. 

The  response-time  histories  are  followed  by  a tabulation 
of  significant  response  parameters  which  include  the  maximum  up- 
ward deflection  of  the  toe  of  the  foundation.  The  toe  is  that 
part  of  the  foundation  resisting  the  overturning  motions  of  the 
structure.  If  the  toe  leaves  the  ground,  this  deflection  will  be 
the  maximum  upward  deflection.  If  the  toe  remains  below  the 
ground,  it  will  actually  be  the  minimum  downward  deflection. 

If  the  structure  reaches  its  peak  rotation,  the  tabulation 
will  include  the  maximum  horizontal  and  vertical  displacements  of 
the  center  of  gravity,  the  overturning  angle  (rotation  of  struc- 
ture at  which  overturning  occurs),  the  maximum  rotation  of  the 
structure,  and  the  ratio  of  the  maximum  rotation  to  the  overturn- 
ing angle.  The  output  includes  the  maximum  upward  displacement  of 
the  center  of  gravity,  the  maximum  bearing  pressure  in  the  soil 
and  the  foundation  design  loads  for  a cantilever  wall  type  barrier. 
If  the  peak  response  has  not  been  attained,  this  output  will  be 
omittedand  a message  printed,  indicating  that  more  time  is  re- 
quired in  the  analysis.  If  this  occurs,  the  number  of  integration 
time  steps  must  be  increased  and  the  analysis  performed  again  un- 
til the  maximum  response  is  reached. 

The  volume  of  output  can  be  limited  by  either  deleting 
whole  sections  or  increasing  the  number  of  integration  time  steps 
skipped  between  printouts. 

6.3  Use  of  Output 

6.3.1  General 


The  output  of  this  computer  program  may  be  used  for  three 
primary  purposes:  (1)  to  evaluate  the  response  of  the  structure; 
(2)  to  detect  errors  in  either  the  input  data  or  the  analysis; 
and  (3)  to  design  foundation  slabs. 


The  first  two  require  no  explanation.  The  third,  however, 
is  discussed  in  the  section  that  follows  and  treated  in  detail  in 
Appendix  C. 
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6.3.2  Foundation  Design 

Design  procedures  for  foundations  of  protective  structures 
vary  according  to  the  configuration  of  the  structure.  Basically, 
all  protective  structures  depend,  to  some  extent,  on  the  soil  be- 
neath the  foundation,  to  resist  the  structure  overturning  and/or 
sliding  motions.  The  extent  of  this  dependence  on  the  soil  is  a 
function  of  the  configuration  and  size  of  the  structure,  and  where 
applicable,  the  configuration  and  size  of  larger  structures  to 
which  the  protective  structure  is  tied. 

Large  multi -cell  barricades,  such  as  the  one  shown  in  Fig- 
ure 5,  will  experience  gross  rotational  and  sliding  motions  under 
the  action  of  the  blast.  However,  the  massive  blast  wall  elements 
and  the  equally  massive  foundation,  combined  with  the  substantial 
width  of  the  foundation,  will  severely  limit  the  gross  motions  of 
the  structure  to  the  extent  that  they  are  no  longer  a factor  to  be 
considered  in  the  design.  This  is  also  true  for  barricades  with  a 
single  row  of  cells.  Foundations  for  these  structures  are  designed 
primarily  to  develop  the  blast  wall  elements.  Long  foundation  ex- 
tensions, projecting  well  beyond  the  extent  of  the  cell,  are  not 
required.  The  length  of  the  foundation  extension,  if  one  is  re- 
quired, is  usually  established  by  the  anchorage  required  for  the 
reinforcing  steel  in  concrete. 

Barrier  walls  and  their  foundations  that  are  integral  with 
larger  structures  (such  as  explosive  manufacturing  or  storage  fa- 
cilities) experience  large  rotational  motions  but  are  generally 
not  susceptible  to  overturning.  The  resistance  to  overturning  is 
generated  by  mobilizing  the  mass  of  the  larger  structure  through 
the  foundation  and  other  components  of  the  larger  structure  (such 
as  the  foundation  walls  shown  in  Figure  4).  In  these  situations, 
the  foundation  design  is  based  on  developing  the  blast  wall  ele- 
ments. In  some  cases,  depending  upon  the  configuration  of  the 
larger  structure  and  the  location  of  the  barrier  wall  within  it, 
some  portion  of  the  resistance  to  the  overturning  motions  may  have 
to  be  generated  in  the  soil  beneath  the  foundation.  Two  examples 
of  this  condition  are:  an  excessively  tall  interior  barrier  wall 
restrained  only  by  the  foundation  of  the  structure  and  a barrier 
wall  located  at  the  end  of  a structural  floor  slab. 

In  these  cases,  the  foundation  is  initially  designed  to 
develop  the  blast  wall  element.  Then,  in  order  to  determine  the 
peak  soil  pressures  acting  on  the  foundation,  an  overturning  anal- 
ysis is  performed.  The  analysis  considers  the  barrier  wall  acting 
together  with  some  effective  segment  of  the  structural  foundation. 
The  foundation  is  then  checked  to  insure  that  it  can  resist  the 
moments  and  shears  produced  by  the  peak  soil  pressures. 
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Simple  cantilever  wall  barriers  and  single  cell  cubicles 
generally  experience  relatively  large  gross  motions,  under  the 
action  of  blast  loads,  which  are  resisted  entirely  by  the  com- 
pression forces  developed  in  the  soil  beneath  their  foundation 
extensions  (see  Figures  1 and  2).  They  are  therefore  highly  sus- 
ceptible to  overturning.  The  compression  forces  in  the  soil  im- 
part shear  and  bending  stresses  in  the  foundation  extension  for 
which  the  foundation  slab  must  be  designed.  For  single  cell  cu- 
bicles, the  foundations  are  designed  both  to  develop  the  blast 
wall  elements  and  to  resist  the  compression  forces  in  the  soil. 

On  the  other  hand,  the  foundation  extension  for  a cantilever  wall 
barrier  is  designed  primarily  to  resist  the  soil  bearing  pressures. 
This  comes  about  because  the  response  of  the  cantilever  wall  bar- 
rier on  the  soil  closely  approximates  a true  rigid  body  response. 
Consequently,  the  only  significant  loadings  on  the  foundation  ex- 
tension are  the  soil  pressures  generated  as  the  structure  responds 
to  the  blast  loads. 

Designing  the  foundation  of, a protective  structure  to 
develop  the  strength  of  the  blast  wall  elements  is  a relatively 
simple  task.  Basically,  the  procedure  entails  designing  the 
foundation  slab  for  the  ultimate  bending  moment  capacity  of  the 
wall . 


Designing  the  foundation  for  the  overturning  motions  is 
accomplished  utilizing  the  computer  program  presented  in  Section  5. 
As  noted  in  Section  6.2.4,  the  bearing  pressure-time  history  is 
included  in  the  printed  output  of  the  program. 

The  foundation  is  designed  to  resist  the  peak  bearing 
pressures  generated  as  the  structure  responds  to  the  blast  loads. 

The  foundation  is  usually  a very  stiff  element  and,  therefore, 
the  time  required  for  the  bearing  pressures  to  build  to  a peak  is 
much  greater  than  the  period  of  the  element.  Consequently,  the 
foundation  must  be  capable  of  resisting  the  full  load  developed 
in  the  soil.  If  the  foundation  were  less  rigid,  it  would  not 
provide  sufficient  strength  to  resist  the  overturning.  Sometimes 
a trade-off  between  the  length  of  the  extension  and  the  mannitude 
of  the  gross  rotations  can  be  made  which  will  result  in  larger 
rotations  but  a thinner  foundation  slab. 

The  procedure  and  criteria  for  designing  the  foundation  to 
resist  the  bearing  pressures  «n  the  soil  are  presented  in  Appendix  C. 
Included  also  are  design  examples  which  illustrate  the  procedure. 
The  procedure  Is  applicable  primcrily  to  the  design  of  protective 
structures  (such  as  cantilever  wall  barriers  and  single  cell  cu- 
bicles) where  the  stability  of  the  structure  is  a critical  factor. 
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Briefly,  the  procedure  calls  for  designing  the  foundation  for  the 
range  of  properties  (specified  in  Tables  1 and  2 of  Section  4)  for 
the  particular  soil  at  the  construction  site.  This  is  done  because 
the  soils  data  normally  available  to  the  designer  is  at  best  a 
gross  estimate  of  the  actual  properties.  Therefore,  the  soil  con- 
ditions producing  the  largest  overturning  motions  and  foundation 
loadings  have  to  be  considered  in  the  design. 

In  some  situations  where  soil  stresses  are  quite  large, 
backup  structural  elements  may  be  required  to  support  the  founda- 
tion. This  condition  usually  occurs  in  the  design  of  very  tall 
structures  where  the  simple  type  foundation  slab  would  be  exces- 
sively long  and  thick.  In  these  cases,  buttress  walls  and  foun- 
dation beams  can  be  provided  to  reduce  the  slab  thickness.  An 
example  of  this  is  shown  in  Figures  C.4  and  C.5  of  the  Appendix. 
These  elements  must  be  designed  for  the  reactions  of  the  foundation 
slab  when  the  peak  bearing  pressures  are  developed  beneath  the 
structure. 

The  design  of  foundations  for  protective  structures  must 
also  satisfy  conventional  design  criteria.  The  conventional  cri- 
teria, applicable  to  the  normal  working  load  condition  (dead  and, 
if  applicable,  live  load),  limits  the  total  settlement  of  the 
structure  and  attempts  to  eliminate  differential  settlements  in 
various  parts  of  the  structure.  To  satisfy  this  criteria,  it  is 
necessary  that  the  structure  foundation:  (1)  transmit  the  load 

of  the  structure  to  a soil  stratum  of  sufficient  strength;  and 
(2)  spread  the  load  over  a sufficiently  large  area  of  that  stratum 
to  minimize  the  bearing  pressure.  If  an  adequate  soil  is  not 
found  immediately  below  the  structure,  deep  foundations  such  as 
piles,  are  used  to  transmit  the  loid  to  deeper  firmer  layers. 

flenerally,  the  foundation  of  a protective  structure  de- 
signed for  the  effects  of  the  blast  on  the  structure,  will  meet 
the  requirements  of  the  conventional  criteria.  However,  situations 
may  arise  in  which  the  conventional  criteria  will  control  the  de- 
sign of  the  foundation.  An  example  of  this  is  the  structure  of 
Example  C.2  of  Appendix  C.  In  this  problem,  the  plan  size  of  the 
foundation  is  dictated  by  the  requirement  for  naintaining  the 
bearing  pressure  below  the  specified  allowable  bearing  pressure 
under  the  weioht  of  the  structure.  As  a result,  the  foundation 
extension  required  is  larger  than  that  which  is  needed  to  prevent 
the  structure  from  overturning. 
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APPENDIX  A 


PROCEDURE  FOR  CALCULATING  THE  AVERAGE  IMPULSE  LOAD 

ON 

THE  FOUNDATION  SLAB  OF  A CANTILEVER  WALL  BARRIER 
A.1  General 


This  appendix  presents  the  data  and  the  procedure  for  com- 
puting the  average  impulse  loads  on  the  foundation  slab  of  a can- 
tilever wall  barrier. 

The  method  of  calculating  the  average  blast  impulses  was 
developed  by  using  a theoretical  procedure  based  on  semi -empirical 
blast  data  and  on  the  results  of  response  tests  of  reinforced  con- 
crete slabs. 

The  parameters  which  are  necessary  to  determine  the  average 
impulse  loads  are  the  structure  size,  charge  weight  and  charge  lo- 
cation. Figure  A.l  shows  the  parameters  necessary  to  determine  the 
average  impulse  loads. 
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Figure  A.l  Cantilever  wall  barrier  configuration  and 
parameters . 
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Because  of  the  wide  range  of  required  parameters,  the  pro- 
cedure for  the  determination  of  the  impulse  loads  was  programmed 
for  solution  on  a digital  computer.  The  results  of  these  calcula- 
tions, presented  in  Figures_A.2  through  A. 13,  give  the  scaled 
average  unit  blast  impulse  ik  as  a function  of  the  parameters  de- 
fined .in  Figure  A.l.  Each  illustration  is  for  a particular  com- 
bination of  values  of  Z/L  and  L/Lp. 

Because  of  the  limitation  in  the  range  of  the  test  data  and 
the  limited  number  of  values  of  the  parameters  given  in  the  impulse 
charts,  extrapolation  of  the  data  given  in  Figures  A. 2 through  A. 13 
may  be  required  for  some  of  the  parameters  involved.  On  the  other 
hand,  the  limiting  values  as  given  in  the  charts  for  other  parameters 
will  not  require  extrapolation.  The  following  are  recoimiended  pro- 
cedures which  will  be  applicable  in  most  cases  for  either  extra- 
polation, or  establishing  the  limits  of  impulse  loads  corresponding 
to  values  of  the  various  parameters  which  exceed  the  limits  of  the 
charts : 

a.  Plot  curve  of  values  of  k versus  Zp  for  constant 
values  of  R/\/Lp.  Lp/h,  L/Lp  and  ?/L.  Extrapolate 
the  curve  to  include  the  value  of  corresponding 
to  the  value  of  Zp  required. 

b.  Extrapolate  given  curve  for  constant  values  of  Zp, 

Lp/h,  L/Lp  and  ?/L,  to  include  value  of  i|j  corres- 
ponding to  the  value  R/^/Lp  required. 

c.  Plot  curve  of  values  of  T.  versus  Lp/h  for  constant 
values  of  Rfl/Lp,  Zp,  L/Lp\nd  :/L.  Extrapolate 
curve  to  include  value  of  ij^  corresponding  to  the 
value  of  Lp/h  required. 

d.  Plot  curves  of  values  of  T.  versus  L/Lp  for  constant 
values  of  R/\/Lp,  2p.  Lp/h  and  ;/L.  Extrapolate 
curve  to  include  value^of  i^,  corresponding  to  the 
value  of  L/Lp  required. 

e.  Values  of  corresponding  to  values  of  l/L  less 
than  0.10  snail  be  taken  as  equal  to  those  corres- 
ponding to  l/l  - 0.10. 
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A. 2 Computation  Procedure  and  Sample  Problem 

The  computation  of  the  average  impulse  load  proceeds  as 

follows: 

Problem;  Determine  the  average  impulse  load  on  the  foundation 
slab  of  a cantilever  wall  barrier. 

Procedure: 


I ^ 


Step  1 . Determine  the  following: 

a.  Charge  weight 

b.  Structure  dimensions*  L,  Lp 

c.  Charge  location  parameters  h,  I 

Step  2.  Apply  a 20%  safety  factor  to  the  charge  weight. 

Step  3.  Calculate  the  chart  parameters,  7/L,  L/Lp, 

RA/Lp,  Lp/h  and  the  scaled  distance  Zp. 

Note:  Use  of  the  average  impulse  charts  may 
require  interpolation  in  many  cases. 
Interpolation  may  be  achieved  by  inspec- 
tion for  the  scaled  distance,  Zp,  and  by 
a graphical  procedure  for  the  chart 
parameters  Lp/h,  L/Lp  and  l/l  using 
2 cycle  by  2 cycle  logarithmic  graph 
paper.  The  following  proc  edure  will 
illustrate  the  interpolation  of  all 
three  chart  parameters. 

Step  4.  Determine  and  ta^bulate  the  values  of  the  average 
scaled  impulse  ik  from  the  impulse  charts  for 
the  required  Ra/^c  and  Zp  and  the  followino 
variables; 

Lp/h  = 1.25,  1.75,  2.5  and  4.0 
L/Lp  = 1.00,  3.00,  6.50  and  10.0 
l/l  = 0.10,  0.25  and  0.50 
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FOR  REOUIRED  R^^/Lp  AND  Z 


1.00  3.00  6.50  10.0  1.00  3.00  6.50  10.0  1.00  3.00  6.50  10.0 


Figure  A.  2 A.3 
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Step  5 a.  Prepare  three  2-cycle  loq-lon  charts  with 
Lp/h  as  the  lower  abscissa,  L/Lp  as  the 
upper  abscissa,  and  i^  as  the  ordinate  (one 
chart  for  each  of  the  VI  ratios).  On  each 
chart  for  constant  Ifl  and  Zp,  plot  iK  versus 
Lp/h  for  all  L/Lp  values. 

b.  Usinq  chart  for  ?/L  = 0.10,  read  values  of 
versus  L/Lp  for  required  Lp/h.  Tabulate 
results. 


FOR  r/L  = n.lO  AND  REOUIRED  Lr/h 


4 


c.  Repeat  Step  5b  for  charts  Z/L  = 0.25  and  0.50. 
Tabulate  results. 

d.  On  each  IIL  chart,  plot  Tj^  versus  L/Lp  from 
Steps  5b  and  5c. 

e.  On  each  Z/L  chart,  read  T[j  for  required 
L/Lp  ratio.  Tabulate  results. 


FOR  CALCULATED  L/Lp 


f.  On  a fourth  chart,  plot  i^  from  Step  5e  versus 
r./L. 

Step  6.  For  required  7/L  ratio,  read  Tk  from  chart  of 
Step  5f.  Calculate  the  unit  blast  impulse  load 
on  the  fojndation: 

<b  " 

Example  A.l:  Computation  of  Averaoe  Impulse  Load  on  Foundation 
Slab  of  Cantilever  Wall  Barrier. 

Required:  Average  impulse  load  on  the  foundation  slab  of  a can- 
tilever wall  barrier. 

Step  1 . Riven; 

a.  Charge  weight:  833.3  lbs 

b.  Structure  dimeiisions  (Figure  A.ia): 

L = 4B  ft  Lp  = 8.5  ft 

c.  Charge  location  parameters  (Figure  A. 14) : 

= 20  ft,  h = 3.5  ft,  I “ 20  ft 
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step  2. 
Step  3. 


Step  4. 


Step  5. 


W - 1.20(833.3)  = 1 ,000  lbs 
Calculate  the  chart  parameters: 

Z/L  = 0.42,  L/Lp  = 5.65,  Lp/h  = 2.43 


R;^/Lp  = 2.35 

Zp  = h/W^/^  = 3.5/(1,000)^''^  = 0.35  ft/lb^^^ 

Interpolation  is  required  for  Zp,  Lp/h, 

L/Lp  and  Z/L. 

Determine  and  tabulate  the  values  of  i^  from 
Figures  A. 2 through  A.  13  for  R/\/Lp  = 2.35, 

Zp  = 0.35  (interpolate  by  inspection)  and 
for  values  given  for  Z/L,  L/Lp  and  Lp/h. 

The  tabulation  of  the  impulses  is  provided  in 
Table  A.l. 

a.  Plot  versus  Lr/h  for  the  values  of 
L/Lp  and  constant  Z/L  (Figure  A. 15). 

b.  Determine  T^j  for  Lp/h  = 2.43,  Z/L  = 0.10 
and  various  L/Lp  ratios  by  entering 
Figure  A. 15(a)  with  Lp/h  = 2.43. 


L/Lp 

'b 

1.0 

105 

3.0 

92 

6.5 

66 

10.0 

55 

c.  Repeat  above  step  for  Z/L  = 0.25  and  0.50 
by  entering  (b)  and  (c)  of  Figure  A.l 5 
with  Lp/h  = 2.43  (tabulation  of  results 
not  shown). 

d.  On  each  Z/L  chart,  plot  Tjj  (Steps  5b  and 
5c)  versus  L/Lp  [upper  abscissa  of  (a) 
through  (c)  of  Figure  A. 15]. 
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tillHli 


e.  Determine  Tu  for  L/Lr  = 5.65  on  each  l/l 
chart  by  entering  (a;  through  Cc)  of 
Figure  A. 15  with  L/Lp  = 5.65  and  reading 
curves  plotted  in  Step  5d. 


Figure 

m 

70 

A. 15(a) 

78 

A. 15(b) 

B 

86 

A. 15(c) 

f.  Plot  i{j  (Step  5e)  versus  l/l  (Figure  A. 16). 

Step  6.  For  l/l  = 0.42  read  Tf^  = 84  ps-ms/lfa^^^  cn 
Figure  A.  16.  Compute  ij^: 

ij,  = = 84(1,000)^/^  = 840  psi-ms 
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APPENDIX  B 


PROCEDURE  FOR  CALCULATING  ARRIVAL  TIME  AND  DURATION  OF 
BLAST  LOADS  ON  THE  STRUCTURE 

B.l  Computation  Procedure  and  Sample  Problem 

This  appendix  contains  the  procedure  (described  in  Sec- 
tion 3.3)  utilized  fay  the  computer  program  to  calculate  the  arrival 
time  and  duration  of  the  blast  loading  on  the  structure.  The  pro- 
cedure is  based  on  the  methods  and  empirical  data  presented  in 
Chapter  4 of  Reference  1 . 

The  computation  proceeds  as  follows: 

Problem:  Determine  the  arrival  and  duration  times  of  blast  loads 
on  structure. 


Procedure : 


Step  1 . Determine  the  following: 

a.  Charge  weight 

b.  Structure  dimensions 

c.  Distance  from  charge  to  back  wall, 

d.  Height  of  charge  above  foundation  slab,  h 

Step  2.  Apply  a 20  percent  safety  factor  to  the 
charge  weight. 

Step  3.  Determine  the  minimum  scaled  distance  from  the 
charge  to  each  loaded  surface.  With  these 
values,  enter  Figure  8.1  and  read  off  the 
curve  labeled  "WWV3'',  the  scaled  times  of 
arrival  of  the  blast  wave  on  each  surface. 

Step  4.  Determine  the  scaled  distances  from  the  charge 
to  the  four  corners  of  each  loaded  surface. 
Enter  Figure  B,1  and  read  off  the  curve 
labeled  "t/\/W‘/3",  the  sc'^^d  tire:  of  arrival 
of  the  blast  wave  at  these  points.  At  the 
same  ti®e.  read  off  the  curve  labeled 
“to/M*'3",  the  scaled  positive  phase  durations 
at  these  points. 


H^EC1':DI^0  SLA^aUiJOT  if'ILMSD 
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Step  5.  Calculate  the  average  time  for  the  wave  to 
fully  engulf  each  surface  and  the  average 
of  the  positive  phase  duration  times  at  the 
farthest  points  on  each  surface. 

4 

C(to)l^vg  “ ^^o^Fl  (^o^F2  ■*■  ^^o^F3  ^^o^F4 

4 

where: 

(t/\)p  = arrival  time  at  one  corner  of 
surface 

(to)p  = positive  phase  duration  at  one 
corner 

^^A^avo  ""  average  time  for  wave  to  fully 
engulf  surface 

(to)ava  ^ average  of  positive  phase  dura- 
tions at  the  farthest  points 

Step  6.  Calculate  the  duration  of  the  loading  on  each 
surface  using  the  following  equation: 

^0  ^'^^A^^avg  ” ^*^A^A  ^ ^^^o^-avg 

where: 

(t^);^  = arrival  time  of  the  blast  wave  at 
the  point  on  the  element  nearest 
to  the  explosion,  defined  by  the 
normal  distance 

Example  B.l:  Cantilever  wall  barrier  with  foundation. 

Required:  Arrival  time  and  duration  of  blast  on  cantilever  wall 
and  foundation. 

Step  1 . Given: 


a.  Charge  weight:  833.3  lbs 

b.  Configuration  of  structure  (Figure  B.2) 
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Figure  B,  1 Shock  wave  parameters  for 


I 

I 

'V 

V' 

step  2 

II 

? Step  3 


Step  4 


c.  = 12  ft 

d.  h = 4 ft 

. W = 1.20(833.3)  = 1,000  lbs 

. Determine  the  minimum  scaled  distance  from 
the  charge  to  each  loaded  surface. 

= 10  15^/3 

Minimum  scaled  distance  to  wall: 

= 12/10  ft(lb)^/^  = 1.2  ft(lb)^/^ 

From  Figure  B.l,  = .1  ms/(lb)^^^ 

(tA)A  = (10)(.l)  = 1.0  ms 

Minimum  scaled  distance  to  foundation; 

Zp  = (4)^  + (4)^  /lO  = .566  ft/(lb)^/3 

From  Figure  B.l  (tA/W^/^)  = .0285  ms/(lb)^^^ 

(tA)A  = 10(.0285)  = .285  ms 

. Determine  Z to  the  four  corners  of  each 
surface. 

Wal  1 : 

= (4)2  + (8)2  + (12)2 

= 1 ,496  ft/(lb)^/^ 

Z3  = (12)2  + (3)2  + (12)2  /lo 

= 1.875  ft/(lb)^/^ 

From  Figure  B.l  read: 

= .155  ms/(lb)l/2 
fl's/(lb)^^^ 
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(ywV3)3  = .23  ms/(lb)l/3 
(to/W^''^)3  = •■'3  ms/(lb)V3 
Foundation: 

= Z2  = 1.496  ft/(lb)l/3 
Z5  = Zg  = V(4)2  + (4)2  + (8)2" /lO 
= .98  ft/(lb)V3 
From  Figure  B.l  read: 

(VWl/3),  = .155  ms/(lb)V3 
(to/wl/3)i  = .090  ms/(lb)l/3 
(tA/W^/3)5  = .073  ms/(lb)^/3 

(to/W^/3)5  = .063  ms/(lb)^/3 

Step  5.  Calculate  and  (to)g^g  for  the  wall 

and  foundation.  ^ 

Wal  1 : 

(tA)^  = 1.5545  ms  (tA)3  = 2.3  ms 

rft«l  1 = 2(1 .55)  + 2(2.3)  = 1.925  ms 


Foundation: 

(tA)i  = 1.55  ms  (tA)s  = .73  i:’? 
C(tA)p]avg  = 2(1.55)  * 2(.73l  = 1 .14  ms 
(tp)^  = .9  ms  {to)^  = .63  ms 

r(t-).-1  = 2(.9)  + 2(.63)  = .765  ms 


.63  ms 


i' 

i'‘ 

t 

i. 


I 

I? 


i 


Step  6.  Calculate  the  durations  (tg)  for  the  wall 
and  foundation. 

Wall: 

fWavg='-®^'“  (V^  = 1.0n,s 

f'‘oUvg"  "" 

(to)w  = 1.925  - 1 .0  + 1 .1  = ?..025  ms 
Foundation: 

f^^^A^F^avg  ^ ^hh  ^ 

■^(^o^F^avg  = ms 

(to)p  = 1.14  - .285  + .765  = 1.62  ms 
Results: 

Wall: 


Arrival  time 

1..000  ms 

Duration  of  Loading 

2.025  ms 

Foundation: 

Arrival  time 

0.285  ms 

Duration  of  Loading 

1 .620  ms 

i 
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Figure  B.2  Example  B.l : Dimensions  of  structure  and  charge 

location  parameters. 


APPENDIX  C 


FOUNDATION  DESIGN  METHOD  FOR 
PROTECTIVE  STRUCTURES  SUSCEPTIBLE  TO  OVERTURNING 

C.l  General 


This  appendix  presents  the  criteria  and  procedure  for 
designing  the  foundations  of  protective  structures  susceptible 
to  overturning.  The  design  procedure  Is  Illustrated  by  two  sample 
probl ems . 


The  design  of  protective  structures  susceptible  to 
overturning  and/or  sliding  motions  proceeds  in  two  stages.  In  the 
first  stage,  the  walls  of  the  structure  are  designed  using  the 
criteria,  methods  and  data  provided  In  Chapters  5 and  6 of  Manual 
TM  5-1300  (Reference  1).  The  methods  In  the  manual  treat  the 
design  of  each  blast-resistant  wall  of  the  structure  Individually. 
The  design  Is  based  on  the  assumption  that  the  supports  along  the 
periphery  of  an  element  are  completely  fixed  against  translation 
and  rotation  and  are  capable  of  fully  developing  the  strength  of 
the  element. 

In  the  second  stage,  the  structure  foundation  Is  de- 
signed for  both  the  conventional  working  (dead  and  live)  load 
and  blast  load  conditions.  Design  procedures  for  the  dead  and 
live  condition  are  well  documented  In  textbooks  on  the  design  of 
reinforced  concrete  structures  (see  Reference  7). 

Designing  the  foundation  for  the  blast  load  condition 
requires  performing  an  overturning  analysis  on  the  structure 
using  the  Overturning  Analysis  Computer  Program  presented  In  Sec- 
tion 5 of  this  report.  The  purpose  of  the  analysis  Is  to  deter- 
mine the  peak  response  of  the  structure  and  the  time  history  of 
the  bearing  pressures  acting  on  the  structure  foundation.  This 
data  is  used  to  determine  the  length  of  the  foundation  extension 
(see  Figures  C.l  and  C.2)  required  to  prevent  the  structure  from 
overturning  or  sliding  (large  distances)  and  the  thickness  of  con- 
crete and  area  of  flexural  reinforcement  required  to  resist  the 
bearing  pressures  developed  in  the  soil  beneath  the  foundation. 

C.2  Preliminary  Tasks 

C.2.1  Introduction 

Two  preliminary  tasks  must  be  accomplished  before  the 
overturning  analysis  of  the  structure  is  performed.  These  tasks 
are:  (1)  estimating  an  initial  size  of  the  foundation  extension 

for  use  in  the  analysis  and  (2)  determining,  with  the  test  data 
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available  for  the  soil  at  the  construction  site,  the  soil  proper- 
ties to  be  used  in  the  analysis. 

Since  some  aspects  of  these  tasks  are  rather  complex, 
a detailed  discussion  of  them  is  provided  before  the  design 
method  is  presented.  These  tasks  are  included  in  the  procedure 
for  designing  the  foundation. 

A discussion  of  the  preliminary  tasks  is  presented  in 
the  following  sections. 

C.2.2  Initial  Estimate  of  Foundation  Size 

In  order  to  perform  an  overturning  analysis,  the  com- 
plete geometry  of  the  structure  must  be  supplied  to  the  computer 
program.  As  described  in  Section  5 of  this  report,  the  informa- 
tion required  consists  of  the  configuration  of  the  structure  des- 
cribed in  terms  of  the  number  and  sizes  (length,  height,  thickness) 
of  each  of  the  blast-resistant  wall  elements  of  the  structure. 

The  sizes  of  the  blast-resistant  walls  are  determined 
using  the  criteria,  methods  and  data  of  Reference  1.  Initially, 
the  size  of  the  foundation  is  not  known  since  the  foundation  de- 
sign, to  some  extent  depending  upon  the  type  of  structure,  is 
governed  by  the  results  of  the  overturning  analysis.  To  proceed 
with  the  analysis,  some  reliable  estimates  of  the  foundation 
dimensions  are  required.  To  this  end,  initial  estimates  of  the 
foundation  dimensions  are  presented  in  the  next  several  paragraphs. 
These  estimates  are  based  on  the  results  of  several  actual  design 
studies  of  protective  structures  susceptible  to  overturning. 

Cantilever  wall  barriers,  lacking  massive  sidewalls 
which  help  prevent  the  structure  from  overturning,  depend  entirely 
on  the  foundation  extension  and  the  soil  beneath  it  to  limit  the 
motions  of  the  structure.  Consequently,  a relatively  large  foun- 
dation extension  is  required.  The  initial  foundation  thickness  to 
be  used  in  the  overturning  analysis  should  be  approximately  1.25 
times  the  wall  thickness  and  tfie  initial  extension  length  utilized 
should  be  approximately  45  percent  of  the  height  of  the  wall  (see 
Figure  C.l).  Sloping  the  bottom  face  of  the  foundation  as  shown 
in  Figure  C.l  is  permissible;  but  the  angle  between  the  bottom 
face  and  the  horizontal  should  be  limited  to  5 degrees.  Larger 
slopes  will  significantly  decrease  the  moment  am  of  the  resultant 
of  the  soil  pressures  about  the  center  of  gravity  of  the  structure. 

Generally,  cantilever  barrier  foundations  should  be 
symmetric  about  the  centerline  of  the  wall  unless  building  con- 
straints dictate  otherwise.  If  the  analysis  indicates  that  the 
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Figure  C.l  Cantilever  wall  barrier  - Estimated  foundation 
dimensions. 
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estimated  length  of  the  extension  is  too  short  to  prevent  over- 
turning, the  foundation  should  be  extended  an  equal  amount  on  both 
sides  of  the  wail.  This  is  more  efficient  than  extending  the 
foundation  outward  from  one  side  of  the  wall  only. 

Single  cell  barriers,  as  a rule,  do  not  require  long 
thick  foundation  extensions  to  prevent  overturning.  The  founda- 
tion thickness  is  established  by  providing  the  foundation  with 
sufficient  bending  strength  (in  both  directions)  to  completely 
develop  the  ultimate  strength  of  each  blast  wall  (backwall  a.nd 
sidewalls,  see  Figure  C.2).  Generally,  this  is  accomplished  by 
making  the  "d^."  of  the  foundation  equal  to  the  maximum  value  of 
“d(.‘‘  for  any  of  the  blast  walls  ("d^."  is  the  distance  between 
centroids  of  the  tension  and  compression  reinforcement  in  an  ele- 
ment, see  Figure  C.2).  With  "dj."  established,  the  areas  of  rein- 
forcing steel  (in  both  directions)  are  computed  (using  Equation 
C.4)  by  equating  the  ultimate  moment  capacities  of  the  foundation 
with  the  moment  capacities  of  the  blast  walls  (see  Figure  C.3). 
Using  the  computed  areas  of  steel,  the  sizes  of  the  actual  rein- 
forcing bars  are  chosen  and  the  thickness  of  concrete  cover  (top 
and  bottom)  is  determined  according  to  the  provisions  of  Section 
7.14  of  the  ACI  Code  (Reference  8).  The  total  thickness  of  the 
foundation  is  the  sum  of  “d^."  (for  the  foundation),  the  diameter 
of  the  outeimiost  reinforcing  bar  and  the  thickness  of  the  top  and 
bottom  concrete  ccver.  The  length  of  the  foundation  extension  is 
established  by  providing  the  anchoiaye  required  for  the  reinforcing 
steel  in  the  concrete.  As  the  foundation  design  is  finalized, 
after  several  overturning  analyses,  the  required  moment  capacity 
of  the  foundation  (and  area  of  reinforcing  steel)  within  the  cell, 
can  be  reduced^  as  illustrated  in  Figure  C.3,  by  subtracting  from 
the  moment  capacity  of  the  wall  the  maximum  moment  developed  in 
the  foundation  extension  during  the  overturning  t'esponse  to  the 
blast  loads. 

In  some  protective  structures,  the  simple  type  founda- 
tion extension  required  to  prevent  overturning  of  the  structure, 
w'ill  be  excessively  long  and  thick.  TMs  usually  occurs  in  the 
following  situations:  (1)  cantilever  walls  with  heights  exceeding 

20  feet  and  (2)  single  cell  barriers  in  which  the  ratio  of  the 
cell  height  to  the  length  of  tine  foundation,  within  the  cell, 
exceeds  1.  Examples  of  these  conditions  are  illustrated  in  Fig- 
ures C.4  and  C.5.  In  these  situations,  backup  stiiictural  elestsents 
(such  as  the  buttress  walls  shown  in  Figures  C.4  and  C.5)  are  re- 
quired to  support  the  foundation.  For  the  pm*pose  of  performing 
the  overturning  analysis,  an  estimate  of  Uie  foundation  thickness 
can  be  made  by  following  the  guidelines  suggested  in  the  preceding 
paragraph  for  single  cell  barriers.  In  this  situation,  however, 
the  moment  capacity  provided  for  the  foundation  is  one  half  of  the 
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Figure  C.2  Single  cell  barrier  - Estimated  Foundation 
dimensions 
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PRELIMINARY  ESTIMATE 


SECTION  THRU  WALL  ELEMENT  AT  FOUNDATION 

FINAL  DESIGN 


SECTION  THRU  WALL  ELEMENT  AT  FOUNDATION 


Figure  C.3  Single  cell  barrier  - Moment  balance  at  junction  of 
backwall  and  foundation  slab 
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ELEVATION  OF  CANTILEVER  WALL 
WITH  SIMPLE  TYPE  FOUNDATION  EXTENSION 


ELEVATION  OF  CANTILEVER  ELEVATION  A- A 

WALL  WITH  BUTTRESS  WALLS 


Figure  C.4  Cantilever  wall  barrier  with  foundation  supported 
by  buttress  walls. 
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nrament  capacity  of  the  backwall  element,  since  the  rigid  supports 
provided  for  the  foundation  extension  (see  Figures  C.4  and  C.5) 
produce  a condition  in  which  the  ultimate  bending  strength  of  the 
foundation  extension  slab  can  be  developed.  Generally,  the  length 
of  the  foundation  extension  required  will  be  approximately  40  per- 
cent of  the  height  of  the  wall. 

Using  the  estimated  dimensions  of  the  foundation  pre- 
sented in  the  preceding  paragraphs,  an  analysis  should  be  performed 
to  determine  the  bearing  pressure  in  the  soil  beneath  the  structure 
for  the  working  load  condition  (dead  and  live  loads).  Decisions 
regarding  the  manner  in  which  these  loads  are  to  be  transferred  to 
the  supporting  soil  (i.e.,  through  foundation  bearing  on  soil  or 
through  piles)  should  be  made  prior  to  performing  the  overturning 
analysis  in  order  to  insure  that  the  finalized  foundation  design 
is  adequate  for  this  load  condition. 

C.2.3  Correlation  of  Soils  Data  and  Overturning 
Design  Criteria 

As  discussed  in  Section  4.3  of  this  report,  the  soil 
data  usually  available  to  the  designer  is  limited  to  the  results 
of  a minimum  of  shallow  test  borings  together  with  a visual  des- 
cription of  the  soil  encountered  and  the  blow  count  from  standard 
penetration  tests.  In  Tables  1 and  2 of  Section  4,  a correlation 
of  soil  properties  with  the  above  data  was  presented.  This  cor- 
relation provides  the  designer  with  the  means  of  determining  the 
properties  of  the  soil  (at  the  construction  site)  for  use  in  the 
overturning  analysis. 

The  properties  of  the  soil  grossly  affect  the  response 
of  the  structure.  Therefore,  these  properties  have  a large  impact 
on  the  foundation  design.  Tables  1 and  2 of  Section  4 provide 
for  a particular  soil,  the  properties  in  the  soft  or  loose  condi- 
tion and  the  compact  or  hard  condition.  The  actual  condition  of 
the  soil  at  a given  site  will  be  somewhere  between  these  extremes. 

In  order  to  account  for  the  most  severe  design  condi- 
tions for  both  overturning  and  strength,  the  structure  is  analyzed 
for  both  conditions  of  the  soil.  The  response  of  the  structure  on 
the  soft  condition  of  the  soil  establishes  the  length  of  the  foun- 
dation extension  required  to  prevent  overturning,  whereas  the  res- 
ponse of  the  compact  condition  of  the  soil  dictates  the  thickness 
of  concrete  and  amount  of  reinforcing  steel  required  for  the  foun- 
dation extension  to  resist  the  bearing  pressures  developed  in  the 
soil  beneath  the  structure. 
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EXPLOSIVE 


Generally,  rotations  of  the  structure  that  approach 
the  point  of  incipient  overturning  (as  defined  in  Figure  C.6)  can 
be  tolerated.  Therefore,  for  the  design  to  be  efficient,  the  peak 
response  of  the  structure  on  the  soft  soil  should  approach  in- 
cipient overturning.  To  insure  that  the  structure  will  not  over- 
turn on  the  soft  soil,  the  peak  rotation  of  the  structure  on  the 
compact  soil  is  limited  to  a percentage  of  the  overturning  angle 
(defined  in  Figure  C.7).  The  results  of  sev;)ral  design  studies 
indicate  that  limiting  the  peak  rotation  of  the  structure  on  the 
compact  soil  to  a value  of  approximately  40  percent  of  the  over- 
turning angle  will  insure  that  the  structure  will  not  overturn, 
but  will  approach  a peak  response  of  incipient  overturning  on  the 
soft  soil . 


The  guidelines  presented  in  the  previous  paragraph  are 
utilized  in  the  design  of  the  foundation  extension  by  performing 
a series  of  overturning  analyses.  After  each  analysis,  the  dimen- 
sions of  the  foundation  extension  ara  modified,  according  to  the 
analysis  results.  This  process  is  repeated  until  the  results  of 
the  analyses,  for  both  soil  condit'ons,  indicate  that  the  structure 
rotates  to  40  percent  of  its  overturning  angle  on  the  compact  soil 
and  does  not  overturn  on  the  soft  soil.  This  procedure  is  generally 
applicable  to  cantilever  wall  barriers  only.  However,  it  can  be 
applied  to  single  cell  barriers  provided  it  does  not  alter  the 
foundation  dimensions  to  the  extent  that  the  following  minima  are 
not  nwintained:  C)  the  minimum  length  of  the  foundation  extension 

required  for  anchorage  of  the  reinforcement  in  the  concrete  or  (2) 
the  minimum  required  plan  size  to  transfer  the  conventional  (dead 
and  live)  working  loads  to  the  supporting  soil  without  exceeding 
the  allowable  bearing  pressure  for  the  soil.  In  the  event  that 
piles  are  utilized,  Ite'  (2)  can  be  ignored. 

The  criteria  presented  in  this  section  and  the  data 
presented  in  Section  4.3  are  intended  to  be  used  when  more  reliable 
soil  data  are  not  available.  If  more  reliable  data  are  available, 
the  structure  should  be  analyzed  for  the  specific  properties  de- 
rived from  the  data.  In  this  situation,  the  structure  can  be  per- 
mitted to  lotate  to  the  point  of  incipient  overturning  under  the 
action  of  the  blast. 

C.3  Design  Criteria  for  Foundation  Extensions 

C.3.1  Introduction 


As  the  structure  responds  to  the  blast  loads,  high 
bearing  pressures  are  developed  in  the  soil  beneath  the  foundation 
extension.  Tfiese  pressures  impart  shears  and  bending  moments  in 
the  foundation  extension  for  which  the  foundation  slab  must  be 


115 


designed.  Since  shear  reinforcement  is  not  utilized,  the  thick- 
ness of  the  member  will  generally  be  dictated  by  the  shear  design. 
Basically,  the  design  involves  determining  the  thickness  of  con- 
crete required  to  resist  the  peak  shear  load  and  the  area  of  re- 
inforcing steel  needed  to  resist  the  maximum  bending  moments. 

There  are  two  types  of  foundation  extensions  utilized 
with  protective  structures.  Most  structures  require  only  the 
simple  type  extension  shown  in  Figures  C.l  and  C.2.  The  extension 
behaves  essentially  as  a cantilever  beam.  Generally,  a short  deep 
member  will  be  required  to  resist  the  applied  shears.  In  most 
cases,  the  design  of  the  simple  type  extension  will  fall  under  the 
deep  beam  provisions  of  Section  11.9  of  the  ACI  Code  (Reference  8). 

Tall  barriers  will  generally  require  a more  complex 
foundation  system  such  as  those  shown  in  Figures  C.4  and  C.5.  In 
these  situations,  the  foundation  behaves  as  a two-way  member  with 
a portion  of  the  applied  soil  bearing  pressures  reacted  by  the 
buttress  walls.  The  two-way  action  will  generally  result  in  a 
thinner  foundation  slab  which  is  designed  using  the  methods  and 
data  of  Section  5,  Reference  1. 

In  the  following  sections,  the  equations  for  the  allow- 
able concrete  shear  stresses  are  presented  for  both  types  of  foun- 
dation extensions.  Presented  also  are  the  equations  for  computing 
the  ultimate  unit  resisting  moments  of  the  foundation.  In  addi- 
tion, minimum  areas  of  flexural  reinforcement  are  specified. 

C.3.2  Thick  Foundation  Extension 

The  term  "thick  foundation  extension"  applies  to  mem- 
bers having  a 1^/^,  ratio  of  less  than  5 where  the  term  "1^"  is 
the  clear  span  of  the  member  to  the  face  of  the  support  (see 
Figure  C.8)  and  "d"  is  the  distance  from  the  extreme  compression 
fiber  to  the  centroid  of  the  tension  reinforcement.  Generally, 
tlie  simple  type  foundation  extension  will  usually  fall  into  this 
category. 


The  allowable  shear  stress  carried  by  the  concrete  is 
computed  according  to  the  following  equation: 

Vc  = ^[3.5  - (2.5Hc^)/Vud][1.9v'fJ  + 2,500f\y(Vud/Mj.^)] 

(C.l) 

except  that  the  tenii  "[3.5  - (2.5M(.y.)/Vyd]"  shall  not  exceed  6vT^. 
The  terms  in  the  equation  are  defined  on  tJie  following  page. 
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4)  = capacity  reduction  factor 

= 0.85  for  all  sections 

Vq  = nominal  permissible  shear  stress  carried  by 
the  concrete,  psi 

= applied  design  load  moment  at  the  critical 
section,  in-lb 

= total  applied  design  shear  force  at  critical 
section,  lb 

d = distance  from  extreme  compression  fiber  to 
centroid  of  tension  reinforcement,  in 

f^  = specified  compressive  strength  of  concrete,  psi 

Py,  = ratio  of  area  of  flexural  reinforcement  to  area 
of  concrete  within  depth  "d"  and  width  "b" 

= Ag/bd 

b = width  of  the  compression  face,  in. 

This  equation,  with  the  exception  of  the  4>  factor,  is 
presented  in  Section  11-9  (Equation  11-22)  of  the  ACI  Code.  To  be 
consistent  with  the  convention  utilized  in  Reference  1,  the  capac- 
ity reduction  factor  is  applied  to  the  allowable  stress  instead  of 
to  the  applied  stress  (which  is  the  method  of  Reference  8). 

The  applicability  of  Equation  (11-22)  is  restricted  by 
the  ACI  Code  to  deep  beams  (Ip/d  < 5)  loaded  at  the  top  or  com- 
pression face  (as  in  simply  supported  beams).  In  this  report,  the 
applicability  of  this  equation  is  extended  to  the  simple  type 
foundation  extensions  (of  protective  structures)  which  have  an 
Ip/d  ratio  of  less  than  5.  These  members  are  cantilevered  from 
tne  backwall  elements  and  hence,  loaded  at  their  tension  faces. 
However,  the  same  load-carrying  mechanism  (arch  action)  Is  devel- 
oped In  the  deep  cantilevered  foundation  extension  as  exists  in 
the  deep  beam  loaded  at  its  compression  face  (see  Reference  9). 
Therefore,  until  further  test  data  indicate  otherwise,  it  is  con- 
sidered appropriate,  in  this  instance,  to  use  this  equation  for 
determining  the  permissible  shear  stress  carried  by  the  concrete. 

The  critical  section  for  shear  Is  taken  as  15  percent 
of  the  clear  span  (0.15  l^}  of  the  member  measured  from  the  face 
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of  the  support.  This  value  is  specified  in  Section  11.9.3  of  Ref- 
erence 8 and  it  applies  to  deep  beams  loaded  by  uniformly  distrib- 
uted loads. 


C.3.3  Thin  Foundation  Extension 

In  situations  in  which  the  clear  span  of  the  foundation 
extension  is  greater  than  5 times  the  "d"  distance,  the  provisions 
of  Section  5.3  of  Reference  1 are  utilized  to  determine  the  per- 
missible shear  stress  carried  by  the  concrete.  Generally,  founda- 
tion slabs  which  behave  as  two-way  members  will  fall  into  this 
category. 


The  permissible  shear  stress  carried  by  the  concrete  is 
computed  using  Equation  5-10  of  the  reference.  This  equation  is 
provided  below  (Equation  C.2).  The  parameters  in  the  equation  are 
defined  in  Section  C.3.2. 

V(.  = 4>(1.9/fJ  + 2.500p^J  < 2.28^/fJ  (C.2) 

In  a thin  member,  the  critical  section  for  shear  is  taken  at  a 
distance  of  "d"  from  the  face  of  the  support. 

C.3.4  Ultimate  Resisting  Moment 

When  designing  the  foundation  extension  to  resist  the 
build-up  of  soil  bearing  pressures  beneath  the  structure,  the 
ultimate  unit  resisting  moment  of  the  member  is  computed  utilizing 
the  following  equation: 

”u  " (As^^s^td  - (a/2)]/b  (C.3) 

wherein: 


= ultimate  unit  resisting  moment,  in-lb/in 

Aj  = area  of  tension  reinforcement  within 
the  width  b,  in2 

fj  ' static  design  stress  for  reinforcement,  psi 

d - distance  from  extreme  compression  fiber 
to  tension  reinforcement,  in 

a = depth  of  equivalent  rectangular  stress 
block,  in 


b = width  of  compression  face,  in 

f(l  = specified  compressive  strength  of 
concrete,  psi. 

fienerally,  in  blast  design,  an  equal  amount  of  flexural  reinforce- 
ment is  provided  at  both  faces  of  the  member.  In  most  cases, 
however,  the  ultimate  moment  capacity  of  the  member  can  be  ap- 
proximated within  acceptable  limits,  by  disregarding  the  rein- 
forcement in  the  compression  face.  Equation  C.3  reflects  this 
approximation. 

When  designing  the  foundation  to  develop  the  strength 
of  a blast  wall  element,  the  ultimate  unit  resisting  moment  of  the 
member  is  computed  using  the  following  equation: 

^ = \^ds"^c/^  (C.4) 

The  parameters  in  this  equation  are  defined  in  the  preceding  para- 
graph with  the  exception  of  tne  quantities  "f^js"  and  "d^"  which 
are  defined  below: 

f^5  = dynamic  design  strength  for  the 
reinforcement  which  is  determined 
according  to  the  provisions  of 
Section  5-6  of  Reference  1 , psi 

d(.  = distance  between  the  centroids  of  the 

compression  and  tension  reinforcement,  in. 

C.3.5  Minimum  Flexural  Reinforcement 

In  order  to  insure  the  proper  structural  behavior  and 
also  to  prevent  excessive  cracking  and  deformations  under  conven- 
tional loadings,  the  minimum  areas  of  flexural  reinforcement 
listed  in  Table  C.l  are  recommended.  The  areas  specified  are  the 
minimum  areas  of  reinforcement  at  the  tension  face  of  the  member. 
An  equal  amount  of  steel  should  be  provided  at  the  compression 
face  also.  The  quantities  specified  in  tl\e  table  were  extracted 
from  Table  5-1  of  Reference  1. 

Table  C.l  is  presented  on  the  following  page. 
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TABLE  CJ 


MINIMUM  AREA  OF  FLEXURAL  REINFORCEMENT 


One-Way 

Two-Way 

Reinforcement 

Elements 

Elements 

Main 

A^  = 0.0025bd 

As  = 0.0025bd 

Other 

Aj  - 0.001  Ob T(. 

Ag  = O.OOlSbd 

The  parameters  in  the  table  are  defined  in  Section  C.3.3  with 
the  exception  of  "T^."  which  is  defined  as  the  total  thickness  of 
the  member. 


C.4  Foundation  Design  Procedure 

C.4.1  Introduction 

The  method  outlined  below  utilizes  the  overturning 
analysis  computer  program  to  compute  the  response  of  the  structure 
on  its  supporting  soil.  The  results  of  the  computerized  analysis 
(the  time  history  of  the  soil  bearing  p»*essures)  are  utilized  to 
design  the  foundation  extension  of  the  structure. 

The  presentation  of  the  design  procedure  is  divided  into 
two  parts.  In  the  first  part  (Steps  1 through  9),  the  emphasis  is 
on  establishing  the  length  of  the  foundation  extension  required  to 
prevent  overturning  and/or  sliding  of  the  structure.  This  is  ac- 
complished by  a series  of  overturning  analyses  in  which  the  lenoth 
of  the  foundation  extension  is  either  increased  or  decreased  de- 
pending on  the  results  of  the  previous  analysis.  The  overturning 
analyses  are  continued  until  a satisfactory  result  is  achieved. 

The  second  part  of  the  procedure  (beginning  with  Step  10),  deals 
with  the  methods  for  designing  the  foundation  extension  to  resist 
the  peak  hearing  nres.sures  in  the  soil.  As  the  overall  proceduv-e 
is  intended  to  be  applicable  to  structures  with  either  the  simple 
(cantilevered)  or  the  more  complex  type  (two-way  member  su(hx)rted 
on  three  or  four  sides)  foundation  extension,  the  second  part  of 
the  procedure  is  divided  into  two  separate  presentations.  One 
presentation,  consisting  of  Steps  10a  through  15a,  deals  exclu- 
sively with  the  design  of  the  simple  type  (cantilevered)  extension. 
A separate  presentation,  consisting  of  Steps  10b  t.r''';cli  17b,  is 
provided  which  treats  the  design  of  the  more  comjjlex  (two-way  mem- 
ber) type  foundation  extension.  The  material  is  presented  in  this 
manner  because  the  procedures  for  designing  both  types  of  founda- 
tion extensions  follow  different  paths  to  a final  solution.  The 


desitjn  of  the  simple  type  extension  proceeos  in  a straightfonvard 
manner  in  which  the  member  is  proportioned  and  the  reinforcing 
steel  is  supplied  to  resist  the  applied  shears  and  bending  moments, 
which  are  calculated  using  the  soil  bearing  pressure-time  history 
from  the  computerized  overturning  analysis. 

The  tv/o-WoV  member  is  designed  in  an  indirect  manner. 
Briefly,  the  procedure  involves  determining  the  amount  of  flexural 
reinforcement  required  for  the  foundation  extension  to  have  an 
ultimate  resistance  in  bending  which  is  equal  to  or  greater  than 
the  peak  design  load.  The  ultimate  resistance  of  the  element  in 
bending  is  computed  using  yield  line  theory.  After  the  airount  of 
reinforcement  is  determined,  the  element  is  checked  for  shear  and, 
if  necessary,  its  thickness  is  increased  until  the  ultimate  re- 
sistance of  the  element  in  shear  equals  or  exceeds  the  applied  load. 
The  applied  shears  and  bending  moments  are  not  utilized  to  design 
the  member.  A direct  computation  of  these  quantities  is,  for  a 
two-way  element,  an  extremely  complex  '■^sk  requiring  the  utiliza- 
tion of  elastic  plate  bending  theory  and  numerical  methods  (such 
as  the  finite  element  technique). 

V 

C.4.2  Analysis 

The  procedure  contnenc'S  after  the  design  of  the  blast- 
resistant  wall  elements  has  been  completed.  /»t  this  point,  the 
designer  has  available  the  follow  •'q  data*  (1)  tf  ^ configuration 
of  the  structure  that  is  required  , confi.ne  explosion,  {?,) 
the  sizes  and  design-  details  of  all  the  ‘^last-resistant  wall 
elements  (includinvg  thicKness,  amount  of  reinfo»'cinq  steel)  and 
(3)  the  average  unit  impulse  loads  on  the  wall  elements.  From 
this  point,  the  analysis  proceeds  as  follows: 

Probl&'ii:  Design  the  foundation  extension  of  a protective  struc- 
ture susceptibU  to  overtuniing  and/or  sliding. 


Procedure : 


Step  1 . Fstablish  the  design  parameters: 

a.  C<  nfiguration  of  structure  and  design 
details  of  blast  wall  elements 

U.  Quantities  and  locations  of  explosives 

c.  Available  test  and  descriptive  data  for 
soil  at  construction  site 
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d.  Design  strength  of  concrete  and  rein- 
forcing steel. 

Step  L.  Based  on  the  configuration  of  the  structure 
and  the  guidelines  of  Section  C.2.2,  estimate 
the  dimensions  of  the  foundation  evtension  to 
be  utilized  in  the  overturning  analysis.  For 
single  cell  barriers  of'  barriers  supported  by 
buttress  walls,  determine  the  area  of  rein- 
forcement required  for  tlie  foundation  to 
develop  the  strength  of  the  backwall  and  side- 
walls  of  the  structure. 

Step  3.  Determine  the  soil  bearing  pressures  beneath 

the  foundation  (using  the  foundatioc  dimensions 
estimated  in  Step  2)  for  the  working  (dead  and 
live)  load  condition. 

The  foundation  must  have  sufficient  plan  size 
to  transfe'-  the  dead  and  live  loads  to  the  sup- 
porting soil  without  exceeding  the  allowable 
bearing  pressure  for  the  soil.  If  the  allow- 
able bearing  pressure  is  exceeded,  the  length 
and,  where  feasible,  the  width  of  tne  founda- 
tion Should  be  increased  in  order  to  provide 
the  plan  size  reotiired.  If  an  excessively 
large  plan  site  is  required,  piles  should  be 
iiseo.  In  any  event,  the  plan  size  should  not 
be  decreasfsi  unless  the  results  of  the  subse- 
quent over turn inq  analysis  indicate  othe»vise. 

Step  d.  Apply  a ?('  oercort  safety  factor  to  the  charge 
weights  and  determine  tne  average  unit  is^lse 
loads  on  tl>e  •’■-niRdation  slab  within  t«?e  cell. 

To  deterrn'ne  the  loads  on  the  foundations  of 
cantilever  wcill  barr'ers,  utili.’e  the  impulse 
charts  and  thg  i nterpoldtiou  i-ethod  of  Appendix 
.A.  For  sinqip  tsss  barriers,  use  the  ."methods 
arhi  -fata  of  Ref'-’-etKC  1 or  the  csssputer  s'rograsr 
of  Reference  ? to  dr-fer-ine  th-  impulse  loads 
on  the  founaation  s'suh.  If  the  structure  coo= 
figurdtion  does  no?  tonfor*?^  to  any  of  those 
shOivn  in  Tigure  ii,  the  arrival  lirte  and  dura- 
tion of  the  Mast  loads  on  each  loade?!  surface, 
produced  by  each  dosive  = ' arge  H-ave  to  be 
cor^puted  by  hand  uti lining  the  procedure  and 
dat/J  oresenterf  in  Ar-penais  B. 


Si. 


step  5.  Correlate  the  available  test  and  descriptive 
data  for  the  soil  at  the  construction  site 
with  the  data  of  Tables  1 and  2 of  Section  4.3 
and  establish  the  range  of  critical  soil  prop- 
erties to  be  utilized  in  the  analysis.  The 
structure  is  analyzed  for  both  the  soft  and 
compact  conditions  (specified  in  the  tables) 
of  the  actual  soil.  In  the  event  that  more 
accurate  data  are  available,  only  the  actual 
soil  condition  need  be  considered  in  the 
analysis. 

Step  6.  Prepare  the  input  data  for  the  computer  program 
according  to  the  instructions  of  Section  5. 
Generally,  two  data  decks  are  required,  each 
containing  the  properties  (determined  in  the 
previous  step)  for  one  condition  of  the  soil. 

Step  7.  Run  the  two  analyses  (if  two  are  required) 
utilizing  the  overturning  analysis  computer 
program. 

Step  8.  Inspect  the  results  of  both  analyses  and  de- 
termine for  each: 


a.  If  the  structure  reached  its  peak  response: 

In  the  event  that  the  structure  failed  to 
attain  its  peak  response,  rerun  the  analysis 
utilizing  more  integration  time  steps. 

b.  If  the  structure  overturned: 


If  the  structure  overturns,  it  will  do  so 
on  the  soft  condition  of  the  soil.  If 
this  occ'-'s,  the  length  of  the  foundation 
extension  will  have  to  be  increased  and 
both  analyses  rerun. 

c.  If  the  structure  experienced  excessive 
horizontal  (sliding)  displacements  under 
the  action  of  the  blast: 


The  horizontal  displacements  become  a 
factor  when  explosives  are  stored  nearoy. 
In  this  situation,  there  is  danger  of  the 
structure  s’’iding  into  the  explosives  and 
detonating  them,  thereby  propagating  the 
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explosion.  Generally,  large  sliding 
motions  will  occur  on  the  cohesive  (clay) 
soils.  This  condition  is  remedied  by 
adding  mass  to  the  structure  foundation 
in  order  to  increase  the  friction  forces 
between  the  structure  and  the  soil.  The 
added  mass  will  also  lower  the  center  of 
gravity  of  the  structure  which  causes  the 
toe  of  the  foundation  to  penetrate  further 
into  the  soil  thereby  decreasing  the 
sliding  motions.  However,  this  will  in- 
crease the  rotations  of  the  stv'-ucture  and 
therefore  the  foundation  extension  may 
have  to  be  lengthened,  depending  on  the 
results  of  the  previous  analysis.  The 
revised  structure  is  always  re-analyzed 
for  both  conditions  of  the  soil. 

Further  evaluation  of  the  results  is  de- 
ferred until  both  analyses  indicate  that 
the  structure  attained  its  peak  response 
and  will  neither  overturn  nor  translate 
large  distances  under  the  action  of  the 
blast. 

Step  9.  Once  the  conditions  of  Step  8 are  satisfied, 
evaluate  the  results  of  the  analysis  in  order 
to  determine  if  any  further  modification  of 
the  foundation  dimensions  is  required.  These 
added  modifications  are  required  when  the 
analysis  indicates  that  the  structure  attains 
peak  rotations  which  are  far  less  than  what 
is  permitted.  In  most  cases,  analyses  of  the 
structure  supported  by  both  a soft  and  a com- 
pact soil  are  required  and  generally,  the 
foundation  extension  can  be  shortened  until 
the  results  of  the  analyses  indicate  that  the 
structure  rotates  to  40  percent  of  its  over- 
turning angle  on  the  compact  soil,  provided 
the  decreased  length  of  the  foundation  e.xten- 
sion  dees  not  violate  the  design  criteria  in 
Section  C.2.3. 

In  the  event  that  more  accurate  soils  data  are 
available,  the  structure  can  be  allowed  to 
rotate  to  the  point  of  incipient  overturning 
(see  Figures  C.6  and  C.7). 
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C.4.3  Design:  Simple  Type  Foundation  Extension 

Once  satisfactory  results  are  obtained  from  the  analysis, 
the  detailed  design  of  the  foundation  extension  commences.  The 
procedures  utilize  the  bearing  pressure  time  history  printed  out 
by  the  computer  program.  In  most  cases,  analyses  are  perfot^ed  for 
the  structure  supported  cn  both  the  soft  and  compact  conditions  of 
the  actual  soil.  The  foundation  extension  is  designed  for  the 
bearing  pressures  developed  on  the  compact  soil. 

This  section  treats  the  design  of  simple  type  (canti- 
levered) foundation  extensions  of  the  type  normally  utilized  with 
cantilever  wall  barriers  and  single  cell  barriers.  In  this  sec- 
tion, the  instructions  (steps)  are  labeled  with  the  suffix  "a" 

(i.e..  Step  10a.)  to  distinguish  them  from  the  instructions  pre- 
sented in  the  next  section. 

The  computation  proceeds  as  follows: 

Step  IQa.  Determine  tne  location  of  the  critical 

section  for  shear  according  to  the  provi- 
sions of  Section  C.3. 

Step  11a.  Determine  from  the  printout  of  the  soil 
bearing  pressure  time  history,  the  peak 
shear  (and  corresponding  bending  moment  for 
thick  sections,  l^/d  < 5)  at  the  critical 
section  for  sheav*  and  the  peak  bending 
moirent  at  tiae  face  of  the  support. 

These  quantities  (Vy,  and  My)  are  com- 
puted by  the  program  for  cantilever  wall 
barriers.  For  other  structural  configura- 
tions, these  quantities  must  be  computed 
manually. 

Generally,  the  bearing  pressure  distributions 
at  several  tin<e  stations  are  investigated  to 
determine  the  peak  shear  and  bending  moments. 
When  the  computation  is  performed  by  the 
cOft^uter  program,  the  bearing  pressure  dis- 
tribution at  every  integration  time  station 
is  investigated. 

The  peak  shear  usually  occurs  when  the  point 
on  the  foundation  with  ze*'o  bearing  pressure 
approaches  the  critical  section  for  shear  on 
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the  extension.  Figure  C.8  shows  the  config- 
uration of  the  bearing  pressure  distribution 
curve  and  identifies  the  design  parameters 
and  critical  sections  for  shear  and  bending. 
The  figure  also  illustrates  the  order  in 
which  the  bearing  pressures  are  printed  out 
by  the  computer  program. 

Figure  C.9  shows  the  free  body  diagrams  for 
computing  the  peak  shear  and  bending  moment 
on  the  foundation  extension.  The  shear  is 
determined  by  computing  the  area  within  the 
bearing  pressure  distribution  curve.  The 
bending  moment  is  determined  by  computing 
the  moment  of  the  area  within  the  bearing 
pressure  distribution  curve  about  the  de- 
sired location  on  the  extension. 

Step  12a.  Determine  the  allowable  shear  stress  that 
can  be  carried  by  the  concrete  using  the 
equations  provided  in  Section  C.3. 

Step  13a.  Determine  the  thickness  of  concrete  ("d") 
required  to  carry  the  peak  applied  shear 
load, 

d = V,/v, 

where : 

d = distance  from  extreme  compression 
fiber  to  tension  steel,  in 

Vy  = peak  applied  shear  load,  lb/ in 

Vj.  = permissible  shear  stress  carried 
by  concrete,  ,^si 

Step  14a.  Determine  the  area  of  the  flexural  (main) 

reinforcement  required  using  the  peak  applied 
bending  moment  (at  the  face  of  the  support) 
in  Equation  (C.3).  At  ttie  same  time,  compute 
the  aroa  of  reinforcement,  perpendicular  to 
the  main  reinforcement,  according  to  the 
provisions  of  Table  C.l. 
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ELEVATION 


notes: 

1.  NS  DENOTES  NUMBER  OF  SOIL  ELEMENTS 

2.  THE  NUMBERS  ON  THE  BOTTOM  FACE  OF  THE  FOUNDATION 
CORRESPOND  TO  THE  NUMBERS  PRINTED  OUT  BY  THE 
COMPUTER  PROGRAM  AT  THE  TOP  OF  EACH  PAGE  CONTAIN- 
ING THE  SOIL  BEARING  PRESSURE-TIME  HISTORY. 


Figure  C.8  Design  paraineters  - Simple  type  foundation 
extension. 
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Mu  (PEAK  BENDING  MOMENT) 


<n-fcr 


FREE  BODY  DIAGRAM  FOR  SHEAR  COMPUTATION 


FREE  BODY  DIAGRAM  FOR  BENDING  MOMENT  COMPUTATION 

Figure  C.9  Free  body  diagrams  of  simple  type  foundation 
extension  for  computation  of  peak  shear  and 
bending  moirent, 
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step  15a.  Determine  the  actual  thickness  of  the  foun- 
dation using  the  following  equation; 

Tc  = d + db/2  + c (C.5) 

where: 

= thickness  of  foundation,  in 

d = distance  from  extreme  compression 
fiber  to  centroid  of  tension 
reinforcement,  in 

dfj  = diameter  of  tension  reinforcement 
bar,  in 

c = thickness  of  bottom  concrete 
cover  specified  in  Section  7.14 
of  Reference  8 (always  3 inches) 

Depending  on  the  configuration  of  the  struc- 
ture, a significant  decrease  in  the  founda- 
tion thickness  could  result  in  a substantial 
increase  in  the  peak  rotation  of  the  struc- 
ture. This  is  generally  the  case  with  can- 
tilever wall  barriers;  therefore,  if  the 
foundation  thickness  computed  is  much  less 
than  the  thickness  used  in  the  overturning 
analysis,  repeat  the  analysis  with  the  re- 
vised foundation  thickness  in  order  to 
verify  the  final  design.  The  verification 
analysis  is  generally  not  required  for 
single  cell  barriers. 

This  completes  the  design  of  the  simple  type  foundation 
extension. 

C.4.4  Design:  Two-Way  Foundation  Extensions 

This  section  treats  the  design  of  foundation  extensions 
which  behave  as  bvo-way  members.  In  this  section,  tlie  instructions 
(steps)  are  labeled  with  the  suffix  "b"  (i.e.,  Step  10b.)  to  dis- 
tinguish them  from  the  instructions  presented  in  the  previous 
section. 


The  computation  commences  upon  t)ie  satisfactory 
completion  of  the  overturning  analyses  (Step  9). 
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FNDN.BEAM  SUPPORTING 
FNDN.  SLAB  — r 


FNDN.  EXTENSION  SUPPORTED  ON  FOUR  SIDES 


FNDN.  EXTENSION  SUPPORTED  ON  THREE  SIDES 


Figure  C.IO  Plan  views  of  foundation  extensions  supported 
on  three  and  four  sides. 


The  computation  proceeds  as  follows; 

Step  10b.  Determine,  from  the  output  of  the  computer 
program,  the  design  loading  for  the  foun- 
dation extension.  The  design  loading  is 
defined  as  the  soil  bearing  pressure  dis- 
tribution corresponding  to  the  maximum  total 
load  on  the  portion  of  the  surface  area  of 
the  foundation  extension  which  is  beyond  the 
critical  section  for  shear  ("d"  from  support) 
as  shown  in  Figure  C.IO. 

Step  lib.  Establish  the  initial  amount  of  flexural 

reinforcement  (for  the  foundation  extension) 
to  be  used  as  a starting  point  in  the  design 
computations.  The  thickness  of  the  founda- 
tion and  the  amount  of  flexural  reinforce- 
ment in  the  direction  perpendicular  to  the 
backwall  element  (see  Figure  C.ll)  are  de- 
termined in  Step  2 on  the  basis  of  providing 
the  ultimate  moment  capacity  required  to  de- 
velop the  strength  of  the  blast  walls.  For 
the  reinforcement  parallel  to  the  backwall 
element,  utilize  the  minimum  area  of  rein- 
forcement specified  in  Table  C.l.  The  same 
;-mount  of  reinforcement  should  be  supplied 
at  both  faces  of  the  member.  If  the  area 
of  reinforcement  perpendicular  to  the  back- 
vyall  is  less  than  the  minimum  quantity  spec- 
ified for  the  main  reinforcement,  the  area 
of  reinforcement  parallel  to  the  backwall 
must  be  greater  than  the  specified  minimum 
for  the  main  reinforcement. 

Step  12b.  Compute  the  ultimate  positive  and  negative 
bending  moment  capacities  in  both  directions 
for  the  foundation  extension  using  Equation 

C.3. 

Step  13b.  Utilizing  yield  line  theory,  compute  the 

ultimate  resistance  in  bending  of  the  foun- 
dation extension.  The  computation  is  per- 
formed for  a load  distribution  which  has  the 
same  characteristics  as  the  design  load 
acting  on  tiie  foundation.  Section  5-10  of 
Reference  1 presents  a discussion  of  yield 
line  theory  and  illustrates  the  computation 
of  the  ultimate  resistance  for  a member  sub- 
jected to  a uniform  load.  The  same  methods 
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PARTIAL  PLAN  OF  FNDN. 


FLEXURAL  REINF. 

TO  DEVELOP  BACK- 
WALL  AND  RESIST 
SOIL  BEARING 
PRESSURES 


FLEXURAL  REINF 
TO  RESIST  SOIL 
BEARING  PRESaiRES 
( ASSUME  MIM  TO 
START  ) 


SECTION  A-A 

Figure  C.ll  Single  cell  barrier  with  buttress  walls  - 
Moment  balance  at  junction  of  backwall  and 
foundation  slab. 


are  utilized  in  this  problem  for  a member 
loaded  by  a trapezoidal  load  distribution. 

Step  14b.  Compare  the  computed  resistance  with  the 
actual  design  load.  If  the  computed  re- 
sistance is  greater  than  the  applied  load, 
go  on  to  the  shear  design  (Steps  15b  and  16b). 

If  the  computed  resistance  is  less  than  the 
applied  load,  increase  the  area  of  reinforce- 
ment in  the  direction  parallel  to  the  back- 
wall  element  and  recompute  the  ultimate  re- 
sistance. Repeat  this  computation  until  the 
computed  resistance  of  the  foundation  exten- 
sion equals  or  exceeds  the  applied  load. 

Step  15b.  Based  on  the  final  location  of  the  yield 
lines  (computed  in  Step  14b.),  compute  the 
shear  stress,  produced  by  the  design  load, 
at  the  critical  section  for  shear  for  each 
sector  of  the  foundation  (see  Figure  C.12). 

The  design  shear  load  for  a sector  is  simply 
the  applied  load  on  the  portion  of  the  sector 
beyond  the  critical  section  for  shear.  The 
applied  stress  is  computed  by  dividing  the 
design  shear  load  for  the  sector  by  the 
product  of  the  width  of  the  sector  at  the 
critical  section  and  the  depth  of  concrete 
("d"). 

A discussion  and  illustration  of  the  shear 
stress  computation  is  presented  in  Section 
5-20  of  Reference  1 . 

Step  16b.  Compute  the  allowable  shear  stress  for  the 
concrete  using  Equation  C.2.  Compare  the 
allowable  shear  stress  with  the  applied 
shear  stress. 

K the  applied  shear  stress  is  less  than  the 
allowable  shear  stress,  pt*oceed  to  the  next 
step. 

If  the  applied  shear  stress  is  greater  than 
the  allowable  shear  stress,  increase  the 
thickness  of  the  foundation  and  proportion 
the  areas  of  the  reinforcement  such  that  the 
revised  foundation  has  >‘oughly  the  same 
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C,5 

Example  C. 
Required: 


bending  moment  capacities  as  the  original 
foundation.  If,  in  oroportioninq  the  areas 
of  reinforcement,  the  area  of  the  steel  in 
one  or  both  directions  falls  below  the  minima 
specified  in  Table  C.l,  use  the  minimum  areas. 

Recompute  the  following  for  the  revised 
foundation: 

1.  Yield  line  location 

2.  Ultimate  bending  resistance  of  the 
member 

3.  Shear  stress  at  the  critical  section 
for  shear  at  both  sectors. 

Repeat  the  above  computation  cycle  until 
the  results  indicate  that  the  applied  shear 
stress  is  less  than  the  allowable  shear 
stress  for  the  member. 

Step  17b.  Compute  for  all  sectors  the  shear  load  at  the 
face  of  the  supports.  Using  these  shear 
loads,  design  the  elements  Uuch  as  buttress 
walls  or  foundation  beams)  supporting  the 
foundation  extension. 

This  completes  the  design  of  the  two-way  foundation 
extension. 

Example  C.l:  Cantilever  Uall  Barrier 

1:  Design  of  the  foundation  extension  for  a cantilever 
wall  barrier. 

Determine  the  length,  Sickness  and  amount  of  reinforcing 
steel  required  for  the  foundation  extension  of  a canti- 
lever wall  barrier. 

Step  1 . Given: 

a.  Configuration  of  the  structure  and  details 
of  the  wall  which  is  designejl  to  the  in- 
cipient failure  condition  (f'^Rure  C.13). 

b.  Quantity  of  explosive;  tliree  charges  each 
weighing  1,900  1bs  (TNT)  and  located  as 
Shown  in  Figure  C.13. 
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Figure  CJ3  Esassle  CJ:  nissnslofts  of  structure,  design 
detaMs  cf  l>ackv'an  am!  charge  loctitions. 


c.  Soil  data  available: 

Description:  Gravel 
Compaction:  Medium 
Blow  count:  40 

Allowable  bearing  pressure:  5 tons/ft^ 

d.  Design  strength  for  building  materials: 

1 . Concrete  = 4,000  psi 

2.  Steel  fy  = 60,000  psi 

Step  2.  Estimate  the  dimensions  of  the  foundation 
extension  to  be  used  in  the  overturning 
analysis. 

Thickness  of  wall:  (TW)  ==  6.25  ft 

Height  of  wall:  (HW)  = 16.25  ft 

Estimated  thickness  of  foundation  extension: 

TS  = 1.25(TW)  = 1.25(6.25)  = 7.81  ft 

Use  7.83  ft  = 7 feet  10  inches 

Estimated  length  of  foundation  extension: 

Lp  = .45(HW)  = .45(16.25)  = 7.31  ft 

Use  7.5  ft  = 7 feet  6 inches 

For  an  efficient  design,  the  foundation  should 
be  symmetrical  about  the  centerline  of  the  wall. 

Step  3.  Determine  the  soil  bearing  pressure  for  the 
weight  of  the  structure. 

Estimate  weight  of  structure: 

W " 06.25(6.25)  [2(7.5)  ^ 6. 251(7. 83)^52) (150) 

2,000 

= 1 ,045  tons 
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Surface  area  of  foundation: 

A - 52[2(7.5)  + 6.25]  = 1,105  ft^ 

Allowable  bearing  pressure  = 5 T/ft^ 

Bearing  pressure  = 1 ,045/1  ,105 

=0.95  T/ft^  < 5 T/ft^ 

The  foundation  is  adequate  for  the  dead  and 
live  load  condition. 

Step  4.  Determine  the  average  impulse  loads  on  the 
foundation  slab. 

Design  charge  weight  = 1.20(W) 

= 1.20(1,900)  = 2,280  lbs 

Using  the  procedure  and  data  provided  in 
Appendix  A of  the  report,  the  following 
average  impulse  loads  on  the  foundation 
slab  are  computed: 

Wi : Tfj  = 4,400  psi-ms 

^2-  “ A, 800  psi-ms 

W3:  Tfj  = 4,400  psi-ms 

Since  the  structure  is  a cantilever  barrier 

with  a configuration  similar  to  tha  one 
shown  in  Figure  12,  the  reiiwining  loading 
data,  consisting  of  the  times  of  arrival 
and  load  durations,  will  be  computed  by 
the  program. 

Step  5.  Establish  the  range  of  critical  soil  properties 
to  be  utilized  in  the  analyses. 

The  field  description  of  the  soil  [Item  (c) 

Step  1]  and  the  results  of  the  penetration 
tests  indicate  that  the  soil  is  a medium  com- 
pact gravel.  Therefore,  the  structure  is 
analyzed  for  the  properties  of  the  loose  and 
very  compact  gravels  provided  in  Table  1 of 
Section  4.3.  The  properties  utilized  in  the 
analyses  are  presented  on  the  following  page. 
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Modulus  of 
Elasticity  (psi) 

Poisson's  Ratio 

Friction  Factor 


Loose 

Very  Compact 

3,000 

20,000 

0.2 

0.15 

0.6 

0.70 

Step  6.  Prepare  the  input  data  decks  for  the  computer 
program. 

Since  the  structure  is  a cantilever  barrier, 
the  "Normal  Option"  mode  of  the  computer 
program  is  utilized  to  analyze  the  structure 
(see  Section  5).  The  input  data  deck  for 
the  analysis  of  the  structure  on  the  compact 
gravel  is  presented  in  Appendix  D. 

Step  7.  Run  the  analyses  on  the  CDC  6600  computer 
using  the  overturning  analysis  program. 

Step  8.  Inspect  the  results  of  the  analyses. 

A portion  of  the  printed  output  for  the 
analysis  of  the  structure  on  the  compact 
gravel  is  presented  in  Appendix  0. 

The  following  is  a summary  of  the  peak  res- 
ponse parameters  for  the  structure  on  both 
the  loose  and  very  compact  gravel. 

loose  nravel: 

Maximum  rotation  of  structure  = 36.20^ 

Maximum  horizontal  displace- 
me?it  of  foundation  =>  10.00  inches 

Ratio  of  maximum  rotation  to 
overturning  angle  =0.70 

Very  compact  gravel: 

Maximum  rotation  of  structure  = 20.10° 

Maximum  horizontal  displace- 
ment of  foundation  = 3.80  inches 
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Ratio  of  maximum  rotation 
to  overturning  angle  = 0.39 

Inspection  of  the  above  tabulation  of  the 
results  indicates  that,  in  both  analyses, 
the  structure: 

a.  Has  reached  its  peak  response. 

b.  Did  not  overturn. 

c.  Did  not  experience  excessive  horizontal 
(sliding)  displacements. 

Therefore,  the  design  proceeds  to  the  next 
step. 

Step  9.  Inspection  of  the  results  (peak  rotations), 
sumnarized  in  Step  8,  indicates  that  no 
further  modifications  of  the  foundation 
dimensions  for  overturning  and/or  sliding 
are  required;  therefore,  the  design  can 
proceed  to  the  next  step. 

Step  10a.  Determine  the  location  of  the  critical  sec- 
tion for  shear. 

1^  = 6.50  ft  = 78.0  inches 

TS  = 7.81  ft  = 94.0  inches 

Assume  4 inches  for  the  bottom  cover  and 
tension  reinforcement. 

d = 94.0  - 4.0  = 90.0  inches 

(n/d  = 78.0/90.0  = .87  < 5.0 

According  to  the  ptx)visions  of  Section  C.3.2, 
the  foundation  is  considered  a thick  section 
and  the  critical  section  for  shear  is  0.15  1^ 
from  the  face  of  tlie  support. 

1(.^  = 0.15  In  = 0.15(78)  = 11.7  inches  from 

the  haunch 
(see  Figure  C.14) 


140 


step  11a.  Determine  the  peak  shear  (and  the  corresponding 
bending  moment)  at  the  critical  section  for 
shear  and  the  peak  bending  moment  at  the  face 
of  the  support.  These  quantities  are  computed 
for  the  response  of  the  structure  on  the  com- 
pact gravel . 

Although  these  quantities  are  computed  by  the 
program,  the  hand  calculation  is  presented  to 
illustrate  the  procedure. 

The  soil  bearing  pressures  at  three  time  sta- 
tions are  investigated.  Figure  C.14  shows  the 
location  of  the  critical  sections  for  shear 
and  bending  on  the  foundation  extension.  The 
locations  of  the  soil  element  attachment  points 
are  also  shown  in  the  figure.  A portion  of 
the  printed  output  of  the  bearing  pressure 
time  history  is  provided  in  Appendix  D.  The 
pressure  distributions  for  which  the  shears 
and  bending  moments  are  computed  are  shown  in 
Figure  C.15.  The  shears  and  bending  moment 
are  computed  for  a 1-inch  wide  segment  of  the 
foundation  extension. 

The  computations  are  presented  below; 

at  t = 0.04266  second: 


P_  = (479.5  - 233.6)(46.9)  + 233.6 

" — mrStsr — 

= 335.5  psi 

P.  = (479.5  - 233.6)(35.2)  + 233.6 
= 310.1  psi 

Vy  = (335.5  ^ 479.5)(66.3) 

= 27,017  Ib/in 


«u  = C310.1|(78)^  + (479.5  - 3l0.1)(7H 


7(j  ^ 2 


= 1,286,867  in-lb/in 
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SOIL  BEARING  PRESSURE  DISTRIBUTION  AT  t=0.0589t  SEC. 


Figure  C.15  Example  C.1:  Foundation  extension  design  loadings. 
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at  t = 0.05356  second; 


= (549.7  - 3n.4)fl8.6)  + 311.4 

rsTfjor 

= 363.6  ps1 

P.  = (549.7  - 311.4U6.9)  + 311 .4 
= 330.6  psi 

Vu  = (363.6  -I-  549.71(66. 3) 

2 

= 30,276  Ib/in 

= (330.6)(78)^  + (549.7  - 330.61(78)^ 

2 3 

* 1,450.020  in-lb/in 

at  t = 0.05891  second: 

Per  = (18.6/28.3)(396.9) 

= 260.9  psi 

Ps  = (6.9/28.3)(396.9) 

= 96.8  psi 

Vu  = [(572.7  > 396.9}(2)(28.3)1 
2 

[(396.9  ^ 260.9)(9.7)3 
2 

= 30.630  Ib/in 

” 96.8(21.4)2  4-  (396.9  - 96.8)(21.4)2 
2 3 

+ (396.9)[(2)(28.3)][(2)(28.3)/2  + 21.4] 

+ (572.7  - 396.9}(28.3)C(4)(23,3)/3 

+ 21.4] 
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= 1.478,660  in-lb/in 

The  peak  shear  occurs  at  t = 0.05891  second. 

The  corresponding  bending  moment  at  the 
critical  section  for  shear  is  computed  as 
follows: 

M_  = 260.9(9.7)^  + (396.9  - 260.9)(9.7)2 
cr  2 3 

+ (396.9)[(2)(2S.3)][2(23.3)/2  + 9.7] 
+ (572.7  - 396.9)(28.3)[4{28.3)/3 
+ 9.7] 

= 1,106,179  in-lb/in 

The  peak  shear  and  corresponding  bending 
moment  at  the  critical  section  for  shear  are: 

Vu  = 30,630  Ib/in 

Mcr  = 1.106,179  in-lb/in 

The  peak  bending  moment  at  the  face  of  the 
support  is: 

My  = 1.478,660  in-lb/in 

Step  12a.  Determine  the  allowable  shear  stress  for  the 
concrete. 

Since  ly/d  < 5,  the  allowable  shear  is  computed 
using  Equation  (C.l). 

Vy  “ 30,630  Ib/in 

= 1,106,179  in-lb/in 

d = 90  inches 

fc  = 4,000  psi 

= assume  minimum  from  Table  C.l 
= 0.0025 
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v^,  = 0.85[3.5  - 2.5(1, 106, 179;/(30,630)(90)] 


X [1.9^'4OT 

+ (2,500) (0.0025)(30,630)(90)/l ,106,179] 

= 0.85(2.5)(135.7)  = 288.5  psi 

Step  13a.  Determine  the  thickness  of  concrete  required  to 
carry  the  shear. 


Vy  = 30,630  Ib/in 


Vj.  = 288.5  psi 

d = 30,630/288.5  = 106.2  inches 

Step  14a.  Determine  the  amount  of  flexural  reinforcement 
required. 

Main  reinforcement: 

Assume  “a"  block  = 6 inches 

As  =_3L_ 

fsfd  - a/2) 

My  = 1,478.660  in-lb/in 
b = 12  inches 
d = 106.2  inches 
a = assumed  value  of  6 inches 
f$  = 60,000  psi 


(M.6od)(106.2  - 6/2) 
= 2.87  in^/ft 
a = (Asf5)/0.85bf‘ 


f^  » 4.000  psi 

a = (2.87)(60.000)  = 4.22  inches 

r^8t)tl2)(4,0&)) 
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Ac  = (1  ,478,660)02) 

C60,000K106.2  - 4.22/2) 

= 2.84  in2/ft 

Amin  = 0-0025bd 

= 0.0025(12)006.2) 

= 3.19  in^/ft 

Use  minimum  steel:  = 3.19  in^/ft 

Area  of  #11  bar  = 1.56  in^ 

Use  two  #11  bars,  top  and  bottom,  at  12-inch 
spacing.  These  bars  should  extend  from  one 
end  of  the  foundation  to  the  other  end  (see 
Figure  C.13). 

Reinforcement  in  other  direction: 

A3  = O.OOlbTj. 

T(.  = 110  inches 

A3  = 0.001(12)010)  = 1.32  in^/ft 

Use  one  #10  top  and  bottom  at  12-inch 
spacing.  These  bars  should  be  placed  in 
both  extensions  of  tiie  foundation  Oee 
Figure  C.13). 

Step  ISa.  Determine  the  actual  thickness  of  the 
foundation. 

d = 106.2  inches 

djj  = 1.41  inches 

c = 3.0  inches 

Tc  = 106.2  1.41/2  ^ 3.0  = 109.9  inches 

Use  Tj.  = 110  inches 
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C.6 


Example  C.2:  Single  Cell  Barrier  with  Buttress  Walls 

Example  C.2:  Design  the  foundation  extension  of  a single  cell 
barrier  with  buttress  walls.  This  problem  is 
based  on  the  actual  design  of  the  structure  shown 
in  Figures  C.16  through  C.18. 

Required:  Determine  the  length,  thickness  and  amount  of  rein- 
forcing steel  required  for  the  foundation  extension 
of  a single  cell  barrier. 

Step  1 . Given: 

a.  Configuration  of  the  structure  and 
details  of  the  backwall  element  which 
is  designed  to  the  incipient  failure 
condition  (see  Figures  C.16  through  C.IR). 

b.  Quantity  of  explosives:  six  charges 
located  within  the  cell  as  shown  in 
Figures  C.16  and  C.17.  The  weights  of 
the  charges  are: 

W-j  and  W2  = 54  lbs  each 

W3  and  = 146  lbs  each 

W5  and  Wg  = 146  lbs  each 

All  charges  are  TNT. 

c.  Soil  data  available: 

The  soil  at  the  construction  site  is  a 
medium  sand.  In  this  problem,  extensive 
test  data  are  available  to  more  accurately 
determine  tiie  properties  of  the  soil  for 
use  in  the  analysis.  The  test  data  are 
omitted  here  for  brevity  and  the  proper- 
ties are  supplied  in  the  subsequent  com- 
putations as  needed. 

d.  Design  strengtii  for  building  materials: 

1-  Concrete  f^  = 4,000  psi 

2.  Steel  fy  = 60,000  psi 
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step  2.  Estimate  the  dimensions  of  the  foundation 
extension  to  be  used  in  the  analysis.  In 
addition,  determine  for  the  foundation,  the 
area  of  reinforcing  steel  required  to  develop 
the  strength  of  the  backwal''. 

In  this  problem,  the  guidelines  of  Section 
C.2.2  were  not  utilized.  The  thickness  of 
the  foundation  chosen  is  24  inches. 

The  computation  of  the  area  of  reinforcement 
required  for  the  foundation  to  develop  the 
backwal 1 is  given  below: 

Assume  #9  bars  (bar  diameter  = 1.13  inches, 
bar  cross-sectional  area  - 1.00  in2)  required 
for  foundation  to  develop  strength  of  back- 
wall  and  sidewalls. 


= 15.61  inches 


Mpu  = Mwu/2 

Oross-sectional  area  of  #10  bar  = 1.27  in^ 


As  d “ Ag  d /2 

F BW  ^BW 

A.  = (1.27)(24.6)/(2)(15.61) 

= 1.0  in'^  = Area  of  #9  bar 


Use  #9  bars  at  9-inch  spacing. 

Step  3.  Determine  the  soil  bearing  pressure  for  the 
weight  of  the  structure. 

The  weight  of  the  structure  is  2.96  x 10^  lbs. 
The  weight  computations  are  omitted  for  brevity. 

The  allowable  bearing  pressure  for  the  soil 
is  1.25  tons/ft^. 

The  area  of  the  foundation  in  contact  with 
the  soil  is: 

A = (44.5)(38.5)  = 1 ,713.3  ft^ 


Bearing  pressure  = 2,960,000 


2,0000  ,713.3) 

= 0.85  T/ft^  < 1.25  T/ft^ 

The  foundation  Is  adequate  for  the  dead  and 
live  load  condition. 

Step  4.  Determine  the  average  Impulse  loads  on  the 
floor  slab  within  the  cell. 

Design  charge  weights: 

and  = 1.20(54)  = 65  lbs 

W3  thru  Wg  = 1 .20(146)  = 175  lbs 

The  Impulse  loads  on  the  floor  slab  are  com- 
puted using  the  computer  program  of  Reference  2. 
The  computed  impulse  loads  are  listed  below: 


W 


1 


'2* 


W: 


w 


4' 


W 


6- 


Tg  = 449  psi-ms 
Tg  = 408  psi-ms 
= 916  psi-ms 
Tg  = 766  psi-ms 
Tg  = 873  psi-ms 
If,  = 744  psi-ms 


Step  5. 


Since  the  configuration  of  the  structure  Is 
unusual,  the  arrival  times  and  load  durations 
must  be  computed  by  hand  using  the  method  and 
the  data  presented  in  Appendix  8.  The  computa- 
tions are  omitted  for  brevity.  The  quantities 
computed  are  Included  In  the  Input  data  deck 
for  this  problem,  which  Is  presented  In 
Appendix  D. 

Establish  the  soil  properties  to  be  used  in 
the  analysis.  As  discussed  previously,  a 
large  quantity  of  test  data  Is  available  to 
determine  the  properties  of  the  soil.  The 
properties  derived  from  the  data  are  listed 
on  the  following  page. 
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Description:  medium  sand 

Modulus  of  Elasticity:  5,600  psi 

Poisson's  Ratio:  0.33 

Friction  Factor:  0.70 

The  structure  is  analyzed  for  these  soil 
properties  only. 

Step  6.  Prepare  the  input  data  deck  for  the  computer 
program.  Since  the  configuration  of  the 
structure  is  unusual,  the  "General  Structure" 
and  "Special  Loading  Options"  must  be  uti- 
lized to  analyze  the  structure  (see  Section 
5).  The  input  data  deck  is  presented  in 
Appendix  D. 

Step  7.  Run  the  analysis  on  the  CDC  6600  computer 
using  the  overturning  analysis  program. 

Step  8,  Inspect  the  results  of  the  analysis.  A por- 
tion of  the  printed  output  of  the  analysis 
is  presented  in  Appendix  D. 

A summary  of  the  peak  response  parameters  for 
the  structure  is  presented  below: 

Maximum  rotation  of  structure  = 1.58° 

Maximum  horizontal  displacement 
of  foundation  =6.34  inches 

Ratio  of  maximum  rotation  to 
overturning  angle  =0.04 

Inspection  of  the  above  tabulation  of  the 
results  indicates  that  the  structure: 

a.  Has  reached  its  peak  response 

b.  Did  not  overturn 

c.  Did  not  experience  excessive  horizontal 
displacements. 

Therefore,  the  design  proceeds  to  tie  next 
step. 
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step  9.  The  peak  rotation  of  the  structure  is  4 per- 
cent of  the  overturninn  anqie,  whereas  rota- 
tions to  incipient  overturning  (approximately 
40  degrees)  could  be  tolerated.  The  bearinn 
pressure,  computed  for  the  dead  load  condition, 
is  approximately  70  percent  of  the  allowable 
bearing  pressure  for  the  soil.  Based  on  these 
figures,  the  plan  size  of  the  foundation  could 
be  reduced  to  approximately  30  percent.  How- 
ever, this  particular  structure  contains  mas- 
sive steel  vessels  which  are  used  in  the 
manufacture  of  the  explosive.  The  vessels  are 
supported  by  an  equally  massive  steel  frame- 
work. At  the  time  this  structure  was  designed, 
the  size  and  weight  of  the  vessels  were  unknown. 
The  estimated  weight  of  these  items  was  500 
tons;  therefore,  a generous  margin  was  provided 
in  the  plan  size  of  the  foundation. 

Therefore,  no  change  in  the  plan  size  of  the 
foundation  is  made. 

Step  10b.  Determine  the  design  loading  for  the  founda- 
tion extension. 

The  locations  of  the  soil  element  attachment 
points  on  the  foundation  and  the  soil  bearing 
pressure  distributions  at  several  time  sta- 
tions are  shown  in  Figures  C.19  and  C.20. 

The  computations  of  the  applied  loads  on  the 
portion  of  the  foundation  extension  beyond 
the  critical  section  for  shear  (approximately 
20  inches  from  haunch)  are  presented  below. 

The  loads  are  computed  for  a 1-inch  wide 
strip. 


At  t = 0.18586  second: 


FNDN.  EXTENSION 


LOCATION  OF  SOIL  ELEMENTS  ON  FNDN. 


SOIL  BEARING  PRESSURE  DISTRIBUTION:  t=  0.18586  SEC. 


66.1  p$i 


SOIL  BEARING  PRESSURE  DISTRIBUTION  = t*  0.20812  SEC. 


Figure  C.19  Example  C.2:  Locations  of  soil  elements  on 
foundation  and  design  loadings. 
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The  nomenclature  of  Section  5-10  of  Reference 
1 is  used  to  identify  the  various  parameters 
of  the  problem.  The  parameters  are  defined 
as  illustrated  in  Fi^gure  C.21. 

The  area  of  steel  in  the  vertical  direction 
(see  Figure  C.21)  was  initially  determined  in 
Step  2: 
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p 

#9  bar  at  9-inch  spacing:  = 1.0  in 

Pw 

As  = 1.0(12)/9  = 1.33  in2/ft 
b = 12  inches 
d = 20  inches 

Pwv  = 1.33/{20)(12)  = 0.0055  > 0.0025 

Since  > 0.0025,  use  minimum  steel  in  the 
horizontal  direction. 

As  = (0.0018)(20)(12)  = 0.432  in^/ft 

Ajj  = (0.432) (9/1 2)  = 0.324  in^ 

\ 

Begin  the  computations  with  ^*6  bars  at  9- inch 
spacing  in  the  horizontal  direction. 

As  = (0.44)(12/9)  = 0.59  in^/ft 

Step  12b.  Compute  the  ultimate  positive  and  negative 
bending  moment  capacities  in  both  directions 
using  Equation  (C.3). 

Vertical  direction: 

Referring  to  Figure  C.18: 

dy^  = 24  - 3 - 0.75  - 1.128/2 

= 19,7  inches 

dyp  = 24  - 2 - 0.75  - 1.128/2 
= 20.7  inches 

ay  = Asfg/0.85f^b 

° i .33(60,000)  = 1.96  inches 

0785(4, 000) (12) 

Mu  = Asf5(d  - a/2)/b 
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Mvn  " 0.33)(60.000)(19.7  - 1.96/2) 


= 124,500  in-lb/in 
Myp  = n.33)f60.000)f20.7  - 1.96/2) 


= 131,100  in-lb/in 
Horizontal  direction: 

= 24  - 3 - 0.375  = 20.6  inches 
dnp  = 24  - 2 - 0.375  = 21.6  inches 
au  = 0.59(60.000)  = 0.87  inch 

" 0.l5(i^o5Tri2T 

* (Q>59)(6Q.OOO)(20.6  - 0.87/2) 

= 59,500  in-lb/in 

" (0.59)(60.QOO)(21.6  - 0.87/2) 

= 62.400  in-lb/in 

Step  13b.  Utilizing  yield  line  theory,  compute  the 

ultimate  resistance  in  bending  of  the  fouiKla- 
tion  extension.  The  computation  is  perfoYTiied 
for  a trapezoidal  loading  having  the  same 
relative  proportions  as  the  design  load.  The 
computations  are  omitted  here  for  brevity. 

The  location  of  the  yi$ki  line  and  the  ulti- 
mate resistance  are  given  below  and  defined  as 
illustrated  in  Figure  C.22  (Trial#!). 

X “ 119  inches 

Pu  = 49  psi 


Step  14b.  Compare  the  computed  resistance  (py)  with  the 
actual  design  load.  For  Trial  #1.  the  com- 
puted resistance  of  49  psi  is  less  than  the 
applied  lead  of  6G.1  psi.  Therefore,  the  area 
0*'  horizontal  reinforcement  is  increased  until 
the  required  resistance  is  achieved. 
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YIELD  LINE  LOCATIONS  RESISTANCE 


TRIAL  ^1  (STEV*  NO.  13  b) 


YIELD  LINE  LOCATIONS  RESISTANCE 


trial  *2  (STEP  NO.  14  b) 


yield  line  locations  resistance 
TRIAL  *3  (STEP  NO.  ISb) 


Figure  C.22  Esas^ple  C.2:  Results  of  design  computations  - 
Yield  line  locations  and  ultiwate  resistance  of 
foundation. 
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Trial  #2  of  #8  bars  at  9-1nch  spacing  in  the 
horizontal  direction  yields  the  following  re 
suits  (see  Figure  C.22,  Trial  #2): 

X = 135.5  inches 

p^j  = 66.3  psi  > 66.1  psi 

Step  15b.  Based  on  the  location  of  the  yield  lines 
determined  in  Step  14b,  compute  the  shear 
stresses  produced  by  the  design  load  at  the 
critical  section  for  each  sector  of  the 
foundation  extension. 

The  shear  stresses  computed  at  the  critical 
section  for  each  sector  are: 

Sector  I:  Vy  = 180  psi 

Sector  U:  = 166  psi 

Step  16b.  Compute  the  allowable  shear  stress  for  the 
concrete  using  Equation  (C.2). 

Vertical  direction: 

Vc  = 4(1.9/fJ  + 2,500p„)  < 2.28i/f^' 

= O.asEl.S^ToM  + 2,S00{.0055)] 

< 2.23(0.85)47600 

= 114  psi  < 122.6  psi 

Vy  = 180  psi  > 114  psi 

The  thickness  of  concrete  has  to  be  substan- 
tially increased  to  carry  the  applied  shear. 

Try  T^-  = 38  inches: 

Maintain  sa«ae  moment  capacities; 

Vertical  direction: 

Aj  = 1.33  in^/ft 

d_(initial)  = 15.6  inches  for  T^  = 24  inches 


dj. (final)  =29.6  inches  for  T^.  = 38  inches 
1.33(15.6)  = Aj  (29.6) 

Ag  =0.7  in^/ft 

Use  ?*7  bars  at  9-inch  spacing 

Aj)  = 0.6  in2 

Ag  = 0.8  in^/ft 

dy  = 34  inches 

= 0.8/(12)(34)  = 0.0020  < 0.0025 
Horizontal  direction: 

= 0.0025  since  = 0.0020  < 0.0025 
d^i  = 35  inches 

Aj  = (0.0025){12)(35)  1.05  in^/ft 

Use  #8  bars  at  9-inch  spacing 
Aj,  = 0.79  in^ 

A^  = 1.05  in2/ft 
Kecompute  moment  capacities: 

Vertical  direction: 

dy,p^  - 38  - 3-1.0  - 0.44  = 33.6  inches 

dyp  = 34.6  inches 

ay  = 0.8(60,000)  = 1 .13  inches 

0,^5r4,66oKl2) 

“ ^'-gCgO.QOOH33.6  - 1.18/2) 

12 

= 132,000  in-lh/in 
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Mvp  = 0.8(60,000)(34.6  - 1 .18/2) 

12 

= 136,000  in-lb/in 
Horizontal  direction: 

= 38  - 3 - 0,5  = 34,5  inches 
d^p  = 35,5  inches 

3u  = 1 ,05(60,000)  = 1 .54  inches 

0.35(4,000)  (127 

M,,,,  - 1 .05(60,000)(34.5  - 1.54/2) 

= 177,000  in-lb/ in 

M.p  = 1. 05(60, Q00)(35.5  - 1.54/2) 

12 

= 182,332  in-lb/lb 

Comoute  yield  line  location  and  ultimate 
resistance  (see  Figure  C.22,  Trial  #3). 

y = 148  inches 

Py  = 94.6  psi  > 66.1  psi 

The  shear  stresses  computed  for  each  sector 
are: 

Sector  I:  Vy  = 74  psi 

Sector  II : = 98  psi 

The  allowable  shear  stresses  are: 

Vertical : v^.  = 106  psi 

Horizontal:  = 108  psi 

The  applied  shear  stresses  for  botli  sectors 
are  less  than  the  allowable  shear  stresses 
for  the  concrete.  Therefore,  the  design  of 
this  portion  of  the  foundation  extension  is 
cotrpl  ete . 


165 


At  this  point,  check  the  lateral  extensions 
(see  Figure  C.16)  of  the  foundation  extension 
to  insure  that  they  can  withstand  an  applied 
load  of  approximately  two-thirds  of  the  peak 
bearing  pressure  at  the  end  of  the  extension. 

The  lateral  extensions  are  designed  in  the 
same  manner  as  the  simple  type  foundation 
extension  (see  Section  C.4.2).  Since  this 
procedure  is  illustrated  in  Example  C.l,  the 
computations  are  omitted  from  this  problem. 

Step  17b.  Compute  the  load  acting  on  the  buttress  wall. 

The  load  on  the  buttress  wall  consists  of  the 
total  applied  load  on  Sector  II  of  the  interior 
portion  of  the  foundation  extension  (see  Fig- 
ure C.2)  plus  the  reactions  of  the  lateral 
extensions. 

Using  this  load,  determine  the  thickness  of 
concrete  and  the  amount  of  reinforcing  steel 
required  for  the  buttress  wall.  The  computa- 
tions are  omitted  from  this  problem. 


APPENDIX  D 

FORTRAN  LISTING  OF  COMPUTER  PROGRAM 


P.l  General 


This  appendix  contains  the  FORTRAN  Listing  of  the  Over- 
turning Analysis  Computer  Program.  Included  also  are  two  sample 
problems. 

D.2  Computer  Program 

The  computer  program  is  written  in  FORTRAN  IV.  It  can  be 
run  on  the  COC  6600  computer  using  either  the  Extended  (FTN)  or 
the  Regular  FORTRAN  compilers.  A central  memory  field  length  of 
150,000  words  (octal)  is  required  for  compilation  and  execution 
of  the  program  on  the  CDC  6600  computer. 

The  program  consists  of  a main  routine  and  ten  subroutines. 
The  operations  performed  by  each  are  summarized  below: 

Main  Routine  BLASST  - The  main  routine  initiates  the  execution  by 
calling  Subroutine  MULTB . 

Subroutine  MULTB  - This  subroutine  reads  in  the  title  card,  the 
problem  specification  card,  the  soil  properties  and  the  structm'e 
geometry.  Using  tlie  structure  geometry  and  the  soil  properties, 
the  subroutine  computes  the  elastic  constants  of  the  soil  elements. 
Following  this,  the  arrival  times  and  load  durations  are  comouted 
(in  the  "Normal  Option")  by  calling  Subroutine  ADTIM.  If  either 
the  "Special  Loading"  or  "General  Structure  Option"  is  used,  the 
calling  of  this  subroutine  is  bypassed  and  the  arrival  time  and 
minimum  load  duration  are  read  in  instead.  Next,  a summary  of  the 
input  parameters  used  in  the  analysis  is  printed  out.  If  the  res- 
ponse time  of  the  backwall  element  is  required,  the  subroutine 
reads  in  the  design  details  of  the  fcackwall  and  calls  Subi-outine 
WALL  which  performs  the  desired  computation.  The  problem  solution 
proceeds  with  tlie  calling  of  Subroutine  C6TH2. 

Subroutine  ADTIM  - This  subixiutine  reads  in  the  quantities  and 
locations  of  the  explosives  and  computes  the  arrival  times  and 
load  durations  on  every  loaded  surface  for  each  explosive  charge. 
This  subroutine  is  used  in  the  "Noi'mal  Option"  mode  of  the  program 
only.  The  TNT  data  required  for  the  computation  are  contained 
within  the  subroutine.  Interpolation  of  this  data  is  accomplished 
utilizing  Function  AR^IM. 


Function  ARTIM  - This  function  performs  a geometric  interpolation 
of  the  TNT  data  contained  in  Subroutine  ADTIM. 


Subroutine  WALL  - Subroutine  WALL  computes  the  maximum  response 
time  of  the  backwall  element  designed  to  either  the  incipient  fail- 
ure or  the  post  failure  fragment  conditions.  Four  support  condi- 
tion options  are  considered: 

1.  One  side  fixed,  three  sides  free. 

2.  Two  adjacent  sides  fixed,  two  sides  free. 

3.  Three  sides  fixed,  one  side  free. 

4.  Four  sides  fixed. 

Subroutine  CGTH2  - This  subroutine  computes  the  following  quanti- 
ties  in  the  "Normal  Option"  n»de  of  the  program: 

1.  Weight,  mass  and  mass  moment  of  inertia  of 
structure. 

2.  Location  of  center  of  gravity  of  the  structure. 

3.  Areas  of  all  surfaces  (perpendicular  to  the 
plane  of  motion  of  the  structure)  which  are 
directly  exposed  to  the  blast. 

4.  Locations  of  the  centroids  of  the  loaded  areas 
relative  to  the  center  of  gravity  of  the 
structure. 

5.  Horizontal  and  vertical  components  of  unit 
vectors  normal  to  every  loaded  surface. 

In  the  "General  Structure  Option",  this  computation  is 
bypassed  and  the  data  listed  above  are  read  in  on  punched  cards. 

The  next  operation  performed  by  this  subroutine  is  the 
computation  of  the  load  history.  In  the  "General  Structure"  and 
"Special  Loading  Options",  the  surface  loading  data  (average 
impulse,  arrival  time  and  load  duration)  are  read  in  on  punched 
cards. 


After  the  completion  of  the  load-history  computation,  the 
load  history  is  printed  out  and  the  program  proceeds  to  the  res- 
ponse computation  by  calling  Subroutine  GRMT2. 
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Subroutine  GRMT2  - Subroutine  GRMT2  computes  the  response  of  the 
structure  to  the  applied  loads.  This  involves  computing  the  re- 
sisting forces  in  the  soil  and  solving  the  equations  of  motion  of 
the  structure  (see  Section  2.2)  to  determine  the  accelerations, 
velocities  and  displacements  of  the  structure.  During  the  compu- 
tation, several  auxiliary  subroutines  are  called  to  perform  spe- 
cialized calculations.  These  subroutines  are: 


1 . 

Subroutine 

DEDE 

2. 

Subroutine 

UPCHK 

3. 

Subroutine 

R3 

4. 

Subroutine 

FDN 

After  the  numerical  integration  is  completed,  this  subroutine 
prints  out  the  response  of  the  structure  in  the  manner  described 
in  Section  6. 

Subroutine  DEDE  - This  subroutine  monitors  the  displacement  time 
hTstory  in  each  of  the  vertical  soil  elements  to  determine  whether 
the  element  is  in  either  the  loading  or  unloading  condition  (see 
Section  4.2).  When  the  element  is  in  the  loading  condition,  the 
subroutine  computes  the  elastic  resisting  force  in  the  element  and 
the  moment  of  the  resisting  force  about  the  center  of  gravity  of 
the  structure. 

Subroutine  UPCHK  - This  subroutine  computes  the  friction  force 
acting  on  the  foundation  for  the  purpose  of  detemining  whether 
the  foundation  of  the  structure  has  separated  from  the  soil. 

Subroutine  R3  - Subroutine  R3  computes  the  monient  of  the  horizontal 
resisting  forces  in  the  soil  about  the  center  of  gravity  of  the 
structure. 

Subroutine  FOH  - This  subroutine  computes  the  shear  and  correspond- 
ing bending  moment  at  the  critical  section  for  shear  on  the  foun- 
dation extension  of  a cantilever  wall  barrier.  The  bending  moment 
at  the  face  of  the  support  is  also  computed.  The  computation  is 
performed  at  every  integration  time  station. 

The  following  pages  contain  the  FORTRAN  Listing  of  the  Overturning 
Analysis  Program  for  the  CDC  6600  computer. 
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CUi  5U|9CUtl5r  99  TO  r|llO  I«c«fll9 

C'7  MCt)7CI*t9l  50U  ■C5|5t5lftCt 

9534 

< 

rc«o<5  itcut  CC5II*  or  ••5»it«  cr 

9431 

C 

1915 

999599 

cut  •IU5,5*»9.m»U,9|J»««55I 

M35 

999911 

CC  TO  799 

9551 

99951? 

H5 

W«ltJl|V  «u«5t5<lt Jim 

9951 

•99511 

WU*MJ»1  tt  >*5«5t5l4UM|mi 

9557 

99U?5 

1559 

«5t5  tl*! 

• 953 

1995M 

00  745  t<l,55 

• 555 

5999?1 

50  l0155.59ltTM7« 

9955 

9994l« 

59 

11  t5»«ni5W5|«5| 

99H1V 

41 

•«5l7lTl*9«5 

9951 

•19451 

53  IQ  144 

5959 

999551 

51 

9955 

•99559 

to  to  54 

94V9 

999^55 

55 

9«'59*l(l*«MV*5Pttl«55»< 

9551 

999959 

54 

(MUVWQQkt  lltl741»74A«H5 

9957 

591951 

141 

9«M9tl»«o9«5«Mtl 

9453 

<99599 

?5« 

-«Bll  5 9ptt»MI 

M55 

991949 

754 

59«5 

999U? 

5V*U5t(Uvl  ^•95IW^t|IJv<i9 

• 959 

■91519 

t«  tf«KVt71l9»lf  l9«7ill 

9lt4M 

?7ji 

|71tV»77|  1.7711,1117 

• 959 

919519 

11*1 

tl  t5tt71t7,77l5«71|5 

5H5 

991115 

7M7 

55*-5» 

9959 

991415 

10  7M1 

5451 

999411 

7H 1 

If  IWjl711l*?M7,7H5 

9547 

9MMI 

7715 

t«  tvVMl 77 15  *7715, ?7i r 

• 551 
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Mini  ««  nitt  • ••IIMIXMIIII-Vt.l  MM 
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Mit;<  «)t  If  •;?,•;?  M»r 
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Mini  Iff  If  i»iiMf n -i»»i»i»M.«.*».»t» 
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M««»l  ir  to  •!•  «MI 

Mlltl  »t«  Mltirt{»««  •'••f 

MIM«  to  10  Ml  «M1 

ttii«i  •?»  nifmi«M  «M* 

MIMI  Minifl*  ••lIMIUMff'^O:!  »M» 

Minr  *!•  •ltl••••f••ll<«lltl  * oil*  MM 

Mint  n to*iu.'*  Mit 

Miirt  IM»Mnift.i;*i*.it*i«  «»*• 

«n»M  i,-*n  n «»•» 

MlMl  lt«M  •*-»••.»  «‘*n 

Mius  e,-  i»t;»  11*1, M •'ll 

•tiiM  If  n»miiiMJf,it«r*,i;*?»  Mit 

itiiM  i.-Mi  xn*wu***ie*  OMi 

Miiii  iji*i 

Mini  t*  14  |J*%»  «»M 

ifinr  it*M  iMi *if«w*»tM  Mil 

MUM  IM»«  It  l♦l♦l•••M»•^«l,l.••^,^»•ll  **>* 

M«lt«  INfl  *»»•*• 


1-1 


NNvfAea*  iu4t  «}•>  )M  iv<4/» 


lt*i 

ttlt 

IIIKf 

carl 

fiC  UASt 

lt2i 

•tu 

Uic«tl«t 

8t;s 
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•V» 
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ur 
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MU 

IHI 

t«  u»«tnir%t»iMi»iHr 

• Mt 

mi 

its% 

Hurt 

CC  tc  Itj9 

• 9M 

tiura 

tut 

l«ll 

tittr; 

CC  IQ  Itu 

•Mr 

MUM 

trit 

•MUiM  lt•l•ft«l1»C«9att• 

ttu 

Miirr 

irj« 

l«St 

tiutr 

eCaUM'Jt 

•Mt 

MU1« 

t*  utt(9( 

Ital 

««IMI 

IM 

t«  iiiiiMi.u9t,mi 

IMf 

itirti 

mi 

M*IJ 

• ♦*> 

MUt« 
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CC  1C  ur 

ftu 

tuc.t 

• Mt 

Muir 

e««t 

IHI 

tttrrt 

19. t 

•Hr 

U\tU' 

W .‘rji 

• HI 

•tUf? 

n*uit»*c.t 

• HI 

Ml??) 

UtcliKt 

•Ml 

ftlM« 

'j|*«CUtl*9.t 

mi 

tlUT* 

• M? 

r»?i 

•«M 

MMM 

• Ml 

< 

• M« 

c 

•»  tc  t*  «c*txe»t«b  MU  Htntucr 

im 

c 

MKCt  M;v«  «•*!*  ^ «?# 

•Mt 

t 

•Mi 

• Ml 

ts  ts  ;m 

• lit 

ttlM.* 

Mt 

fit  t.i,  a 

mi 

Hi;«t 

C( 

me 

Mif*r 

tut 

time 

5«U 

• Ml 

f«]H« 

tt:  cell 

•III 

ttir«* 
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tut 

H)K« 

CO  tf 

mi 
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mi 

f 

•Ml 

r 
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001267 

2836 

0976 

001273 

FFdJsftJ 

6975 

00127% 

IFU-L5WI  2«l2,28t3*2»12 

6976 

C 

••  0977 

C 

COMPUTE  OSIPLACfHEMTS  OF  STROCTORf 

USING  RECURSION  FORHULA 

0978 

C 

• » 0979 

001276 

28«  > 

KIJP1U1«7IJ*U*?.0  -XSAVf(t»  ♦ 6JitJTn*0ELT»0fLT 

0980 

I0l3«% 

CO  TO  2816 

0981 

00130% 

7812 

XlJPtm  = XlJC!l»2.6  -K!j**un  *AXIJ«XI  •DFlT*OELT 

098? 

C 

••  0983 

C 

COMPUTE  AVERAGE  VELOCITY  OVER  TIME 

increment 

0986 

c 

••  0985 

06131? 

28U 

MYdjaixijpun-tuaiJ/oEu 

0986 

001316 

U-1  ;aVF»?8(,28U»201 

0987 

061320 

2611 

xSAVE(t)3Xi3i;n 

0988 

oin?2 

2<1 

CONTINUE 

0989 

01132$ 

IFU-LTIHEl  302?»3Ol9»3622 

0990 

001327 

3U19 

ANUMP*NUNPT 

0991 

601331 

MUMP T*NUHPT/NCEL2 

0992 

80133^ 

3022 

JNUNsjNUN  ♦ 1 

0993 

061336 

TE(JNUNoMUNFm630»  3826  >3020 

• 996 

0113%; 

3020 

099$ 

6013%! 

KNUHsKNUH  «t 

1996 

66t3%3 

TFUNUH-100  0U6?l>30ri>2A3 

0997 

0113%$ 

36?) 

l IME  IKNLHTaTj 

1998 

II13%7 

00  96  Ktt,3 

8999 

001351 

R^SP4K»XNLNI»X1J<K: 

1000 

06135% 

RS3tSTlX,XNWT*R«KI 

1001 

06U57 

ACCEl  (K>KMUM)iiXT3CK» 

1012 

60136? 

96 

VELOC(KtKMUH|sRl(Kl 

1003 

611367 

00  92  x«l,h5 

1066 

001376 

92 

STRESJiX,KNUMI*SN(X» 

1015 

601376 

3630 

IFUUPU2n  S03>  363>  306 

10B6 

I81%I6 

303 

CANA  « AnCCT  * XlJPlIS) 

1007 

601%62 

306 

TMXIJPlun  36lt  301,  362 

loot 

t01%6% 

302 

IF(XIJPU3)  •XUtSII  305,  501,  501 

1689 

Q4U07 

561 

ALPHA  « ANCCT  f XUPU3) 

1010 

Q61%11 

TLsO 

toil 

001%12 

IF  (ILPMA  > t.571>  502>  503,  $63 

1812 

6ll%15 

$02 

tF(XI3Pt(3i'INAXmQ5»365,S36 

ms 

C 

••  161% 

c 

PRINT  MESSAGE  INQICATIIIC  STRUCTURE 

OVERTURNS 

101$ 

c 

••  1816 

llt%2t 

563 

WRITE  (5,5061  fj 

1017 

I61%26 

-HAYT»xr»PU3l 

1016 

661636 

U»1 

1019 

I61631 

$66 

FORMAT 1////,?9M  STRUCTURE  OVERTURNS 

AT  TINE  >F10.51 

1020 

I61%31 

GO  TO  ?I3 

1021 

I61%S1 

305 

U ■ -1 

U22 

6fl%3? 

tF(XtJPlT3}-AKAXT)3et,36t,336 

I Jt  \ 

661%35 

336 

AMAxT  « X IJPH3) 

1026 

611637 

361 

TFINMAIU  356,356,351 

1029 

661%61 

351 

7FnWAUl  356,356,5$« 

1026 

601%%3 

356 

lF(rj«TWALLt3$6*352,352 

1027 

661666 

352 

IF(VFI3$3t353,S55 

1026 

CRKT? 
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IIU5I 

QU««( 

MINS7 

QILNM 

0IU73 

MU7¥ 

CI15I0 

fidist? 

eoi9a« 

ttl«67 

MtaU 


8I1S13 

UX%tk 

sit^7r 

MlSll 

dq\953 

iitau 

Q41958 

dll^S? 

omN? 

eits^t 

iiiv^a 


311573 
l5Ua  1 
MU?3 

dti*?« 

iiiara 

II3R«8 

111858 

fliun 

iiirtf 

«U?1f 

i%%rr^ 


iiirtF 

•I1M5 

Mini 


M5  MltTC{5«t5ttdlTJ 
tt*-l 

CO  to  784 

183  tdPtT(<5«tM8|Tj 
XMHl«l 
158  JP  • J*l 

If  UP  - HJl  ZnTt  781,  7tJ 
«82  00  71  !«t,l 

irtvtjpitt)  •ONiim  ) 77,}7,?3 
?3  OHitrti  « luputi 
?t  tniUPUtl  -ONlMttU  75»75»?| 
28  OMIMIII  • VUPUa 
21  CONflNUC 

tr  tl>fONI28, 26,2t 


C CALL  8ueP0UtTHC  FOM  10  COttPuTC  SHEARS  AND  BCNDT^O  HOMCMTS  FOR 

C FOUKOATIQN  FlUKSION  OF  CANTKCtfCt  NALL  8AR02EP 

C 

77  CALL  FOMAS, LI, teP,llH«SM»lMAAt»1NJlV<XMA«,Ul 
28  tF(SHINS)  - MNAf>  12,12 ill 
33  PRHAV  s IMKSI 
it  JUM 

00  U7  1«t,.1 
vUPtl1»«iUlt> 

tit  Atijpm«>iuui 

TT0(*ltJ{2l  *U(ASI*S|Nt«tlUII  • tOtlC 

tr  urof  •TiMAmm,tm»iu 

lilt  TtH*8«Tf0f 
ftMltvtJ 
00  10  lit 

281  0€Lf  • OtLt  2R0CU 
MUHPTaANUAP 


C PRlHt  S'RUCtUIAL  8H0  tCSPCNSC  PA»AM(1{»$  rci  RACK  NFLl  UtKfMI 

5 

trihVALU  288,288,788 
288  tfAMCt^itmOA'lC 

00  TQ(28F,288,288«288I»I. 

212  MlTE1«,l8|llQC,PV.rm,RU,VU,VV8U. 

CQ  IQ  2881 

288  32111285, 288,211 

781  N«Srn8,t197IO€,PV,PM,T,ros,tALllu,9UNP,tU,lL^ 

CO  to  2888 

711  N8|121i*lUll0C,PV^Pii,l«f01,«<Lt)U»SUM,»^,tUR 
7888  N«nri1.1188lt|,llR,IW8U 
7181  irt>riM8,?A8,28S1 
7811  NtUCI1,t3l4l8F 
?A8  tFUetlllt, 857,114 

t 

c P8Tm  oimiciNfttts  of  imlciuri  imq  icu  ttsisiAMtr  fcpc;i 

418  vtl12ff ,«1l4CfU 
NRfTf  (1,777) 

¥Riir  «4, 'll  I 


1025 

tlJ8 

1831 

1832 
1893 
1018 
1031 
1038 

1037 

1038 
1815 
1880 
1881 
1082 
1181 
1088 
1081 
1888 

1887 
1088 
1885 
1810 
toil 
1812 
1813 
1118 
1811 
1818 

1817 
1811 
1815 
1880 
1181 
1082 
1181 

1818 
1881 

1888 
1887 
1848 
1181 

1871 
IIM 

1872 
1171 
1878 
1171 
1878 

1877 

1878 
1875 
1818 
1881 
1882 
1182 


itmi 

CRHT?  MHV740N  RUH  V7.S  PSt  SAO  0O/74/7S 

MMfr(1»tM7| 

17.17.31 

1004 

KltNl^O 

1041 

6I17«7 

DO  20'  J • 1.  KNUH 

1084 

• Q17M 

KLlNl«NLIIt(  * 1 

1087 

•It79« 

tr<4LtKe-44l7004,70Q«,70IS 

108.8 

901740 

7003 

M.TN7>\ 

1081 

•01741 

MmtS.lDKCCIC 

1010 

•01741 

HRm(4(70?l 

toil 

••1773 

UVtTC (1,7041 

1017 

••1777 

UftlTKI, 11071 

101s 

• •7003 

7004 

«fONiRf1P(t,JI>K»1tMtRC1»(S,J)» 

1114 

007Q13 

•F1F(J,JI«17*?14«ReSP(S,JI 

1111 

007I14 

70S 

N»ttE(1»104lTlltC(j|,  (Rllf  (K,4t,K«l,31  » tl»e«tS7  (K,  Jl , K9 1 » SI  »XF  ON 

1014 

••70oe 

TMNVFL  11707,1707,1701 

1017 

C 

••  1018 

c 

HTIKT  ACCFlFftOTQONS  4H0  VElOCIMCS  Of  StURCUTRC 

1011 

c 

••  1100 

• Q7Q41 

1701 

WltfC(1»t1)KOC1C 

1101 

••7147 

M»tfCt1,7«tl 

1107 

•17043 

WRITE <1,7041 

1183 

••7147 

wPtTEtl, 11041 

1104 

•«7«73 

tioi 

••7074 

00  ?••  J«  l,KNU«t 

1104 

••7074 

M.TME»HKt4r  » 1 

1107 

C07t00 

tMN|.tNC-4RI7007, 7007,7044 

1111 

107107 

7011 

MUhEal 

1101 

»071«3 

wtttEditiiKoeie 

1110 

• •7111 

utm(4,70ii 

im 

•Q7114 

l»1TC(1»?44» 

1117 

•17171 

MtlTEd, 11041 

1113 

•17174 

7107 

VTC!N«vH0C(l,Jl*u*VUQC(S,J)«tC01(tCSF(S,Jl/17.714ll 

1114 

•17141 

?o« 

W»tTE(1,miTtMCUI,UCCU<t,JI,««l,Sl,  lveLO€(K,JI,I*l,  SUWFDM 

ms 

•17174 

701 

rOWMft  |1I,77N  OCCUFROftOH  07  C»S.  »1SX» 

1114 

•17174 

1707 

uiH  vucem  cr  c.c«  /i 
trden  4i7,i?os,  417 

1117 

1118 

C 

*•  1U1 

e 

PWlliT  ion  9MR1NC  PUCSSUOfS  St  sou  EUHFKT  STtSCHHCHT  POINTS 

1170 

e 

••  U?l 

••7171 

17«) 

UN|Tr(1,i1IKOMC 

1177 

•4770  7 

00  7704  J«1»N1 

1173 

• 12??Q1 

77«4 

TTUI»J 

1174 

•4771 • 

UBMC  t1,70lMtV(J»,>|,Ntl 

1171 

•4»777 

MnC(1,77QM 

1174 

••7177 

Mim«o 

ll?r 

•0777t 

00  704  J»t,  «NU« 

1174 

•09717 

HimctMi rwe  *i 

1171 

• t)7M4 

ir  I44(MC*41I7»H,7004,7091 

mo 

•O17J0 

7«04 

H t4E«l 

mi 

• OPIIF 

W«ITT(1,t«IKOCSC 

1137 

•C(74l 

•RTTf  (1,7011  mill  »K«l, 411 

ms 

tiimo 

Mm(4, 77011 

IU4 

•07744 

70«4 

V»1t(  (1,7071 1l»tui,  title  1SU,JI»««1, MSI 

mi 

• II73I7 

r«4 

CONtlNV^ 

IIU 

• II7M4 

• 47 

If  IvtPlii  tl••,t|•l•u•s 

iisr 

407H7 

ItOO 

Tf«40«-W«4B 

IU1 

'fm. 
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tieJil 

l«lTC(4|Ut3>rrfl4X»TTM47 

1139 

••S32« 

1U3 

70RM4T(//t7|99MN4l  VCfIT  ICAL  0ISPL4C£IICIiT  QP  TOC  ABOVE  TNC  MOUND  t 

1140 

|Ft.3,t|H  AT  TtMC  •fl4.«/» 

1191 

M?3}| 

GO  TO  1149 

1143 

H<32t 

till 

MtTEI5.tl44)TT*UXtTTN4l 

1193 

liZSll 

1144 

70KH4T(//tt,«4HHt9XNUH  VERTICAL  OIlRLACCKeHT  07  TOC  BELOH  THC  CAOU 

1144 

IHO  •fA.ItlQK  AT  TINE  »Ptl.S/l 

1149 

M2331 

1149 

Ifl  tOTCI  49l*444»444 

1144 

113333 

414 

CUL  MULTB 

1147 

mss^ 

414 

IFIlLt  19t  Ut  317 

1144 

M23SI 

347 

If  IICll3lf,3M*44 

1149 

H33%l 

19 

tfltCtl7Ql,7gi«91|fi 

1191 

II3393 

7U 

1M  ton  7e9,7«9t7«3 

1141 

ltS3«« 

7«3 

OTM  ■-TVOT  *TU 

1193 

Il33«ft 

QTVTU  ■ QTM  / tMOT 

1143 

««33%7 

irtOfH«3.0)407»917,707 

1199 

IC3SS3 

717 

If  ( 0THT»-4.aMI4«7t  447,  744 

im 

••3399 

70I 

0TN3  ■ OTm  /3* 

im 

••?397 

TW  • TVCT  • QTW3 

1147 

••334t 

T5  > C9  •TM 

1144 

••2393 

CALL  C4TM3 

1149 

••23«« 

709 

lOT  • -1 

lUI 

• •?349 

TMOT  • TV 

1141 

• •3347 

TV  • Ty/3. 

1113 

••3373 

T3  ■ fV*C9 

1143 

••3373 

If  i TV«t?.  1 7|4,  704,  784 

1144 

• •337« 

7«4 

Wttfe  1 4»71B1 

1144 

• 03MI 

711 
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D.3  Sample  Problems 


Two  sample  problems  are  presented.  The  first  problem  is 
the  cantilever  wall  barrier  that  was  analyzed  in  Appendix  C.  The 
configuration  of  the  structure  is  shown  in  Figure  C.13,  The  input 
data  are  given  in  Figures  D.l  through  D.4.  Immediately  following 
the  figures  is  a portion  of  the  printed  output  of  the  analysis. 

The  second  problem  is  the  single  cell  barrier  that  was 
analyzed  in  Appendix  C.  The  configuration  of  the  structure  is 
shown  in  Figures  C.16  and  C.17.  The  input  is  given  in  Figures 
D.5  through  D.22.  A portion  of  the  printed  output  of  the  analysis 
is  presented  on  pages  236  through  249. 
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Figure  0.2  Input  data  Kheet  - Example  0.1:  Card  Types  3 and  4 


CARO  TYPE  5 -Charge  No. 
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Figure  D.:?  Input  data  sheet  - Txample  n.I:  Card  Type  F 


Figure  D.4  Input  data  sheet  - Example  D.l:  Card  Types  5 and  7. 


eXAHPte  PROBLCH  O.t  CAKTILCVf*  H4U  BAKRUK  ON  COMPACT  CRAVEl 


AUMfiER  OF  CHARGES  3 
AUMER  OF  HALLS  1 
AUM8CR  OF  SOIL  SPRINGS  lA 
AUMOCR  OF  TINE  inCReACHTS  SIM 
LEAST  TIME  Of  ARRIVAL  OF  BLAST  TL  NEARCSt  POlHTISCCt  *miA 
LEAST  DURATION  TXHElSECI  .IISSS 
PRINT  FREQUENCY  21 
IMTfCRATlON  TINE  INTERVAL  itStCI  2.ANSAE-IA 
tNT£ CRAY  ION  T|AC  INTERVAL  2(S€CI  E«ATAR4E*I4 


SOU  PROPERTIES 


NOOULUS  OF  ELASTICXTT  EUPSt)  EIIM.M 

POQULUS  OF  ELASTtCITT  EKPStI  tllM.II 

AOiVAl  stress  at  HHtCK  NQOULWS  CHAUCCSTPST)  2MQ.M 

COEFFTCIGNT  Of  FRICTtON  .71 

PCtSSOAS  RATIO  .tS 

SKAFE  FACTOR  BX  l.M 

SHAPE  FACTOR  IT  2,11 

NHl/SPRtMCiLA/lHI  TSTTRA.M 

KVt/SPRfKQUO/tKt  AlilSI.SI 

AH2/SPIUCTLQ/XHI  MfYSA.M 

AV?/SPRlNQUa/IN|  lUMA.ll 

MHt  TOTAL  Ua/TMf  m7>vA»lS 

IVl  TQ1ALTL8MMI  miaiS.M 

HH2  TOTIULO/fAl  fSTTSAl.lS 

cv?  ToTALua/iNi  itAiaas,!} 


STRWCTLUE  WONCT3V 


6ACRVAU  TNCIKACSSTINI 

■ATTC  OF  SLAl  THICKNESS  TO  BACRHALL  TRTCKSS 
lEMCTH  OF  iACKHALLIFTI 
HTTAHT  Of  lACKVALLTFIt 

UTIO  OF  HIQTH  Of  FLCQR  TO  TOTAL  HIOtM  Of  BASE 

•ATIQ  Cf  HflCMT  Of  HALL  TO  BASE  Of  STRUCTURE 

HtOTH  Of  RISE C tut 

NTOTN  Of  LQAOCO  ARfR  ON  BAVUMI 

SLAB  TMlCKHESSttHI 

RATIO  Qf  hall  TMICtmSS  TO  HlOTH  OF  STtVICTURE 
MiCtITQHTAL  OISTAHCE  fRQN  CC  TO  REAR  FACE  OF  HALLttMt 
VCtItCAL  QIS1ANCE  FCCN  TCP  OF  fiCCR  SLAB  TO  CCttNI 
QTSTtMCC  ftQN  CS  TO  REAR  CllO  Of  SLAB<|M 
FASSTLB<SCC**2fVNI 

MSS  NONEMT  of  tH|ttUUB'-l«l*SEe««?l 
KIMiTKBtl 


TR.III 

l,2S3 

R2.III 

IA.2SI 

,4fl 

,TAS 

BS4.RVIA3 

BS.SM 

RA.MI 

.2SA 

•sr.fti 

T,fU 

•IBY.ASR 

VAtS.B 

B,N7||TE*TT 
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CXARHC  MOBtrM  0*1  CAMTUCVCR  MALL  BARMXeft  ON  COMPACT  CBAVEL 


LOAOXNC 


charge  NUMEEA  1 


HIH  OtSTAKCE  FROM  CHARCC  TO  NALLtFTI  G4A79 
CHARGE  MElGHTaOSl  EElA.t 
tRPUlSe  ON  MALLtPSt^RSl  3A||«I 
IRPULSe  ON  ELCQR4PSI-MSt  AAOfi.D 
HIM  OtSTAMCC  FROM  CHARGE  TO  EDGE  OF  MALLCFTl  i4.SH 
HEIGHT  OF  CHARGE  ABOVE  FLOOR  SLAOtFTI  9*100 

AREA  OF  FAC£tlN**TI  lEUIO.O 
HORTTOHTAL  DISTANCE  FROM  CENTROID  CF  AREA  TO  CSHN)  -S7.S 
VERTICAL  DISTANCE  FROM  CENTRQXO  OF  AREA  TO  CGtlNT  -RS*7 


PRESSURE  TINE  CURVE 


POINT  1 PRESSURECPTtl  RSS.t 

TINE  tSEC) 

.4419? 

POINT  7 PRESSURE  (PStI  *• 

TINE  TUC) 

•447S1 

TINE 

FORCE  1 

FORCE  S 

HONENT  ABOUT  7 

TSECt 

ILBSI 

ILIST 

(IH-LISt 

7 

.Aom 

t.lR7Qf*M 

• 0 

S.7S4SE4IS 

3 

*tll4f 

t.l4l4C«H 

.9 

S.41I7C^|S 

4 

*HI«7 

1*III9E*H 

.1 

S.944SEMS 

S 

.Hits 

S*?tStE«l7 

*• 

I.7779E4IS 

« 

• tlUl 

S*nS4EAl7 

.1 

4»97eU*IS 

7 

.11147 

I.SSEIEtir 

• I 

•*I99SE»4S 

• 

«H1S4 

S.S4«Sf*R7 

• V 

7.44S7E«4S 

1 

*11771 

•*1RSIE*I7 

• 1 

7.34S9C«4S 

11 

*H74r 

7*lim*R7 

• 1 

7*m9E«IS 

It 

.11774 

r,4tm»Rr 

• 1 

B*4S7lE»IS 

17 

.11911 

7.SJ44£*R7 

.1 

4.S174E»IS 

U 

.Q9977 

«,49ISf*R7 

*4 

S.14I7LMS 

14 

.HSS4 

4«T?4m*tr 

.4 

S.474UAIS 

n 

S,BRS4E*07 

• 1 

S«7A43E*44 

t« 

.M4VI 

S.SSI9EAR7 

«• 

4.S3I7EHS 

17 

.11414 

9.tl47E«V7 

4.SSriE«4S 

U 

.11441 

4,  Mm*RF 

• 4 

4.F474E*|S 

n 

.M44R 

4.|4S7E«R7 

• 1 

9.S|34E»4S 

71 

.lists 

1.SS47C*I7 

• E 

9.SSS4E*|S 

71 

.11441 

9.fR7SE*|7 

.1 

3.7TS7C*IS 

77 

.IIS4I 

l.lSSlfAf? 

.1 

7.B71IEMS 

fS 

.IISSS 

7«MSW«I7 

.t 

?.S744f»|4 

7A 

.11477 

7*4sm»ir 

*1 

7.1A74EHS 
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.11441 

7*R4RSC«i7 

.4 

1.43I4EHS 

74 

••I47S 

I.4SSIEAI7 

*9 

1.4447f«4S 

77 

.11117 

l.ttlUAtF 

*4 

l*144irHS 

7R 

.11774 

I.S7«9r*«4 
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A.4S43M4I 

IS 

.HTSS 

S.t44Vf4|4 

.4 

4.4I44EMI 
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1.174SEH4 

AREA  CF  FACflTN«H»  S»ISR.« 

MCtirCMiAi  QIStANCf  FROM  CENTROID  CF  AREA  10  CGIINI  oH.S 

VCRtUU  OISTRAtCf  FM^  CENTtOlO  MF  ARIA  TO  CGlUl  T.l 
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tt.«  CUtVC 


Pdlkl  1 PICSSMEIPSII  123».< 

TIHnSECl 

POtKT  2 pacssuacipsii  >o 

TIHE(SEC» 

TtHC 

20«C€  X 

80RCE  T 

tseci 

IL8SI 

UBS) 

1 

.eoitt 

.8 

6.9988E»Q7 

2 

.10133 

.8 

8.8899E»87 

3 

• 8IIM 

.8 

8.9237E»07 

% 

• 8IM2 

.1 

8.1828E»I7 

9 

• MU3 

• 0 

9.9819£»87 

e 

.Itl98 

• 8 

9«89e3E*07 

7 

«88U7 

• 8 

9.3792C*Q7 

1 

•tim 

• 8 

9.it8tC»87 

9 

*88228 

.8 

9,8989e»87 

18 

•882%2 

.8 

9.99$n«l2 

it 

. 18279 

• 8 

9,S39IC»8r 

12 

.11311 

• 6 

9.8739r»B7 

13 

.81322 

• 8 

3.  81231*87 

U 

.61399 

• 0 

3.99126*87 

19 

.883t«. 

• 8 

3.29816*87 

u 

.91988 

• 8 

3.82896*87 

17 

*••939 

.8 

2.78776*07 

It 

.88981 

• 0 

2.98886*87 

11 

.88988 

• 8 

2,29996*87 

28 

•88811 

• 8 

1.98936*87 

21 

.88991 

.8 

1.72316*87 

22 

•88981 

• 8 

1.98216*87 

23 

.88999 

•e 

1,28892*87 

29 

.88822 

• 8 

9.39726*18 

29 

.IM98 

«8 

8, 78976*88 

29 

.88879 

«8 

9.17936*88 

?7 

.88282 

• 8 

i.m9<*8« 
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• oftMa 

MOMENT  «eOUT  2 

(iN-ies) 

-5.7J0^£*0^ 

•9.aa^iE*a^ 

•4.a$  32t*0<) 
•^•^3782*09 
-^.2223t*09 
•4.iaai€*a9 
-3.7119E*09 
•3.97Ut*l9 
•3.36iaE»09 
-3a«9ic*i9 
-2.97irr»09 
•.2»7ia2r»B9 

•?,2l3aC*B9 

-2*aa79C*09 

•i.is?ar«e9 

-l*6378t*09 

*uame»«9 

•9.9»7tt*0« 
•7.7f?IC»ll 
-f.99l2C»ta 
• 3«^43«E«M 


CIAHPLF  F0QBLC»(  0*1  CAMTUtVCt  MALL  BAARRA  ON  COMPACT  CRAVCL 


LOAOINO 


CHK0CC  NUMOCR  2 


NtM  QtStAHCE  fRON  CMARCT  TO  MALLIFT)  A«A79 
CHARG€  HetCHTUeSl  ?2AD*0 
IHPOtSF  ON  HALLIPSI-M^I  SIOO.I 
IRPULSF  on  FLOO0IPSt>N^I  AASfi.O 
MtN  OISTANCC  fRON  CHAMGC  TO  EOCE  Of  HALLIFT)  24.009 
HCICHT  OF  CHARGE  ABOVE  FLOCR  SLAStfTI  3.IIB 


AREA  OF  fACeaN**?!  121AM, Q 

HQRWCNTAL  OISTANCC  FROM  CENTROID  Cf  AREA  TO  CCtlN) 

VERTICAL  OISTANCC  FROM  CENTROID  OF  AREA  TO  CCfIN)  •A4»7 


PRESSURE  TIKE  CURVE 


PQThT  1 PRCSSUREIPSn  1424*3 

TINE  tSEC) 

.44032 

POINT  ? PRESSUtCtPSlI  *4 

TINE  ISEC) 

.40439 

TINE 

FORCE  X 

FORCE  T 

HOHEHT  AOOUT  t 

*SEC» 

1L0SI 

(L49T 

tiH-LBS) 

2 

.4IA33 

1.2S43CA44 

.4 

1.1113EA14 

3 

•480AQ 

l,t43AE«4l 

• 0 

1.44221*14 

A 

.94447 

1.1294EA44 

.4 

1.0tS2E*l4 

4 

.44113 

1.47A3EA44 

.4 

9.4A19f*49 

A 

.441A4 

1.4197EA44 

*4 

9.19112*19 

T 

.44147 

9.69AACA47 

• 4 

4.M47F*49 

% 

«44MA 

9.l4A0r«47 

«l 

4.1742f*49 

1 

.14224 

4.99T9E*42 

• 4 

7»A794e«4A 

11 

.I42A7 

«.Qll4rA47 

.9 

7.1I9AE*49 

U 

•0477A 

7.AAA4I«97 

.4 

4.A998E*49 

12 

.44341 

A.9l4lfA47 

.4 

A.2444E*49 

IS 

.44327 

4.37 IAf»47 

.4 

9.71421*09 

lA 

.I439A 

9.4292EA47 

.4 

9.227)1*19 

19 

.I43U 

9.2747EAI7 

.4 

A. 73731*49 

lA 

.49A44 

A.7J23EA47 

• 9 

A.2A|9F»|9 

IT 

.4IASA 

A.l49|f»47 

.4 

3.79491*49 

14 

.4IA41 

3.4S93EA47 

.4 

3.24411*49 

19 

.44A4V 

3.4929C«4? 

• 4 

2.77971*49 

24 

,41919 

7.9A4AF*4f 

.4 

2.21971*19 

21 

• 4I9A1 

t, 99992*47 

>9 

1.79All*|9 

22 

.14944 

|.A939C*47 

.4 

1. T|AA€*49 

23 

,14999 

9.474IEA4A 

.4 

t.l3941*44 

2A 

.44422 

3.44942*44 

• 4 

9.23979*14 

ARFI 

Qf  PaCtltN«»2l 

94199.9 

M04I2CNT4L 

Q19I4NCC  FROM  CfNTttOtO  CF  AREA  TO 

CCtlN)  >42.9 

VfATKAl  OSST4HCI  FRC«  CENTROID 

OF  4REA  TO  COIIHI  7.4 

MESSUtf  TINE  CURVE 


POINT  X PtVSSUNCIPSn  IPRA.A  HMfrSEtt  «MM4 

POTNI  ? PMS1URC(P$ll  .4 
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itNC  roflct  X FOR<:t  t honcni  about  z 

TSECI  UBM  CL6S)  tTN>LBSI 


1 

•filtOA  .1 

•A.3tS1E«l1 

z 

•tQosa  .1 

1«1?S6C»|7 

•7,cii?e*oi 

3 

«tmi  .0 

-7.4C?1E*C1 

4 

•.94l«E*07 

-7.0C6lt*C1 

S 

• cuts  *1 

•.CA?CE*I7 

-4.611tl*l1 

6 

•iotic  .« 

M1UC*C7 

•6.?S14E«C1 

T 

• cm?  .1 

7.  094?Ci|7 

•1»At17E«01 

1 

•tlt14  •« 

f‘.510SI*C7 

•1.4C49r»01 

•CIZZI  .1 

C.04C1E*C7 

•4.1IASE»I1 

11 

•cc?«r  .« 

i.i4?ce»i* 

•4,17?Ct*C1 

It 

»M??4  .C 

1.0SC7e»07 

•4.tnie*ii 

iZ 

•CCSIt  *1 

4»1S41E*C7 

•9»74l?r*0l 

U 

•cos??  *1 

4«0SICC*87 

-S.S?S1f*C1 

14 

•CCS14  .« 

S.1?7ie»l7 

•?.1|11C»I1 

11 

•CCSCl  *1 

S»l?3?f  17 

•?.M4?r»C1 

u 

• M4M  «C 

?.9t14r«Q7 

•?.|7I9C*I1 

ir 

.CC4S%  •• 

?.Btin*D? 

•i.ce?ic*ii 

ti 

•cmi  t 

t.9tUC*|7 

•l.?47tt*l1 

11 

•mti  .0 

1. 0077E»»7 

•C.  U91C*QI 

Zl 

•CCItf  .* 

1. CSI?C»I» 

•IC141  •? 

• •17ZSC-C7 

•C.74?tf>09 

THE  PRINTED  OUTPUT  CONTAINING  THE  LOAD  TIME  HISTORIES 
FOR  CHARGE  NUMBER  3 HAS  BEEN  OMITTED. 
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CXAHPLf  PRO&LCM  0.1  CANTXlCVCR  HAU  eARRICR  ON  COMPACT  CRAVEL 


LOAOINC 


TIME  fO»CF  t FORCE  Y 

fSCCI  uesi  (LBS) 


t 

• 0Q886 

2 

.08833 

3 

.10068 

% 

.00087 

5 

•08113 

6 

.00190 

7 

•00167 

• 

.00199 

9 

.00770 

tl 

.00797 

It 

.10779 

u 

.10301 

IS 

.00377 

tA 

.00389 

19 

.00381 

16 

• 00908 

17 

.00939 

18 

.00961 

19 

.00980 

2Q 

.08819 

71 

.81891 

77 

.80868 

73 

.00598 

79 

.8067? 

79 

•Q869« 

76 

•00678 

77 

.80707 

78 

.00779 

79 

.00788 

38 

•00767 

3.A122F*I)A 

S.2»Q<)E»0A 

3.1A^SF»0a 

3.oift?f*oa 

e.iAFAF^aa 

?.7SS5F*aa 

7.s?Atc*oa 

?.A927F»aa 

2»73iee«ea 

2«0$S7f>0a 

t.5673F»BI 

t.a3sic*ei 

t.70Aef4ia 

i.s733F*ea 

t.AA14E»ia 

t.siosf*aa 

ia7<i?r*aa 

t.OAr^FABI 

^•l*S?FM7 

7*«irr*07 

«>s3t7f4ar 

^•7?A6f»07 

^•8i7Rr*aT 

3.3300C*07 

?,S«?^F»I7 

t,7$Mf*07 

t.07a«F*D7 

7.«17AF*9I 


?.3H<)e*oa 
7.29A3E«0a 
Z* 1917E*0a 
?.Qa<)iCfoa 
i.9assc*oa 
i.aa3<)F*Da 
l,7U3E*0a 
1.678$E*oa 
l»$760F«Oa 
l.A73AC»&a 
u37«aE«Qa 
U2i»7E»0a 

1.0S29E*0t 

9*b033C»l7 

«*S771E*07 

7,SI«9€»07 

6.S?AIE»07 

S.A9nE«87 

%.A7?5E»07 

3.A4S3E»87 

7*97%OE*07 

e.A0i7F*a7 

l.a79AE;D7 

|.3S7U*«7 

••3AaAf*0S 
3.  US7E*0S 

.fi 

• I 

.0 


ro«CF  vccTON  rouAiS  ictc  apter  time  ,ao7L 


MOMENT  ABOUT  2 
(1N*IB$> 

-U9775EM0 

1.1S9AF»1D 

1.13S7E»10 

1.1030E»iO 

i«069ac»io 

1.036SE«1Q 

1.Q033E«1Q 

9.7009E»I9 

9*3Sa7E»09 

9*03S5f4D9 

«.70A3E*09 

a.37?8E»Q9 

a.039IE^09 

7.707SE*0S 

7.375AE*89 

7.DA3U*09 

S.7t89C*l9 

S*37an«09 

B.DA69E»89 

5.71A?E*09 

9.3I78E*89 

A.STA1F»D9 

3.8A67E»09 

3.UA3F*09 

7.M7?E*89 

7.7997F*09 

?»8A7?C»89 

l.bU7F*89 

9.?378F*I8 

?.3A89<*88 


NAiL  ■CACWtS  INCtPUNf  FAliURE  TIME 


87«l 


eiAMPU  PflOeieK  o.t  CtNTtLevCR  uau  sarsur  on  compact  gnavcl 


OCttNT  AA.OO 

SCUfOflCgMCNt  RATIO 

Ife^i  OYNAhtC  OCStCN  STRCSSfpSI)  RODOg*0 

ULTlflTE  Resist ANCCTPSn 

TCTAl  IMPULSC  OK  HAU  IPSt-KSi  9900.0 

ftKf  Af  NMICM  UAt.t.  PAUStSfCI 


o««aaaa»«*««o«oao«o««««aoi»e*»»»*««»«*Baooa*o«oM 

UtliUJliiUjUJWUJUJWIUUtUUiWMJWWWmMltflUWUlItJtUWtitUJIUhiluWUiWlUWWUUUJUUjUiUlU 
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O'  N • K f 
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<tMa#««KrtfNeoa*i^40 
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JUJhlWUUWU 
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LlWUlUJuJWUtWh 


> a a o a c 
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J lu  W U Uj 
^ O K 

0 o ■#  • « 

kH  «« a 
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Figure  O.S  Input  data  sheet  - Example  0.2:  Card  Types  I and  2 


Card  Type  11 


Charge 


Input  data  sheet  - Example  D.2:  Charge  Number  1,  Card  Type  12  for  roof  and 

front  wal"! . 


iput  data  sheet  - Example  D.2:  Charge  Number  2,  Card  Type  12  for  backwall 

d floor. 


Figure  0.14  Input  data  sheet  - Example  0.2:  Charge  Number  2,  Card  Type  12  for  roof  and 

front  v/al  1 . 


Figure  D.15  Input  data  sheet  - Example  D.2:  Charge  Number  3,  Card  Type  12  for 

backv/all  and  floor. 


Input  data  sheet  - Example  D.2:  Charge  Number  3,  Card  Type  12  for  roof  and 

front  v;all . 


Number 


Charne  Number  4,  Card  Type  12  for  roof  and 


eet  - Example  D.2:  Charge  Number 


Input  data  sheet  - Example  0.2:  Charge  Number  5,  Card  Type  12  for  roof  and 

front  v/al  1 . 


Input  data  sheet  - Example  D.2:  Charge  Number 

backv/all  and  floor. 


[nput  data  sheet  - Example  D.2:  Charge  Number  6,  Card  Type  12  for 

'oof  and  front  v/all. 
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P408LEH  SOCUTICK  COHPCtlEO 


APPENDIX  E 


LIST  OF  SYMBOLS  USED  IN  TEXT 

A Surface  area  of  foundation  (ft^) 

Cross-sectional  area  of  steel  reinforcing  bar  (in^) 

A,^in  Minimum  area  of  flexural  reinforcement  (in^/ft) 

Ar  Area  of  tension  reinforcement  within  a width  "b" 

(in2/ft) 

a Depth  of  equivalent  rectangular  stress  block  (in) 

B Width  of  foundation  along  axis  of  rotation  for 

rocking  or  normal  to  direction  of  horizontal 
force  (in) 

b Width  of  compression  face  of  flexural  member  (in) 

c Thickness  of  bottom  concrete  cover  (in) 

Cr,  c..  Distance  from  resultant  of  applied  loads  to  axis  of 

^ rotation  for  Sectors  I and  II,  respectively  (in) 

d Distance  from  extreme  compression  fiber  to  centroid 

of  tension  reinforcement  (in) 

d[j  Diameter  of  tension  reinforcement  bar  (in) 

d Distance  between  centroids  of  the  compression 

and  tension  reinforcement  (in) 

f^  Friction  factor  for  cohesive  soils 

f'  Static  ultimate  compressive  strength  of  concrete 

at  28  days  (psi) 

Dynamic  design  stress  for  reinforcement  (psi) 
fg  Static  design  stress  for  reinforcement  (psi) 

G Shear  modulus  for  soil  (psi) 

g Acceleration  of  gravity  (in/sec^) 
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Distance,  in  vertical  direction,  between  supports 
and/or  free  edges  of  foundation  extension  supported 
on  three  or  four  sides  (in) 

Resultant  of  horizontal  blast  loads  acting  on 
structure  at  time  "t"  (lbs) 

Height  of  backwall  (ft) 

Charge  location  parameter  (ft) 

Distance  from  center  of  gravity  of  structure 
to  soil -structure  interface  (in) 

Mass  moment  of  inertia  of  structure  (Ib-sec^/in) 

Scaled  unit  blast  impulse  (psi-ms/lb^/2) 

Unit  blast  impulse  (psi-nis) 

Total  spring  constant  for  soil  for  horizontal 
translation  (lbs/1n) 

Total  spring  constant  for  soil  for  vertical 
translation  (Ibs/in) 

Spring  constant  for  soil  element  for  horizontal 
translation  (Ibs/in) 

Spring  constant  for  soil  element  for  vertical 
translation  (Ibs/in) 

1)  Length  of  rectangular  foundation,  in  plane  of 
rotation  for  rocking  or  in  direction  of  horizontal 
force  (in) 

2)  Length,  in  horizontal  direction,  of  foundation 
between  supports  and/or  free  edges  (in  or  ft) 

1)  Length  of  loaded  area  of  foundation  (ft) 

2)  Length  of  foundation  extension  of  cantilever 
wall  barrier  (in) 

Charge  location  parameter  (ft) 

Distance  from  face  of  support  to  critical  section 
for  shear  for  simple  type  foundation  extension  (in) 


Clear  span  to  face  of  support  for  simple  type 
foundation  extension  (in) 

M Applied  unit  design  load  moment  at  critical  section 

for  shear  for  thick  foundation  (in-lbs/in) 

Mp^  Unit  bending  moment  at  face  of  support  for  foundation 

extension  of  single  cell  barrier  (in-lbs/in) 

Mp^  Ultimate  unit  bending  moment  capacity  of  foundation 

required  to  develop  blast  wall  (in-lbs/in) 

Mum  Ultimate  unit  negative  moment  capacity  in  horizontal 

direction  (in-lbs/in) 

M^p  Ultimate  unit  positive  moment  capacity  in  horizontal 

direction  (in-lbs/in) 

M(t)  Moment  of  resultant  of  blast  loads  about  the  z axis 

at  the  center  of  gravity  of  the  structure  at  time  "t" 
(in-lbs) 

M^J  1)  Ultimate  unit  resisting  moment  (in-lbs/in) 

2)  Applied  unit  design  load  moment  at  a section 
(in-lbs/in) 

Ultimate  unit  negative  moment  capacity  in  vertical 
direction  (in-lbs/in) 

Myp  Ultimate  unit  positive  moment  capacity  in  vertical 

direction  (in-lbs/in) 

Mu,,  Ultimate  unit  moment  capacity  of  backwall  element 

^ (in-lbs/in) 

m Mass  of  structure  (Ibs-sec^/in) 

N Blow  count  from  standard  penetration  test 

NS  Number  of  soil  elements  used  in  overturning  analysis 

Pg  Soil  bearing  pressure  at  face  of  support  of  simple 

type  foundation  extension  (psi) 

P^^  Soil  bearing  pressure  at  critical  section  for  shear 

for  foundation  extension  (psi) 
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Reinforcement  ratio  in  horizontal  direction  on 
each  face 

Py  Reinforcement  ratio  in  vertical  direction  on 

^ each  face 

Pjj  Ultimate  unit  internal  resistance  of  foundation 

extension  supported  on  three  or  four  sides  and 
subjected  to  a trapezoidal  loading  (psi) 

p^  Reinforcement  ratio  equal  to  A^/bd 

Ry^  Normal  distance  from  charge  to  backwall  (ft) 

R^  Horizontal  resistance  of  soil  (lbs) 

Ry  Vertical  resistance  of  soil  (lbs) 

Rq  Moment  of  horizontal  and  vertical  soil  resistance 

forces  about  the  center  of  gravity  of  structure 
(in-lbs) 

Rj,  Rjj  Total  internal  resistance  of  Sectors  I and  II, 

^ respectively  (lbs) 

Thickness  of  concrete  section  (in) 

TS  Thickness  of  foundation  slab  (in) 

TW  Thickness  of  backwall  (in) 

t Time  (sec) 


ty^  Arrival  time  of  blast  wave  (sec  or  ms) 

tg  Duration  of  positive  phase  of  blast  pressure  (ms) 

u Horizontal  displacement  of  structure  (in) 

li  Horizontal  acceleration  of  structure  (in/sec^) 


Static  deflection  under  weight  of  structure  at  soil 
element  actachment  point  (in) 


V(t)  Resultant  of  vertical  blast  loads  on  structure  at 

time  ''t"  (lbs) 


Design  load  for  foundation  extension  supported  on 
three  or  four  sides  (lbs/ in) 

V^J  Total  applied  design  shear  force  at  critical 

section  (lbs) 

V Vertical  displacement  of  structure  (in) 

V Vertical  acceleration  of  structure  (in/sec^) 

Nominal  permissible  shear  stress  for  concrete  (psi) 

Vu  Shear  stress  at  critical  section  for  shear, 

Sector  II  (psi) 

Vw  Shear  stress  at  critical  section  for  shear. 

Sector  I (psi) 

W Charge  weight  (lbs) 

X Yield  line  location  in  horizontal  direction  (in) 

Horizontal  distance  from  center  of  gravity  of 
structure  to  soil  element  attachment  point  (in) 

y Yield  line  location  in  vertical  direction  (in) 

Scaled^ formal  distance  from  charge  to  backwall 

Zc  1)  Scale  nornal  distance  from  charge  to  foundation 

slab  (ft/lbl/3) 

2)  Minimum  scaled  distance  from  charge  to 
foundation  (ft/lbV3) 

6 Influence  coefficients  for  horizontal  and  vertical 

^ spring  constants  for  soil  elements 

a Distance  between  soil  elements  on  foundation  (in) 

0 Rotation  of  structure  (deg) 

d Angular  acceleration  of  structure  (rad/sec^) 

p Poisson's  ratio  for  soil 
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a Vertical  stress  in  soil  (psi) 

Overturning  angle  (deg) 

<!i  1}  Capacity  reduction  factor 

2)  Bar  diameter  (in) 
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LIST  OF  SYMBOLS  USED  IN  INPUT  TO  COMPUTER  PROHRAM 

AFACE  Area  of  loaded  surface  on  structure  (in?) 

B Length  of  foundation,  in  plane  of  rotation  for 

rocking,  or  in  direction  of  horizontal  force  (in) 

C9  Ratio  of  foundation  thickness  to  backwall  thickness 

Distance  betoeen  the  centroids  of  the  compression 
and  tension  reinforcement  (in) 

El.  E2  Moduli  of  elasticity  of  bilinear  stress-strain  curve 
for  soil  (psi) 

f(.  1)  Friction  factor  for  cohesive  soils 

2)  Adhesion  constant  for  non-cohesive  soils  (psf) 

fjls  Dynamic  design  stress  for  reinforcement  (psi) 

HAUNH  Width  of  backvMall  haunch  (in) 

HB  Ratio  of  height  of  backwall  to  length  of  foundation 

HW  Height  of  backwall  element  (ft) 

h Charge  location  parameter  (ft) 

I Mass  moment  of  inertia  of  structure  (Ib-sec^/in) 

ICl  Computer  program  option  paraireter 

ICAl  Computer  program  option  parameter 

i^  Unit  blast  impulse  (psi-ms) 

(KLM)u  Load-mass  factor  in  the  ultimate  range 

L Length  of  backwall  element  (ft) 

1 Charge  location  parameter  (ft) 

N Number  of  blast  walls  in  structure 

NFON  Computer  program  option  parameter 

NDEL2  Constant  for  changing  integration  time  step  in 

overturning  analysis  computer  program 


1)1 


OAD 


Computer  program  option  parameter 
Number  of  charges 


& 

r- 


NP 

NUMPT 

NUMTM 

NVEL 

NS 

NWALL 

P2/P1 

Ph 

PV 

Ra 

SMTP 

SSR 

TASM 

TOLG 

"AO 

UVECT(i) 

Vf 


Number  of  integration  time  steps  skipped  between 
output  time  stations 

Number  of  integration  time  steps  used  in  over- 
turning analysis 

Computer  program  option  parameter 

Number  of  soil  elements  used  in  overturning  analysis 

Computer  program  option  parameter 

Ratios  of  initial  pressure  to  final  pressure  on 
surface 

Reinforcement  ratio  in  horizontal  direction 

Reinforcement  ratio  in  vertical  direction 

Normal  distance  from  charge  to  backwall  (ft) 

Stress  at  which  modulus  of  bilinear  soil  stress- 
strain  curve  changes  (psi) 

Ratio  of  length  of  loaded  area  of  foundation  to 
total  length  of  foundation 

Time  of  arrival  of  blast  wave  on  structure  (sec) 

Smallest  duration  of  loading  produced  by  any  one 
charge  on  any  surface  of  the  structure  (sec) 

Thickness  of  backwall  (in) 

Arrival  time  of  blast  wave  (sec  or  ms) 

Arrival  time  plus  duration  of  positive  phase  of 
blast  pressure  (sec) 

fieri zontal  or  vertical  com{>onent  of  unit  vector 
nomwl  to  a loaded  surface  of  structure 

Velocity  of  post  failure  fragments  (in/ms) 
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w 

WT 

XB 

XCORD(i) 

XR 

X 

YB 

y 

6x» 

u 

db 


Charge  weight  (lbs) 

Weight  of  structure  (lbs) 

Horizontal  distance  from  center  of  gravity  of 
structure  to  rear  face  of  backwall  element  (in) 

Horizontal  distance  from  center  of  gravity  of 
structure  to  centroid  of  loaded  surface  (in) 

Horizontal  distance  from  center  of  gravity  of 
structure  to  left  end  of  foundation  (in) 

Yield  line  location  in  horizontal  direction  (in) 

Vertical  distance  f’^om  center  of  gravity  of 
structure  to  top  of  foundation  slab  (in) 

Yield  line  location  in  vertical  direction  (in) 

Influence  coefficients  for  horizontal  and  vertical 
spring  constants  for  soil  elenvents 

Poisson's  ratio  for  soil 

Summation  of  unit  blast  impulses  on  backwall 
element  (psi-ms) 


r 
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